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SUMMARY 


This  report  documents  the  development  and 
findings  of  a  computer  model  that  simulates  the 
behavior  of  the  Eaton-Thioat  Valve  Element 
(ETVE)  prototype,  and  describes  a  proposed 
modification  concept  for  the  ETVE.  This  high¬ 
speed  valve  was  designed  for  the  U.  S.  Army 
Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  to  control  the  simulation  of  the 
dynamic  effects  of  a  nuclear  blast 

The  computer  model  reveals  three  main  find¬ 
ings:  (1)  the  ETVE  chatters  during  the  open 
cycle,  (2)  the  chatter  is  caused  by  the  high  gas 
forces  on  the  diding  sleeve  as  the  driver  gas 
passes  through  its  portholes,  and  (3)  the  chatter  is 
aggravated  because  there  is  insufficient  damping 
in  the  system 

The  ETVE  simulation  model  has  proven  to 
be  an  extremely  valuable  tool  in  assessing  the 
qualitative  nature  of  the  current  valve’s  operation 
and  is  indispensable  in  assessing  the  effects  of 
valve  modificatiorrs  or  redesigns.  The  ETVE’s 
complex  nonlinear  behavior  is  impossible  to  pre¬ 
dict  in  any  way  other  than  numerical  simulation, 
even  if  the  numerical  simulations  are  accepted  on 
a  qualitative  basis  only. 

In  addition,  the  model  can  be  used  to  input 
new  Belleville  washer  characterization  data  for 
different  washer  configurations  and  then  assess 
the  qualitative  effects  of  each  configuration  on 
the  overall  valve  behavior. 

The  INEL  recommends  opening  the  ETVE 
by  sliding  the  sleeve  toward  the  downstream  end 
of  the  valve  instead  of  toward  the  upstream  end. 


as  the  ETVE  is  presently  configured,  and  to  pro¬ 
vide  additional  damping  to  the  system.  The  main 

areas  of  redesign  follow. 

1.  Increase  the  length  of  the  sleeve  by  approxi¬ 
mately  1 .6  in.  and  add  two  sets  of  seals  in  the 
added  length  to  seal  the  last  row  of  body 
portholes.  Reduce  the  length  of  the  liner  by 
0.7  in. 

2.  Place  the  sleeve  upstream  (from  its  current 
as-built  location)  as  indicated  in  Figure  5-1 
so  the  pneumatics  can  slide  the  sleeve  for¬ 
ward  to  its  open  position. 

3.  Increase  the  length  of  the  ETVE  body  to 

house  approximately  32  Belleville  washers 
on  the  open  chamber  and  8  washers  in  the 
close  chamber.  This  modification  requires 
retapping  the  pneumatic  actuation  lines. 

4.  Lengthen  the  piston  rod  to  accommodate  the 
added  number  of  washers. 

5.  Increase  the  stop  plate  size  and  increase  the 
number  and  size  of  bolts  used  to  fasten  the 
stop  plate  to  the  valve  body. 

3.  The  close-cycle  actuation  needs  to  be  evalu¬ 
ated  and  redesigned  to  overcome  the  driver 
gas  forces  that  will  make  the  valve  difficult 
to  close. 

4.  Any  general  design  concepts  should  be  eval¬ 
uated  by  using  the  already  developed  valve 
sleeve  dynamic  simulation  model  to  qualita¬ 
tively  predict  and  evaluate  the  design  con¬ 
cept’s  effects  on  the  ETVE  performance. 


Intentionally  left  blank. 
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1.  INTRODUCTION 


1.1  Background 

As  part  of  an  ongoing  effort  to  improve 
techniques  for  simulating  nuclear  blasts,  the 
U.S.  Army  Research  Laboratory  (ARL)  has 
been  studying  the  merits  of  computer-controUed 
valves.  The  valve  studies  have  been  aimed  at 
providing  the  Department  of  Defense  (DOD) 
with  the  capacity  to  conduct  nuclear  blast  and 
thermal  survivability  testing  on  full-scale  tacti¬ 
cal  vehicles.  In  support  of  this  DOD  objective, 
two  new  simulators  are  currently  under  devel¬ 
opment:  (1)  the  Large  Blast/Thermal  Simulator 
(LB/TS),  a  Defense  Nuclear  Agency  facility 
currently  under  construction  at  the  White  Sands 
Missile  Range,  and  (2)  the  ARL  1/6-scale  Test 
bed,  a  technology  demonstrator  currently  under¬ 
going  characterization  testing  at  the  Aberdeen 
Proving  Ground,  Maryland.  If  use  of  computer 
controlled-valves  proves  feasible  and  cost 
effective,  it  is  expected  that  the  computer-con- 
trolled  valves  will  replace  the  more  labor-inten¬ 
sive  rupture  disks  currently  used  to  simulate 
nuclear  blasts.  The  computer-controlled  valves 
would  be  retrofitted  first  to  the  ARL  Test  Bed 
and  then  to  the  LB/TS  at  the  White  Sands  Mis¬ 
sile  Range. 

In  support  of  the  U.S.  Army’s  woik,  the 
Idaho  National  Engineering  Laboratory  (INEL) 
studied  a  number  of  valve  design  concepts  and 
determined  that  a  multi-element  valve  would 
best  control  blast  simulation.  A  multi-element 
valve  is  formed  from  an  array  of  small  valves 
whose  openings  and  closings  are  coordinated  to 
give  the  required  total  output  through  the  valve 
as  a  whole.  The  INEL  subcontracted  Eaton  Con¬ 
solidated  Controls  to  btiUd  the  valve  element 
prototype,  called  the  Eaton  Throat- Valve  Ele¬ 
ment  (ETVE),  which  uses  a  pneumatically-dri¬ 
ven  sliding-sleeve  (see  Figure  1-1).  The  term 
throat  was  given  to  the  prototype  since  the  valve 
acts  like  a  throat  in  that  the  driver-gas  released  is 
funneled  through  an  opening  much  narrower 


than  the  driver  cross  section.  We  tested  this  pro¬ 
totype  ETVE  at  the  INEL  during  1991.  The  test 
results  revealed  segmented  shock  waves,  which 
suggested  the  possibility  that  the  ETVE  chat¬ 
tered  during  the  ETVE  opening  cycle.  We  also 
observed  that  the  duration  of  the  chatter 
increased  with  increased  gas  pressure  in  the 
upstream  pressure  chamber  (hereafter  called  dri¬ 
ver  tube). 

The  primary  objective  of  our  work  per¬ 
formed  during  1993  has  been  to  develop  a  com¬ 
putational  fluid  dynamics  model  of  the  valve  and 
then  use  it  to  propose  a  preliminary  modification 
design  concept  of  the  ETVE  to  meet  the  original 
design  requirements.  If  the  U.S.  Army  approves 
the  design,  and  funds  are  available,  the  INEL 
may  perform  the  final  analysis,  design,  modifica¬ 
tion,  and  testing  of  the  ETVE,  and  then  ship  the 
valve  back  to  ARL  to  perform  full  conditions 
testing  there  (700°F  and  1800  psi  nitrogen  driver 
gas). 

1.2  Valve  Testing  History  at  the  INEL 

The  INEL  tested  the  ETVE’s  performance  in 
two  separate  phases:  the  first  phase  during  1989, 
the  second  during  1991.  The  1989  phase  used 
the  ETVE;  however,  the  effort  was  focused  on 
the  equipment  that  had  been  purchased  from 
Thermal  Sciences  Incorporated  of  St.  Louis, 
Missouri.  The  equipment  consisted  of  the  high- 
pressure  tube  (called  the  driver  tube),  its  interior 
insulation,  and  the  gas  heater.  ^ 

The  1991  phase  focused  on  testing  the  per¬ 
formance  of  the  ETVE,  during  which  testing  the 
chattering  condition  was  diagnosed  as  a  possible 
malfunctioa^  The  main  objection  to  the  ETVE’s 
chatter  is  its  effect  on  the  downstream  shock  sig¬ 
nal;  as  the  valve  chatters,  the  effective  flow  area 
for  the  driver  gas  changes,  causing  a  segmented 
wave  downstream,  which  is  imacceptable  for 
blast  simulation. 
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Figuro  1  -1 .  Photograph  of  ETVE  main  components. 


1.3  Current  ETVE  Configuration 

Figure  1-2  illustrates  the  current  ETVE  con¬ 
figuration  .  The  main  features  of  the  current  con¬ 
figuration  are  as  follows: 

•  The  ETVE  is  designed  with  three  rows  of 
portholes  in  both  the  ETVE  body  and  sleeve. 
When  the  corresponding  portholes  are 
aligned,  the  gas  in  the  upstream  driver  tube 
flows  freely  through  these  portholes. 

•  The  ETVE  is  closed  when  the  sliding  sleeve 
is  in  its  most  downstream  position,  as  indi¬ 
cated  in  Figure  1-2. 

•  The  ETVE  is  actuated  back  to  its  open  posi¬ 
tion  by  an  electronic  signal  sent  to  the  open- 
solenoid,  which  opens  a  2000-psi  gas  supply 
line.  The  gas  pressurizes  the  open-actuating 
line  and  the  open-chamber.  As  the  chamber 
is  pressurized,  the  piston,  which  is  rigidly 
connected  to  the  sliding  sleeve,  will  move  to 
align  the  body  and  sleeve  portholes.  This  is 
the  ETVE  open  position,  illustrated  in  Fig¬ 
ure  1-3. 

•  The  ETVE  is  actuated  to  its  closed  position 
by  an  electronic  signal  that  actuates  both  the 
close-solenoid  to  supply  2000  psi  to  the 


close-chamber,  and  the  open-solenoid  to 
relieve  the  pressure  in  the  open-chamber  and 
its  actuating  line. 

1.4  Purpose 

The  primary  objective  of  the  work  per- 
fomaed  during  1993  has  been  to  develop  a  com¬ 
putational  fluid  dynamics  model  of  the  ETVE 
that  can  be  used  to  propose  a  preliminary  modi¬ 
fication  design  concept  to  meet  the  ETVE’s  orig¬ 
inal  design  requirements.  The  gas  dynamics 
model  uses  the  ETVE  moving  parts  equation  of 
motion  and  associated  parameters  calculated  in 
the  mechanical  dynamics  subroutine.  This  sub¬ 
routine  calculates  these  parameters  by  using  the 
results  from  tests  we  performed  on  the  ETVE  to 
characterize  the  sleeve’s  displacement  history, 
and  from  tests  we  performed  on  Belleville  wash¬ 
er  stacks  similar  to  those  in  the  ETVE.  Note  that 
the  characterization  of  the  ETVE  sleeve  dis¬ 
placement  was  performed  under  ambient  pres¬ 
sure  and  temperature  conditions  on  the  upstream 
end  of  the  ETVE.  This  report  presents  the  results 
from  the  characterization  tests  for  the  sleeve  and 
the  BeUevUle  washers,  and  the  results  from  the 
mechanical  dynamics  subroutine  and  the  gas 
dynamics  model. 
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Figure  1-2.  Current  ETVE  configuration— ETVE  closed. 
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Current  ETVE  configuration — ETEV  opened. 
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2.  VALVE  COMPONENT 
CHARACTERIZATION  AND  VALVE  MODELING 


We  developed  the  computer  model  in  four 
phases:  (1)  characterization  of  the  sleeve’s  dis¬ 
placement  history,  (2)  characterization  of  the 
Belleville  washer  stack  behavior  under  quasi-sta¬ 
tic  loading,  (3)  development  of  the  mechanical 
dynamics  subroutine,  and  (4)  development  of  the 
gas  dynamics  model.  The  sleeve  displacement 
and  BeUeviUe  washer  characterization  tests  were 
necessary  to  complete  the  mechanical  dynamics 
subroutine.  In  turn,  the  mechanical  dynamics 
subroutine  was  necessary  prior  to  finishing  the 
gas  dynamics  model.  Each  of  the  four  phases  are 
described  in  this  section. 

2.1  ETVE  Characterization  Testing  for 
the  Existing  Configuration 

We  performed  two  sets  of  tests  to  character¬ 
ize  the  behavior  of  the  ETVE  under  its  existing 
configuration.  The  first  test  characterized  the 
ETVE  sleeve’s  displacement  behavior  using  a 
linear  variable  displacement  transducer  (LVDT). 
The  LVDT  data  were  used  in  the  development  of 
the  ETVE  moving  parts  mechanical  dynamics 
subroutine.  The  second  test  set  explored  the 
effects  of  varying  the  open-  and  close-actuating 
solenoids  sequencing  to  reduce  or  eliminate  the 
ETVE  chatter.  Appendix  A  describes  this  work 
in  greater  detail. 

2.1.1  Characterization  Test  1  —  LVDT 
Testing.  The  tests  performed  at  the  INEL  dur¬ 
ing  1991  indicated  that  the  ETVE  was  possibly 
chattering  during  the  opening  cycle  before 
achieving  steady-state  open  positioa  At  that 
time,  when  the  ETVE  arrived  at  the  INEL  for  its 
initial  testing  in  1989,^  it  was  not  possible  to 
exactly  confirm  the  chatter,  since  the  only  means 
provided  to  detect  the  sleeve’s  movement  was  a 
Hall  Effect  Sensor.  This  instrument  provided 
only  qualitative  displacement  information.  Prior 


to  characterizing  the  sleeve’s  displacement  histo¬ 
ry,  we  purchased  and  installed  an  LVDT  and  used 
it  to  obtain  a  high-resolution  sampling  of  the 
ETVE’s  sliding  sleeve  displacement  during  the 
open  cycle.  It  is  important  to  note  that  we  tested 
the  ETVE  at  ambient  temperature  and  pressure 
because  the  driver  tube  was  not  available  for  use 
at  the  INEL  facility  during  this  1993  work. 

2.1.2  Characterization  Test  2  — 
Damping  with  Soienoid  Vaives.  We  per¬ 
formed  these  tests  to  explore  the  possibility  of 
reducing  or  eliminating  the  chatter  by  varying 
the  sequencing  of  the  acmating-solenoids.  We 
tried  two  different  actuating  sequences: 

1.  We  used  2000  psi  to  open  the  ETVE  while 
the  close-chamber  was  pressurized  to  1000 
psi;  we  then  vented  the  close-chamber  70 
msec  after  the  LVDT  detected  any  move¬ 
ment  (see  Appendix  A  for  detailed  test  infor¬ 
mation). 

2.  Similar  to  Sequence  1,  we  pressurized  the 
close-chamber  to  1000  psi  and  used  2000  psi 
to  actuate  the  ETVE  to  open.  However,  we 
kept  the  close-chamber  pressurized  at  1000 
psi  for  500  msec  after  the  LVDT  detected 
any  movement  before  venting  it  Figure  2-1 
shows  overlaid  plots  of  the  ETVE  sleeve  dis¬ 
placement  for  the  tests  results  in  the  two 
sequences.  The  initial  LVDT  test  results  are 
described  in  Section  2.1.1. 

These  tests  show  that  the  actuating 
sequences  described  in  sequences  1  and  2  did  not 
change  the  chattering  considerably.  Further¬ 
more,  the  overlaid  plots  clearly  show  that  pres¬ 
surizing  the  close-chamber  before  opening  the 
valve  increases  the  valve  open  time  considerably, 
approximately  60  msec. 
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Figure  2-1 .  Overlaid  LVDT  test  plots. 


2.2  Belleville  Washer 
Characterization  Testing 

We  characterized  the  Belleville  washers 
under  quasi-static  load-deflection  behavior  for 
the  current  washer  configuration  and  for  other 
alternative  configurations.  We  tested  the  wash¬ 
ers  in  compression  and  monitored  the  deflection 
with  respect  to  the  load  applied  for  each  washer 
configuration. 

The  tests  reveal  tiiat  the  Belleville  washers 
experience  a  significantly  different  loading  and 
unloading  force-displacement  behavior,  and  that 
damping  decreases  exponentially  as  the  loading 
approaches  the  design  load.  To  provide  high- 
load  damping,  we  believe  it  is  necessary  to  stack 
several  washers  (eight  to  ten)  in  parallel  and  con¬ 
nect  them  in  series  with  one  or  more  similar 
stacks.  Appendix  B  describes  in  detail  the  test 
results  for  the  different  washer  stack  configura¬ 
tions. 

hi  addition  to  characterizing  the  washers,  we 
developed  a  simple  FORTRAN  program  using 
the  mechanical  dynamics  subroutine  to  study  the 
washer’s  damping  behavior.  We  observe  that 
overloading  the  washers  reduces  the  effective 
damping  significantly;  however,  the  Belleville 
washers  are  an  effective  method  to  provide 
damping  for  the  system  if  designed  correctly. 
Appendix  C  describes  the  FORTRAN  program 
in  greater  detail  and  discusses  the  conclusions 
drawn  from  its  use. 

2.3  Mechanical  Dynamics  Model 

We  analyzed  the  ETVE’s  moving  parts 
(sleeve,  piston,  piston  rod,  and  Belleville  wash¬ 
ers)  and  derived  die  equation  of  motion  for  their 
configuration.  The  resulting  solution  of  the  non¬ 
linear  equation  of  motion  is  represented  as 
explicit  time-difference  equations.  Appendix  D 
details  the  FORTRAN  implementation  of  the 
motion  solution  and  its  supporting  subroutines; 
namely  “belspr,”  “fact,”  and  “motion.” 


2.3.1  General  Code.  Since  all  of  the 
moving  parts  of  the  ETVE  are  rigidly  connected, 
the  mechanical  dynamics  model  developed  uses 
a  single  degree-of-freedom  rigid-body  dynamic 
model,  represented  in  Figure  2-2.  The  gas 
dynamic  model  described  in  Section  2.4  calls  up 
the  mechanical  dynamics  model  as  a  subroutine. 
The  LVDT  test  results  provided  the  sleeve  dis¬ 
placement  information  used  in  the  subroutine. 
The  Belleville  damping  characteristics  were 
obtained  from  the  quasi-static  compression  tests 
described  in  Section  2.2. 

2.3.2  Description  of  Subroutines. 

Subroutine  belspr  —  This  subroutine  com¬ 
putes  the  Belleville  springs’  stiffness  from  the 
force-displacement  characterization  values  spec¬ 
ified  in  BLOCK  DATA.  The  subroutine  is  only 
called  from  the  main  program  once  (prior  to  first 
call  to  motion). 

Subroutine  fact  —  This  subroutine  calcu¬ 
lates  the  current  actuator  force  acting  on  the 
moving  ETVE  parts.  This  subroutine  must  be 
called  prior  to  each  call  of  motion. 

Subroutine  motion  —  This  subroutine  cal¬ 
culates  the  incremental  motion  of  the  ETVE 
movable  parts  and  outputs  their  displacement, 
velocity,  and  acceleration.  Arguments  passed  to 
the  subroutine  are  current  calculation  time,  cur¬ 
rent  time  increment,  net  fluid  dynamic  force  act¬ 
ing  on  the  moving  mass,  and  an  integer  flag.  The 
flag  indicates  whether  a  Belleville  spring  is  load¬ 
ing  or  unloading,  which  is  very  important  owing 
to  the  very  large  changes  in  the  motion  parame¬ 
ters  that  occur  in  the  transition  between  relative¬ 
ly  free  motion  (neither  Belleville  spring  is  com¬ 
pressed)  and  very  stiff  resistance  (either  of  the 
Belleville  springs  are  compressed).  The  remain¬ 
ing  section  of  the  subroutine  consists  of  logic 
blocks  wherein  the  aj^licable  motion  parameters 
are  specified,  and,  finally,  the  updated  motion 
variables  are  computed.  Appendix  D  presents 
this  subroutine  in  Table  D-3. 
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^cb 

Gso  =  Sleeve  opening  gap 
Gsc  =  Sleeve  closing  gap 

M  =  Mass  of  the  moving  parts 

Fes  =  Coulomb  damping  on  sleeve 

Feb  =  Coulomb  damping  on  the  Belleville  washers 

Cs  =  Stiffness  and  viscous  damping  on  the  sleeve 
Ks  =  Stiffness  and  viscous  damping  on  the  coil  spring 

Kbc  =  Closing  Belleville  washers  stiffness 
Kbo  =  Opening  Belleville  washers  stiffness 
Cb  =  Belleville  washer  viscous  damping 

Externally  applied  force  resulting  from  the  time-dependent  actuator  pressure 
and  hydrodynamic  gas  pressure  (time-  and  sleeve  position-dependent)— 
Fp(y,t) 


Figure  2-2.  Mechanical  dynamic  model. 
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2.3.3  Findings.  Once  the  mechanical 
dynamic  model  was  completed,  the  actual  sys¬ 
tem’s  damping  characteristics  were  obtained  by 
running  several  algorithm  predictions  and  vary¬ 
ing  different  parameters  until  the  axial  displace¬ 
ment  history  matched  the  actual  test  data 
obtained  from  the  tests  described  in  Section  2.1 
and  2.2.  We  found  that  the  following  damping 
parameters  produce  a  reasonable  fit  to  the  test 
data: 

Sleeve  Coulomb  damping  force  =  100  Ibf 
Sleeve  viscous  damping  ratio  =  0.0 
Belleville  washers  Coulomb  damping  force = 0.0 
Belleville  washers  viscous  damping  ratio  =  0.08. 

Figure  2-3  compares  the  subroutine-calculat¬ 
ed  ETVE  displacement  response  damping  data 
with  the  test  data  obtained  from  actuating  the 
ETVE  to  open  at  ambient  pressure  and  tempera¬ 
ture,  as  described  in  Section  2.1.1. 

2.4  Gas  Dynamics  Model 

Figure  2-4  shows  the  gas  dynamic  simula¬ 
tion  model  representing  the  ETVE  configuration. 
The  stationary  and  moving  valve  parts  are 
approximated  as  rectangular  solid  blocks;  the 
mesh  spacing  used  is  9x  =  9y  =  1  mm,  in  the 
axial  and  radial  direction,  respectively.  Con¬ 
stant  pressure  boundary  conditions  were  applied 
at  the  inlet  and  outlet  boundaries  of  the  flow 
model.  On  the  high-pressure  boundary 
(upstream  in  Rgure  1-2)  the  pressure  is  set  to 
1800  psi  and  700°F;  the  low-pressure  boundary 
(downstream  in  Figure  1-2)  is  set  to  14  psi  and 
ambient  temperature.  Appendix  E  presents  the 
gas  dynamics  model  work  in  detail. 

2.4.1  Description  of  the  Modei.  We 

used  the  compressible  Euler  equations  formulat¬ 
ed  for  two-dimensional,  axisymmetric  coordi¬ 
nate  systems  to  describe  the  gas  dynamics  of  the 
ETVE  during  the  actuating  open  cycle.  The  code 
uses  a  variable  thickness  or  depth  in  die  radial 
direction  for  each  mesh  cell,  named  A. 


In  two-dimensional  computations,  the  vol¬ 
ume  of  a  cell  of  width  3x  and  height  9y  is  A3x9y, 
giving  a  quasi-three-dimensional  effect.  The 
physical  assumption  imderlying  the  variable 
depth  description  is  that  the  flow  variables  are 
depth-averaged  in  the  depth  direction  to  give 
mean  values  with  variation  only  in  the  x  and  y 
dimensions,  or  equivalently,  that  the  flow  vari¬ 
ables  are  independent  of  displacement  in  the 
depth  direction. 

The  variable-depth  equations  are  used 
because  they  greatly  increase  the  geometrical 
flexibility  of  the  description  without  increasing 
its  complexity.  The  variable.  A,  may  be  used  to 
represent  flow  through  a  two-dimensional  duct 
of  gradually  varying  thickness.  The  use  of  zero 
values  of  A  in  selected  mesh  cells  provides  a 
convenient  means  of  including  internal  obstacles 
in  the  flow  region.  In  particular,  a  cell  with  A  = 
0  will  allow  no  flow  across  its  boundaries.  The 
stationary  obstacles  (valve  internal  parts)  are  rep¬ 
resented  with  this  technique.  Axisymmetric 
coordinates  are  generated  by  having  A  increase 
linearly  with  the  distance  from  the  axis.  Since 
the  ETVE  valve  is  approximately  axisymmetri- 
cal,  this  feature  is  also  used.  As  presently  coded, 
the  A-quantity  must  be  a  constant  in  time. 

Equations  and  Constitutive 
Relations  —  The  Euler  equations  describing 
the  dynamics  of  an  inviscid  gas  have  been  cast  in 
many  fonns.  The  dependent  variables  are  cho¬ 
sen  here  to  be  the  mass  density,  p,  the  gas  veloc¬ 
ities,  u  and  V,  and  the  internal  energy  per  unit  vol¬ 
ume,  E.  In  terms  of  the  chosen  dependent  vari¬ 
ables,  the  basic  two-dimensional  Euler  equations 
used  in  the  code  are  the  continuity,  the  momen¬ 
tum,  and  the  internal  energy  equations. 

In  these  equations,  die  independent  variables 
are  time,  t,  and  coordinates,  x  and  y.  The  gravita¬ 
tional  acceleration  terms,  the  heat  conduction  terms, 
and  tire  viscous  terms  have  been  neglected,  owing 
to  their  insignificance.  To  complete  these  equa¬ 
tions,  a  constitutive  relatirai,  which  relates  pressure, 
p,  and  internal  energy  per  unit  volume,  E,  is  defined. 
The  following  equation  is  based  on  the  assumption 
that  the  gas  obeys  the  perfect  gas  laws. 
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Axial  Displacement  (in) 


Comparison  of  actual  LVDT  test  data  with  the  Mechanical  Dynamics  Model 
output  for  ambient  pressure  and  temperature  driver  gas.  Valve  actuated  to 
open  position  with  2000  psi,  ambient  pressure  in  the  close-chamber. 


Figure  2-3.  Calculated  versus  actual  ETVE  displacement  history. 
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Figure  2-4.  Gas  dynamics  model. 


p  =  (Y-  1)E 

The  relation  y  is  the  ratio  of  specific  heats. 

Dynamical  Coupling  —  To  represent  the 
valve  opening  dynamics,  the  moving  parts  of  the 
valve  must  be  coupled  to  the  subsequent  fluid 
dynamics  of  the  gas  as  it  flows  through  the  valve. 
The  gas  surrounding  the  valve  moving  parts 
exerts  a  force  (the  dominant  force)  on  these  mov¬ 
ing  parts  to  produce  its  motion.  The  ETVE 
sleeve  motion,  in  turn,  influences  the  gas  direct¬ 
ly  with  its  motion  and  indirectly  by  opening  or 
closing  the  portholes  to  alter  the  gas  flow  dynam¬ 
ics.  The  ETVE  sleeve  motion  is  also  influenced 
by  the  dynamic  response  of  other  components 
and  phenomena  such  as  the  coil  spring,  the 
ETVE  sleeve  friction,  the  Belleville  spring  load¬ 
ing  and  unloading,  and  the  actuating  pneumatics. 
The  gas  dynamic  system  and  ETVE  sleeve  mov¬ 
ing  parts  dynamics  are  solved  together  in  a  fiiUy 
coupled  manner. 

2.4.2  Findings  from  the  Current  ETVE 
Configuration.  The  overall  ETVE  model 
clearly  shows  that  the  driver  gas  dynamic  force 
on  the  sleeve  portholes  dominates  and  is  respon¬ 
sible,  when  coupled  with  the  Belleville  washer 
force,  for  the  sleeve’s  oscillations  during  the 
open  cycle.  The  other  dynamical  forces,  such  as 
the  actuating  pneumatics,  the  coU  spring,  and  the 
coulomb  damping  forces  during  the  open  cycle 
have  little  effect  on  the  ETVE  sleeve  motion,  as 
is  seen  in  the  following  discussion.  Examination 
of  the  driver  gas  dynamic  force  on  the  sleeve 
shows  that  this  force  predominantly  acts  on  the 
left  side  of  the  portholes  or  toward  the  direction 
of  valve  closure  for  the  current  configuration. 


Figure  2-5  plots -the  predicted  displacement 
of  the  sleeve  with  respect  to  time.  The  plot  was 
generated  by  rurming  the  model  using  8% 
BeUeville  washer  viscous  damping.  Figures  2-6 
plots  the  predicted  displacement  of  the  sleeve 
using  16%  Belleville  washer  viscous  damping. 

The  model  predicts  oscillation  rates  similar 
to  those  observed  during  the  1991  tests  at  the 
INEL.  The  1991  tests  suggest  average  oscilla¬ 
tions  of  153  cycles  per  second.  The  model  pre¬ 
dicts  average  oscillations  of  approximately  180 
cycles  per  second  for  8%  viscous  damping,  and 
125  cycles  per  second  for  16%  viscous  damping. 
These  results  indicate  that  shifting  the  system’s 
viscous  damping  can  shift  the  system’s  frequen¬ 
cy  response. 

It  is  important  to  keep  in  mind  that  (1)  the 
model  is  an  idealized  situation  with  an  infinite 
supply  of  gas  at  1800  psi,  contrasted  to  the  actu¬ 
al  test  situation  wtere  the  drive  gas  pressure 
decreases  very  quickly,  and  (2)  the  model  uses 
driver  gas  at  700°F,  contrasted  to  the  actual  test 
situations  where  the  driver  gas  temperature  was 
ambient  for  that  time  of  the  year  at  the  INEL 
(approximately  30°F).  Appendix  E  contains 
plots  of  other  variables  such  as  the  sleeve  veloc¬ 
ity,  sleeve  acceleration,  and  gas  pressure  at  the 
pressure  plot  point  with  respect  to  time. 

Since  the  model  shows  that  the  driver  gas 
forces  dominate  the  open  cycle,  we  began 
exploring  different  ETVE  configurations  to  use 
the  driver  gas  forces  to  our  advantage,  or  to 
reduce  their  overall  effect  on  the  valve. 
Section  3.  describes  the  different  configurations 
explored. 
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•  Sleeve  viscous  damping  set  at  8% 

•  Belleville  washers  Coulomb  damping  set  at  100  Ibf 

•  Portholes  are  fully  aligned  at  the  1.12  inch  marie  on  the  displacement  axis 

•  Belleville  washers  are  fully  compressed  at  the  1.2  inch  marie  on  the  displacement  axis:  1.12  inch¬ 
es  to  alignment  of  downstream  edges  of  portholes,  plus  0.05  inches  portholes  width  mismatch 
(washers  begin  to  be  compressed  at  1.17  inches),  plus  0.03  inches  maximum  washer  deflexion 
(form  Appendix  B  data  plots). 

•  Finer  inlaid  plot  of  force  represents  the  net  force  on  the  moving  parts  (sleeve,  piston,  piston  rod, 
and  Belleville  washers)  acting  in  the  axial  direction.  A  positive  value  represents  force  in  the 
upstream  direction. 

Figure  2-5.  Eight  percent  Belleville  washer  viscous  damping-model  plot. 
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Sleeve  viscous  damping  set  at  16% 

Belleville  washers  Coulomb  damping  set  at  100  Ibf 
Portholes  are  fully  aligned  at  the  1.12  inch  marie  on  the  displacement  axis 
Belleville  washers  are  fully  compressed  at  the  1.2  inch  marie  on  the  displacement  axis:  1.12  inch¬ 
es  to  alignment  of  downstream  edges  of  portholes,  plus  0.05  inches  portholes  width  mismatch 
(washers  begin  to  be  compressed  at  1.17  inches),  plus  0.03  inches  maximum  washer  deflexion 
(form  Appendix  B  data  plots). 

Finer  inlaid  plot  of  force  represents  the  net  force  on  the  moving  parts  (sleeve,  piston,  piston  rod, 
and  Belleville  washers)  acting  in  the  axial  direction.  A  positive  value  represents  force  in  the 
upstream  direction. 


Figure  2-6.  Sixteen  percent  Belleville  washer  viscous  damping-model  plot 


16 


Gas  Force  (Ibf) 


3.  EVALUATION  OF  ALTERNATIVE  CONF^GURATIONS 


3.1  Model  Runs 

The  main  objective  of  reconfiguring  the 
valve  is  to  diminish  or  eliminate  the  valve  chat¬ 
ter  effects  on  the  shock  wave  produced  by  the 
opening  of  the  valve.  We  made  many  different 
configuration  runs;  however,  only  those  configu¬ 
rations  that  show  the  most  significant  results  are 
described  here. 

3.1.1  Beveled  Portholes.  The  first  con¬ 
figuration  change  attempted  was  to  bevel  the 
downstream  edge  of  each  porthole  in  the  sleeve. 
This  intuitive  approach  was  to  obtain  the  least 
obstruction  possible  to  the  gas  thereby  decreas¬ 
ing  the  net  force  by  the  gas  on  the  sleeve.  Figure 
3-1  shows  the  configuration  used  in  the  attempt. 
The  sleeve’s  response  frequency  changed  some¬ 
what,  but  the  sleeve  stiU  oscillated  violently  dur¬ 
ing  the  open  cycle  (see  Figure  3-2). 

3.1.2  Backward  Sleeve  Configuration. 

Since  beveling  the  portholes  had  no  net  effect  on 
reducing  the  driver  gas  forces  on  the  sleeve,  we 
changed  the  ETVE  configuration  so  that  the  ini¬ 
tial  sleeve  position  was  shifted  to  the  right 
(upstream),  with  tiie  sleeve  ports  on  the  ri^t  side 
of  the  stationary  ports.  This  configuration  revers¬ 
es  the  opening  motion  to  the  left  (downstream),  as 
illustrated  in  Figure  3-3.  We  term  this  configura¬ 
tion  backward.  The  intuitive  thinking  in  this  case 
was  that,  perhaps,  the  gas  dynamic  forces  on  the 
sleeve  would  assist  the  open  actuating  force  to 
thrust  the  valve  into  the  open  position. 

Running  the  model  with  the  sleeve  backward 
forced  us  to  lengthen  the  downstream  end  of  the 
sleeve,  for  two  reasons:  (1)  to  maintain  the 
required  alignment  of  the  seals  and  portholes, 
and  (2)  to  block  the  gas  flow  through  the  body’s 
most  downstream  portholes  when  the  ETVE  is  to 
be  in  the  close  positioa 

The  backward  configuration  model  runs 
show  that  the  gas  dynamic  forces  thrust  the  valve 


into  a  strongly  open  position  after  3  to  4  small 
bounces  on  the  Belleville  washers.  Beyond  this 
irtitial  opening  transient,  we  observe  that  the 
valve  maintains  the  open  position  with  a  slightly 
growing,  high-frequency  oscillation  imposed  on 
a  fairly  significant  steady  BelleviUe  washer  com¬ 
pression  (see  Figure  3-4). 

The  growing,  high-frequency  oscillation  is 
seen  to  correspond  with  the  gas  dynamic  force 
on  the  sleeve.  However,  the  net  average  gas 
dynamic  force  on  the  sleeve  is  shifting  toward 
the  open  position.  In  the  interest  of  time,  we  did 
not  run  the  model  any  more  time  steps  to  deter¬ 
mine  if  the  sleeve’s  apparent  growing  oscillation 
could  eventually  grow  to  the  point  of  blocking 
enough  of  the  flow  port  to  force  the  sleeve 
motion  into  another  transient  form.  We  recom¬ 
mend  performing  a  longer  run  (at  least  the  dura¬ 
tion  of  the  longest  expected  blast)  to  determine 
the  oscillations  trend  and  prove  or  disprove  the 
possibility  of  another  transient  in  the  open  cycle. 

Figure  3-4  plots  the  predicted  displacement 
of  the  sleeve  with  respect  to  time.  The  plot  was 
generated  by  rurming  the  model  with  the  sleeve 
in  the  backward  configuration  using  8% 
Belleville  washer  viscous  damping.  The  finer 
inlaid  plot  of  force  in  the  plot  represents  the  net 
force  on  the  moving  parts  (sleeve,  piston,  and 
piston  rod)  acting  in  the  axial  direction.  A  posi¬ 
tive  value  represents  a  force  in  the  upstream 
direction. 

3.2  Stress  Analysis 

Since  the  ETVE  was  not  designed  to  run  with 
the  sleeve  in  the  backward  configuration,  we 
performed  a  cursory  stress  analysis  on  the 
ETVE  stop  plate,  and  we  calculated  the  maxi¬ 
mum  allowable  width  for  the  body  portholes  for 
the  original  design  conditions  (700°F  and  1800 
psi).  Appendix  F  presents  this  portion  of  the 
work  in  greater  detail. 
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•  Sleeve  portholes  are  beveled  at  a  45-degree  angle  on  the  downstream  edge  of  each  porthole 

•  Sleeve  viscous  damping  set  at  8% 

•  Belleville  washers  Coulomb  damping  set  at  100  Ibf 

•  Portholes  are  fully  aligned  at  the  1.12  inch  marie  on  the  displacement  axis 

•  Belleville  washers  are  fully  compressed  at  the  1.2  inch  marie  on  the  displacement  axis:  1.12  inch¬ 
es  to  alignment  of  downstream  edges  of  portholes,  plus  0.05  inches  portholes  width  mismatch 
(washers  begin  to  be  compressed  at  1.17  inches),  plus  0.03  inches  maximum  washer  deflexion 
(form  Appendix  B  data  plots). 

•  Finer  inlaid  plot  of  force  represents  the  net  force  on  the  moving  parts  (deeve,  piston,  piston  rod, 
and  Belleville  washers)  acting  in  the  axial  direction.  A  positive  value  represents  force  in  the 
upstream  direction. 


Figure  3-2.  Predicted  displacement  history  for  the  sleeve  with  45-degree  beveled  portholes. 
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Figure  3-3.  Gas  dynamic  model  for  sleeve  in  the  backward  configuration. 


Displacement  (in) 


8%  Viscous  +  100  Ibf  Bellville  Damping 
Backwar(j  Opening  /  Lengthene(j  Spool 


Time  (nsec) 


•  Valve  opens  by  sliding  the  sleeve  in  the  upsteam  direction — backward  configuration 

•  Sleeve  viscous  damping  set  at  8% 

•  Belleville  washers  Coulomb  damping  set  at  100  Ibf 

•  Portholes  are  fully  aligned  at  the  1.12  inch  mark  on  the  displacement  axis 

•  Belleville  washers  are  fully  compressed  at  the  1 .2  inch  mark  on  the  displacement  axis:  1.12  inch¬ 
es  to  alignment  of  downstream  edges  of  portholes,  plus  0.05  inches  portholes  width  mismatch 
(washers  begin  to  be  compressed  at  1.17  inches),  plus  0.03  inches  maximum  washer  deflexion 
(form  Appendix  B  data  plots). 

•  Finer  inlaid  plot  of  force  represents  the  net  force  on  the  moving  parts  (sleeve,  piston,  piston  rod, 
and  Belleville  washers)  acting  in  the  axial  direction.  A  positive  value  represents  force  in  the 
upstream  direction. 


Figure  3-4.  Model  predictions  for  running  sleeve  in  the  backward  configuration. 
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Gas  Force  (Ibf) 


3.2.1  Wider  Portholes  In  Body.  As  the 

sleeve  portholes  oscillate  around  the  body  port¬ 
holes  during  the  opening  cycle,  the  effective  flow 
area  changes.  This  change  in  the  flow  area  gen¬ 
erates  a  segmented  shock  wave  which  is  unac¬ 
ceptable  for  blast  simulation. 

Since  the  real  requirement  is  to  maintain  a 
constant  effective  flow  area,  one  possible  config¬ 
uration  would  be  to  increase  the  width  of  the 
body  portholes  to  allow  the  sleeve  to  oscillate, 
yet  maintain  the  effective  flow  area  constant.  We 
analyzed  the  valve  to  determine  the  largest  width 
the  portholes  could  be  widened  to  and  still  main¬ 
tain  its  integrity  under  operating  conditions 
(700°Fand  1800  psi). 

The  stress  analysis  determined  that  the  port¬ 
holes  in  the  ETVE  body  could  be  widened  to 
0.76  in.  centered  on  the  existing  portholes  cen¬ 
terlines,  as  shown  in  Figure  3-5.  This  means  that 
the  sleeve’s  centerline  can  oscUlate  at  an  ampli¬ 
tude  of  0.155  in.  to  either  side  of  the  body  port¬ 
hole  centerline  and  stiU  maintain  a  constant 
effective  flow  area. 

3.2.2  Stress  analysis  on  stop.  When 
the  ETVE  is  actuated  open  in  the  current  config¬ 
uration,  the  actuating  piston  slides  to  the  back 
(upstream)  end  to  impact  on  the  end  plate.  How¬ 
ever,  in  the  proposed  backward  configuration, 
the  actuating  piston  will  move  to  the  fi’ont 
(downstream)  to  impact  on  the  stop,  exposing  it 
to  much  larger  forces  than  it  was  originally 
designed  for. 


The  analysis  revealed  that  the  stop  would  expe¬ 
rience  stresses  over  400  ksi;  the  six  bolts  fasten¬ 
ing  the  stop  to  the  body  would  experience  224.9 
ksi;  and  the  ETVE  body  around  the  portholes 
would  experience  stresses  in  excess  of  58  ksi. 
All  these  stresses  exceed  the  maximum  allow¬ 
able  stress  of  32  ksi  for  these  components 
(defined  in  the  ASME  Boiler  and  Pressure  Vessel 
Code,  Section  VIII,  Rules  for  Construction  of 
Pressure  Vessels). 

The  large  dynamic  (impulse)  forces  created 
by  the  deceleration  of  the  piston  and  sleeve  are 
aggravated  by  the  stif&iess  of  the  Belleville 
washer  stack;  the  stack  is  clearly  very  stiff  and 
has  a  short  distance  to  travel  before  bottoming 
out.  However,  the  piston  impulse  forces  can  be 
greatly  reduced  if  either  of  the  following  modifi¬ 
cations  are  made; 

•  Arrange  (or  add  to)  the  stack  of  Belleville 
washers  so  that  the  effective  travel  length 
(before  bottoming  out)  is  much  longer.  Note 
that  this  recommendation  is  based  on  the 
observations  from  the  washer  characteriza¬ 
tion  tests;  when  the  washers  bottom-out, 
their  stiffness  increases  tremendously,  which 
in  turn  increases  the  piston/sleeve  decelera¬ 
tions. 

•  Employ  Belleville  washers  (or  other  springs) 
that  have  a  lower  stiffiiess  value  and  a  longer 
travel  length.  This  will  reduce  the  magni¬ 
tude  of  the  piston/sleeve  deceleration. 
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4.  SEALS  STUDY 


The  1991  ETVE  tests  at  the  INH^  revealed 
another  area  of  the  ETVE’s  perfonnance  that  needs 
further  development:  the  main  valve  seals.  Sever¬ 
al  kinds  of  seals  were  used  during  testing,  tnit  none 
were  successful  under  the  operating  conditions 
(700°F  and  2000  psi).  We  performed  a  cursory 
seal  manufacturer  and  material  study  to  determine 
the  available  seal  materials  and  configurations 
capable  of  withstanding  the  operating  conditions. 
Appendix  G  presents  this  work  in  detail. 

The  study  concludes  that  no  conunercial 


seals  are  readily  available  that  will  withstand  the 
operating  conditions.  We  recommend  testing 
different  types  of  seals,  including  split-metal  pis¬ 
ton  ring  type  seals  and  VESPEL  seals.  Appen¬ 
dix  G  presents  a  Draft  Proposed  Seal  Test  Pro¬ 
gram  to  evaluate  existing  seal  technologies. 

All  vendor  information  gathered  will  be 
maintained  in  EG&G  Idaho,  Inc.,  Project  File 
93-713,  and  a  copy  will  be  transmitted  to  the  U. 
S.  Army  Research  Laboratory  technical  monitor, 
Mr.  Richard  J.  Pearson. 
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5.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  ETVE  simulation  model  has  proven  to 
be  an  extremely  valuable  tool  in  assessing  the 
qualitative  nature  of  the  current  valve’s  operation 
and  is  indispensable  in  assessing  the  elfects  of 
valve  modifications  or  redesigns.  The  ETVE’s 
complex  nonlinear  behavior  is  impossible  to  pre¬ 
dict  in  any  way  other  than  numerical  simulation. 
Even  if  the  numerical  simulations  are  accepted 
on  a  qualitative  basis  only,  they  increase  under¬ 
standing  of  behavior,  show  trends  and  sensitivi¬ 
ties,  and  suggest  alternative  designs.  On  the 
other  hand,  the  alternative,  iterative  “cut  and  try” 
approach  to  design  would  be  a  long,  tedious,  and 
expensive  route  to  valve  design/redesign.  In 
addition,  the  model  can  be  used  to  input  new 
BelleviUe  washer  characterization  data  for  differ¬ 
ent  washer  configurations  and  then  assess  the 
qualitative  effects  of  each  configuration  on  the 
overall  valve  behavior. 

For  future  design/analysis  work  and  to  pre¬ 
dict  accurate  quantitative  performance,  we  rec¬ 
ommend  that  the  numerical  technique  developed 
to  simulate  the  gas  dynamics  be  modified  in  four 
main  areas: 

1.  Modify  the  model  to  give  a  partially  implic¬ 
it  time  integration  instead  of  the  explicit  time 
integration  (this  would  reduce  the  CPU  time 
required  to  make  runs  and  should  increase 
the  accuracy  by  reducing  the  laige  number  of 
time  steps  required  for  the  model  runs) 

2.  Examine  the  importance  of  viscous  effects, 
especially  in  the  valve-open  position  where 
the  high-speed  flow  through  the  ports  could 
possibly  place  significant  shear  forces  on 
the  sleeve,  and,  therefore,  influence  its 
motion 

3.  Modify  the  model  to  study  the  ETVE  close 
cycle  under  operating  conditions 

4.  Design  the  valve’s  effective  damping  with 
either  BeUeviUe  washers  or  another  damping 
system  for  both  the  open  and  close  cycles. 


5.1  PROPOSED  MODIFICATION 

We  propose  the  ETVE  configuration  change 
(Figure  5-1)  as  a  preliminary  concept  only,  and 
strongly  caution  that  further  design  and  analysis 
work  must  be  performed  prior  to  executing  any 
actual  modification  to  the  existing  ETVE.  Fur¬ 
thermore,  we  recommend  taking  advantage  of 
the  already  developed  valve  sleeve  dynamic  sim¬ 
ulation  model  to  characterize  the  performance  of 
the  ETVE  under  its  proposed  modified  configu¬ 
ration  or  any  other  configuration  that  may  be 
developed  by  further  analysis  and  research. 

The  study  described  in  Section  3.2.1  analyz¬ 
ing  widened  portholes  in  the  body  is  left  as  a 
matter  of  consideration  for  the  next  redesign 
phase.  Time  precluded  our  running  the  model 
with  widened  portholes.  Since  we  do  not  have 
model  predictions  for  this  configuration,  we  do 
not  include  widened  portholes  as  part  of  the  pro¬ 
posed  modification  concept. 

The  preliminary  concept  makes  reference 
only  to  correcting  the  chatter  in  the  open  cycle; 
however,  no  analysis  was  performed  on  the  close 
cycle.  The  proposed  concept  evolved  from  the 
results  of  the  dynamic  simulation  model  of  the 
open  cycle,  the  stress  analysis  performed  on  the 
stop  plate  and  the  valve  body,  and  from  the 
Belleville  washers  characterization  studies. 

The  modification  concept  is  best  illustrated 
in  Figure  5-1.  Its  main  areas  of  redesign  follow. 

1.  Increase  the  length  of  the  sleeve  by  approxi¬ 
mately  1.6  in.  and  add  two  sets  of  seals  in  the 
added  length  to  seal  the  last  row  of  body 
portholes.  Reduce  the  length  of  the  liner  by 
0.7  in. 

2.  Place  the  sleeve  upstream  (from  its  current 
as-built  location)  as  indicated  in  Rgure  5-1 
so  the  pneumatics  can  slide  the  sleeve  for¬ 
ward  to  its  open  position. 
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Figure  5-1 .  Proposed  preliminary  concept  for  ETVE  configuration. 


3.  Increase  the  length  of  the  ETVE  body  to 
house  approximately  32  Belleville  washers 
on  the  open  chamber  and  8  washers  in  the 
close  chamber.  This  modification  requites 
letapping  the  pneumatic  actuation  lines. 

4.  Lengthen  the  piston  rod  to  accommodate  the 
added  number  of  washers. 

5.  Increase  the  stop  plate  size  and  increase  the 
number  and  size  of  bolts  used  to  fasten  the 
stop  plate  to  the  valve  body. 

A  fiill  set  of  nonreleased  drawings  for  this 
proposed  configuration  is  included  in  the  engi¬ 
neering  design  file  at  EG&G  Idaho,  and  a  copy  is 
being  sent  to  U.  S.  Army  research  Laboratory 
technical  monitor,  Mr.  Richard  J.  Pearson. 

5.2  CONCLUSIONS 

1 .  The  ETVE  sleeve  dynamic  simulation  model 
confirms  that  the  ETVE  does  imdergo  chat¬ 
tering  during  the  opening  cycle.  In  fact,  we 
observe,  from  Figure  2  no  indication  that  the 
total  amplitude  of  the  sleeve  oscillations 
decreases  within  the  first  40  msec.  We  also 

observed,  during  the  1991^  testing  at  the 
INEL,  that  the  duration  of  the  deeve  oscilla¬ 


tion  increased  as  the  driver  gas  pressure 
increased.  From  these  two  observations,  we 
deduce  that  the  sleeve  will  continue  to  oscil¬ 
late  for  as  long  as  there  is  a  relatively  high 
pressure  driver  gas  flow  through  its  port¬ 
holes.  The  model  reveals  that  the  main  dri¬ 
ving  force  causing  the  Chatteris  generated  by 
the  driver  gas  as  it  flows  throu^  the  sleeve 
portholes.  The  model  also  reveals  that  the 
ETVE  does  not  have  sufficient  damping  to 
absorb  energy  from  the  moving  sleeve. 

2.  The  ETVE  should  not  be  used  imder  its  cur¬ 
rent  configuration  in  high-pressure  or  high- 
temperature  gas  conditions.  Furthermore, 
the  current  ETVE  configuration  requires  sig¬ 
nificant  redesign,  analysis,  and  refabrication 
to  meet  its  original  design  requirements . 

3.  The  close-cycle  actuation  needs  to  be  evalu¬ 
ated  and  redesigned  to  overcome  the  driver 
gas  forces  that  will  make  the  valve  diflBcult 
to  close. 

4.  Any  general  design  concepts  should  be  eval¬ 
uated  by  using  the  already  developed  valve 
sleeve  dynamic  simulation  model  to  qualita¬ 
tively  predict  and  evaluate  the  design  con¬ 
cept’s  effects  on  the  ETVE  performance. 


Intentionally  left  blank. 
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Appendix  A 

Eaton  Throat-Valve  Element  Prototype  1993  Testing 


I  perfonned  two  sets  of  tests  in  November  1993 
at  the  Idaho  National  Engineering  Laboratory 
(INEL)  to  characterize  the  behavier  of  the  Eaton 
Throat- Valve  Element  (ETVE)  under  existing  con¬ 
figuration.  The  first  test  characterized  the  ETVE 
sleeve’s  displacement  behavior  with  a  linear  vari¬ 
able  displacement  transducer  (LVDT).  The  LVDT 
data  then  contributed  to  J.  G.  Arendts’  analytical 
development  of  a  computer  model  of  the  mechani¬ 
cal  dynamics  of  the  ETVE  prototype  moving  parts 
(see  Apendk  D).  The  second  test  determined 
whether  different  sequencing  of  the  open  and  close 
actuating  sloenoid  valves  could  reduce  the  valve 
chatter  observed  in  the  initial  testing  of  the  ETVE 
prototype  at  the  INEL. 

This  report  describes  and  documents  the  two 
sets  of  tests.  The  system  setup  was  slightly  differ¬ 
ent  for  each  test,  so  a  block  diagram  of  each  con¬ 
figuration  is  presented  in  the  following  discussiond 
for  each. 

C2.  CHARACTERIZATION  OF  THE  ETVE 

The  initial  ETVE  tests  performed  at  the  INEL 
indicate  that  the  ETVE  was  possibly  chattering  dur¬ 
ing  the  opening  cycle  berfore  achieving  steady-state 
open  position.  However,  the  instrumentation  avail¬ 
able  at  the  time  could  not  give  a  detailed  dis- 
placemnt  history  of  the  ETVE’s  sliding  sleeve. 
Therefore,  in  order  to  subsequently  characterize  the 
suspected  chatter,  a  linear  variable  displacement 
transducer  (LVOT)  was  purchased  and  installed 
(STARRET  Model  254,  Serial  Number  702168), 
and  the  ETVE  characterization  test  was  performed 
to  obtain  a  high  sampling  of  the  ETVE’s  sliding 
sleeve  during  the  open  cycle.  The  testing  was  per¬ 
formed  at  the  INEL’s  North  Holmes  Laboratory  in 
Idaho  Falls,  Idaho.  Since  the  pressure  tube  that  nor¬ 
mally  is  pressurized  during  the  Large 
*  Blast/Thermal  Simulation  tests  was  not  available 
for  use  at  the  INEL  facility,  the  tests  were  per¬ 
formed  only  at  ambient  temperature  and  pressure. 


Figure  A-1  shows  the  hardware  configuratioa 
The  data  system  consisted  of  a  Kaye  MDAS  1000 
data  system  with  a  20-Mbyte  hard  drive  and  a  1.2- 
Mbyte  floppy  disk.  Data  were  collected  to  memo¬ 
ry  then  stored  in  ascii  format  to  hard  disk  and  flop¬ 
py.  The  LVDT  and  pressure  transducers  where  both 
calibrated  (see  Attachment  1),  and  the  data  were 
converted  to  the  corresponding  values  in  engineer¬ 
ing  units  before  being  saved  to  disk.  Two  separate 
tests  were  performed  under  this  configuration. 

2.1.2  CHARACTERIZATION  TEST  1. 

Under  this  test  condition,  1200  psi  were  used  to 
open  the  ETVE.  The  data  acquisition  system  was 
then  used  to  measure  the  sleeve  displacement, 
open  pressure  used  and  fire  signal  from  the  Eaton 
control  box.  When  the  fire  button  was  pressed  on 
the  Eaton  control  box,  a  signal  was  sent  to  the 
data  acquisition  system  to  start  collecting  data  at 
10,000  samples  per  second  for  1.5  seconds. 
Attachment  2  contains  the  results  of  this  test 

2.1.2  CHARACTERIZATION  TEST  2. 

Under  this  test  condition  2000  PSI  was  used  to 
open  the  ETVE.  The  data  acquisition  system  was 
then  used  to  measure  the  sleeve  displacement, 
open  pressure  and  fire  signal  from  the  Eaton  con¬ 
trol  box.  When  the  fire  button  was  pressed  on  the 
Eaton  control  box,  a  signal  was  sent  to  the  data 
acquisition  to  start  collecting  data  at  10,000  sam¬ 
ples  per  second  for  1.5  seconds.  Attachment  3 
contains  the  results  of  this  test. 

2.2  Damping  Testing 

The  characterization  test  confirmed  that  the 
ETVE  chattered  during  the  open  cycle.  This  test  set 
was  performed  to  determine  whether  pressuring  the 
close  side  of  the  valve  would  dampen  the  chatter. 
During  the  testing,  both  the  close  and  openactuating 
pressures  and  the  sleeve  displacemtn  were  mea¬ 
sured.  The  data  system  was  triggered  from  the 
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LVDT  movement  (see  Figure  A-2).  The  system 
was  modified  so  that  the  Eaton  control  box  con¬ 
trolled  only  the  ETVE  opening.  The  close  solenoid 
was  controlled  by  the  Kaye  data  system  through  a 
solenoid  state  relay.  A  28  volt  dc  power  supply 
was  connectted  across  this  relay  to  actuate  the 
close  solenoid.  Wth  this  control  system,  it  was 
possible  to  apply  a  pressure  to  the  close  side  of  the 
valve  and  then  release  that  pressure  an  arbitrary 
time  from  when  the  valve  moved.  Two  tests  were 
run  with  this  configuration. 

2.2.1  Damping  Test  1 .  Test  1  was  set  so 

that  the  close  pressure  would  be  relieved  in  the 


middle  of  the  valve  opening  stroke.  See  appendix 
D  for  these  test  results. 

2.2.2.  Damping  Test  2.  Test  2  was  set  so 

that  the  close  pressure  was  left  on  for  0.5  seconds 
after  the  valve  was  actuated  open.  This  configura¬ 
tion  allowed  the  valve  fuUy  open  and  stabilize 
before  the  close  pressure  was  released.  See 
Attachment  5  for  the  plotted  test  results.  Both  test 
used  a  open  pressure  of  2000  psi  and  a  close  pres¬ 
sure  of  1000  psi.  The  data  was  collected  at  10,000 
samples  per  second  for  1  sec  for  both  of  these  tests. 
Figure  A-3  shows  the  ETVE’s  sleeve  displacement 
overlayed  for  both  damping  tests. 
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re  A-1 .  Characterization  test  setup. 
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Figure  A-2.  Damping  test  setup. 
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Figure  A-3.  ETVE’s  sleeve  displacement  overlayed  for  both  damping  tests. 
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PRESSURE  TRANSDUCERS 
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Attachment  1 


EG&S  Idaho,  Inc. 

STANDARDS  &  CALIBRATION  LABORATORY 
CALIBRATION  DATA  SHEET 

DATE:  08/24/1993 


MANUFACTURER:  NOMENCLATURE:  MODEL  #:  PROPERTY  #: 

UIATRAN  PRESSURE  218  703192 

TRANSDUCER 

SERIAL  #:  PROCEDURE  #:  NOMINAL  INTERVAL:  NEXT  DUE  DATE: 


295S78 

3i 

043 

15 

MONTHS 

1  1/2 

4/1994 

PRESSURE 
HIGH  SIDE 

NOMINAL 

mV 

ACTUAL 

mV 

INCREASING 

ACTUAL 

mV 

DECREASING 

mV 

HYSTERISIS 

LINEARITY 

LOWER  UPPER 

LIMIT  LIMIT 

0 . 00000 

-3.9000 

-3.9000 

“3.2000 

-0.7000 

-15.9064 

8.1064 

300 . 00000 

29B.2B00 

295.8000 

29B.8000 

-1 .0000 

284.2536 

308.2664 

B00. 00000 

59B.4200 

595.8000 

596.9000 

-1.1000 

584.4136 

B08.42B4 

900,00000 

89G.5800 

896.0000 

896.9000 

“0.9000 

884.5736 

908.5864 

1 200 . 00000 

1 19B.7400 

1196.2000 

1197.5000 

“1 .3000 

1184.7336 

1208.7464 

1 E00 . 00000 

1498.9000 

149G.3000 

1497.1000 

“0.8000 

1484.8936 

1508.9064 

1 800.00000 

1797.0600 

179B.6000 

1797.3000 

-0.7000 

1785.0536 

1809.0664 

2100-00000 

2097.2200 

2096.8000 

2097.4000 

“0 . B000 

2085.2136 

2109.2264 

2400 . 00000 

2397.3800 

2397.2000 

2397. B000 

“0 . 4000 

2385.3736 

2409.3864 

2700.00000 

2G97.5400 

2697.5000 

2697.7000 

“0 , 2000 

2685.5336 

2709.5464 

3000.00000 

2997.7000 

2997.7000 

2997.7000 

0.0000 

2985.6936 

3009.7064 

THE  HYSTERESIS,  IN  mV,  MUST  BE  LESS  THAN  12.00B4 


PRESSURE  IN  UNITS  OF  PSI 
FULL  SCALE  PRESSURE:  3000.0000  PSI 

LINEARITY  ACCURACY:  0.4000  % 

HYSTERESIS  ACCURACY:  0.4000  % 

EXCITATION  VOLTS:  24.0000  VOL^ 

TECHNICIAN:  K .  ANDER50N-.-/0-T 

CALIBRATED  FOR:  8. COLSON 

REMARKS:  A-9 

INPUT  1+  2-/OUTPUT  3+ 
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EG&S  Idaho,  Inc. 

STANDARDS  &  CALIBRATION  LABORATORY 
CALIBRATION  DATA  SHEET 

DATE:  08/24/1993 


MANUFACTURER:  NOMENCLATURE:  MODEL  #:  PROPERTY  S: 

UIATRAN  PRESSURE  218  705191 

TRANSDUCER 

SERIAL  S:  PROCEDURE  #:  NOMINAL  INTERUAL:  NEXT  DUE  DATE: 


295478 

3045 

15 

MONTHS 

11/24/1934 

PRESSURE 
HI6H  SIDE 

NOMINAL 

nU 

ACTUAL 

nU 

INCREASING 

ACTUAL 

nU 

DECREASING 

nU 

HYSTERISIS 

LINEARITY 

LOWER  UPPER 

LIMIT  LIMIT 

0 . 00000 

2.3000 

2 .5000 

3.4000 

-1 .1000 

-9. 8880 

14.2880 

500 • 00000 

302 . 0000 

502 . 1 000 

505.4000 

-1 .3000 

290.0120 

313.9880 

600,00000 

601.7000 

602.0000 

B03.5000 

-1 .3000 

589.7120 

613.6880 

900.00000 

901 .4000 

901.9000 

903.2000 

-1 .3000 

889.4120 

913.3880 

1 200.00000 

1201.1000 

1201.8000 

1203.0000 

-1 .2000 

1 189.1120 

1213.0880 

1500.00000 

1500.8000 

1501.5000 

1502.8000 

-1 .3000 

1488.8120 

1512.7880 

1 800.00000 

1800.5000 

1801 . 1000 

1802.3000 

-1 .2000 

1788.5120 

1812.4880 

2100.00000 

2100.2000 

2100.6000 

2101 .6000 

—1 . 0000 

2088.2120 

2112.1880 

2400.00000 

2399.9000 

2400.2000 

2401 .0000 

-0.8000 

2387.9120 

241 1 .8880 

2700.00000 

2699.6000 

2699.7000 

2700.0000 

-0.3000 

2687.6120 

271 1 .5880 

3000 • 00000 

2999.3000 

2999.3000 

2999.3000 

0 . 0000 

2987.3120 

301 1 .2880 

THE  HYSTERESIS,  IN 

nU,  MUST  BE 

LESS  THAN 

1  1.9880 

PRESSURE 

IN  UNITS 

OF  PSI 

FULL  SCALE  PRESSURE:  3000.0000  PSI 

LINEARITY  ACCURACY:  0.4000  % 

HYSTERESIS  ACCURACY:  0.4000  % 

EXCITATION  UOLTS:  24.0000  UOLTS 

TECHNICIAN:  K .  ANDERSON^^^ 

CALIBRATED  FOR:  B. COLSON 

REMARKS: 

INPUT  1+  2-/OUTPUT  5+  4- 

A;  .7/. 
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vitran  703191 


0 

0.0023 

300 

0.3021 

600 

0.602 

900 

0.9019 

1200 

1.2018 

Regression  Output: 

1500 

1.5015 

Constant  -3.40787 

1800 

1.8011 

Std  Err  of  Y  Est  0.685176 

2100 

2.1006 

R  Squared  1 

2400 

2.4002 

No.  of  Observations  22 

2700 

2.6997 

Degrees  of  Freedom  20 

3000 

2.9993 

3000 

2.9993 

X  Coef f icient (s)  1001.186 

2700 

2.7 

Std  Err  of  Coef.  0.154164 

2400 

2.401 

2100 

2.1016 

1800 

1.8023 

1500 

1.5028 

1200 

1.203 

900 

0.9032 

600 

0.6033 

300 

0.3034 

0 

0.0034 

vitran  703192 


0 

-0.0039 

300 

0.2958 

600 

0.5958 

900 

0.896 

1200 

1.1962 

Regression  Output: 

1500 

1.4963 

Constant 

7.544301 

1800 

1.7966 

Std  Err  of  Y  Est 

5.557849 

2100 

2.0968 

R  Squared 

0.999969 

2400 

2.3972 

No.  of  Observations 

22 

2700 

2.6975 

Degrees  of  Freedom 

20 

3000 

2.9977 

3000 

2.9977 

X  Coefficient (s)  997.8521 

2700 

2.6977 

Std  Err  of  Coef.  1.24637 

2400 

2.3976 

2100 

2.0974 

1800 

1.7973 

1500 

1.4971 

1200 

1.1975 

900 

0.8969 

600 

0.5969 

300 

0.2698 

0 

-0.0032 
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Idaho 


STANDARDS  &  CALIBRATION  LABORATORY 

CALIBRATION  DATA  SHEET 


FORM  £G&G-t510B 
{Rev, 


Pa^Qe  j  of  ^ 

0*Tf _ 


manufacturer  L^iAQPiS/ 

nomenclature 

model  no. 

serial  no. 

SCHAEVITZ 

LVDT 

/Do/:>  ncj\~’xi^i 

yj.c>7ci 

PROPERTT  NO. 

PROCEDURE  NO. 

INTERVA  L 

CALIBRATED  FOR: 

7 1  lyciH- 

3SSD 

/  7.  MP. 

CALIBRATION  EOUIPMENT  USED 

STANDARDS  CERTIFIED  TRACEABLE  TO  THE  NATIONAL  BUREAU  OF  STANDARDS 
(Records  on  file  m  the  EG&G  Standards  &  Calibration  Laboratory.  CF-698) 


manufacture  R 

MODEL  HO. 

serial  no. 

recall  oau 

manufacturer 

model  no. 

serial  ho . 

recall  date 

starrett 

7^'XIA^S 

9'  O.S-f3 

FLUKE 

SOSOA 

7c>/c^5a 

IOHJ-J3 

lljjjllllfllllllll 

HI^HH 

inifmm 

i|H|||||^ 

CALIBRATION  DATA 


PROCEDURE 

SUP 

NO. 


rUNCTION  TESTES 


jtM 


t: 


■Support  Arm  ' 
Holding  Rod 

— Core 
-LVDT 


IN 

TOLERANCE 
.  CA  ) 


L  VDT  UJ.±c: 


1£ 


Mil 


He. 


Angle 


-Vee  Block 


JEk 


Angle 


yX2. 


/Ajas  yrk 


fh .  'the. 


our  or 
tolerance 

(READING; 


AO  JUS  TKENT 
tCA) 


C.]jU21 


le. 


tolerances 


LOVER  LIMIT*  UPPER  L IMIT 
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Ll)At  Y1  a cTeH^. 


ver  J±.± 


t  to  AS  -Ec. 


TYPE  Of  DATA; 

□  TEST  REfERRAL 
P.fl .  #  _ 


O  calibration 
Q  certification 


Q  Evaluation 

O  cross  ckecx 
D  SPECIAL  measurement 


teckkic IAn 


MRS  . 

^rC> 


DATE  next  cal  IBRA  I  ION  OUE 


yxs  Loer's 


IC)'^  riF  C. 


OV\ 


P\Vyiy-iriflpy'  jS/AJ- 
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Ivdt  711724 


Attachment  1 


0  10 


0.1 

8.984 

0.2 

7.97 

0.3 

6.955 

0.4 

5.947 

Regression  Output: 

0.5 

■  ■  4.943 

Constant 

1.000416 

0.6 

3.944 

Std  Err  of  Y  Est 

0.002634 

0.7 

2.953 

R  Squared 

0.999983 

0.8 

1.973 

No.  of  Observations 

21 

0.9 

0.993 

Degrees  of  Freedom 

19 

1 

0 

1.1 

-0.986 

X  Coef f icient (s)  -0.10054 

1.2 

-1.972 

Std  Err  of  Coef.  0.000095 

1.3 

-2.952 

1.4 

-3.945 

1.5 

-4.94 

1.6 

-5.941 

1.7 

-6.947 

1.8 

-7.956 

1.9 

-8.961 

2 

-9.9752 
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ETVE  CHARACTERIZATION  TEST  DATA  PLOTS 


TEST1 


Open  actuating  pressure:  1200  psi 
Data  sampling  rate:  10,000  samples  per  second 
Data  sampling  time:  1.5  seconds 


;  Name:  _ FIRE1200  BH  Sample  Rate:  le+004 

I  Number:  1  |H  Maximum:  9.99512 

Name:  12ii^B|  Minimum:  -0.551758 
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A-16 


Dataset  Name:  Unspecified  H|  Sample  Rate:  l&f004 
Uersion  Number:  0  B  M^imum:  1.72108 
Series  Name:  Unspecified  Hj  Minimum:  0.576754 


A-19 


Attachment  3 


ETVE  CHARACTERIZATION  TEST  DATA  PLOTS 

TEST  2 

Open  actuating  pressure:  2000  psi 
Data  sampling  rate:  10,000  samples  per  second 
Data  sampling  time:  1.5  seconds 
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A-22 


Dataset  Hawe^  OPEH2000  ||||  Sample  Rate:  lef004 
Uersion  Number:  1  W  Maximum:  2499.48 
Series  Name:  2000PSI  H|  Minimum:  -1333  18 
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Dataset  Hawe:  Unspecified  lBt004 

Ueiraion  Humber:  0  W  Haxiimm:  1,71961 


■'f'-r  :  ;V'.r, 


nr- 


li- i 

I-:-  ••  -5. 

,  .  ■■■,■■  . 

;¥;-::¥v^¥Sl;'r¥| 
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05  O  5  j 


f¥;¥v|g#; 


Dataset  Natne:  Unspecified  M  Sample  Rate:  lBt004 

Uersion  Humber:  0  ||P||  tlaximum:  1.6931 

Series  Name:  Unspecified  H  Minimum:  ~~  0.567917 
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ETVE  DAMPING  TEST  DATA  PLOTS 
TEST1 


Open  actuating  pressure:  2000  psi 
Close  actuating  pressure:  1000  psi 
Close  actuating  pressure — release  time:  70  msec 
from  start  of  sleeve  displacement 
Data  sampling  rate:  10,000  samples  per  second 
Data  sampling  time:  1  second 
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Dataset  Name:  OPEN  ||||  Sample  Rate:  l&f004 
UersioTfi  Number:  1  |l||  Maximum:  1998.7 
Series  Name:  CIS  PRES  ON  Hi  Minimum:  1201.33 
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Attachment  4 


Dataset  Nawe:  Unspecified  WM  Sample  Rate:  le+004 
Uersioifi  Number!  0  ||||  Maximum!  1.69507 
Series  Name!  Unspecified  Hj  Minimum!  1.6715 
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ETVE  DAMPING  TEST  DATA  PLOTS 
TEST  2 

Open  actuating  pressure:  2000  psi 
Close  actuating  pressure:  1000  psi 
Close  actuating  pressure — release  time:  0.5  sec 
Data  sampling  rate:  10,000  samples  per  second 
Data  sampling  time:  1  second 
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Dataset  Hawe:  UIBT  Wm  Sample  Rate:  _  ^  le+004 

Uersion  Number:  1  mm  Maximum:  1.70341 
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set  Name: _ Unspecified 

ion  Number:  R 
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APPENDIX  B 

BELLEVILLE  WASHER  CHARACTERIZATION  TESTING 

G.  E.  Korth 


B-l 


Intentionally  left  blank. 


B-2 


Objective  of  Tests:  These  tests  were  conducted  to  determine  the  quasi -static 
load-deflection  behavior  of  different  configurations  of  Belleville  washers  used 
or  proposed  to  be  used  in  the  Eaton  Throat  Valve  prototype.  This  loading 
behavior  provides  input  for  the  numerical  modeling  task  of  the  dynamic  valve 
evaluation. 

Description  of  Material:  The  test  specimens  were  made  from  17-7PH  CH900  CRES 
(AMS  5529).  They  were  nominally  1.00  inch  in  diameter  with  a  0.50  inch  hole  and 
0.073  inches  thick.  The  part  number  of  these  washers  is  B1000-073-S.  A  free 
standing  height  of  one  washer  without  load  is  0.088  to  0.090  inches. 

Description  of  Test  Eouipnient:  The  washers  were  tested  in  a  compression  mode 
using  an  Instron  Model  1128  Universal  load  frame.  This  load  frame  is  a  four  post 
screw-driven  machine  with  maximum  capacity  of  100,000  pounds.  The  washers  were 
loaded  between  two  compression  platens.  The  accuracy  of  the  load  cell  is  within 
±  0.5%  of  the  indicated  load  or  ±  0.25%  of  the  load  range  in  use.  Deflection  was 
monitored  between  the  compression  platens  with  a  linear  variable  capacitance 
transducer  with  an  accuracy  of  within  0.02%  of  full  range  (0.50  inches). 

Test  Procedure:  The  washers  were  tested  either  dry  (as  received)  or  lubed  with 
DN-600  Polar  Start  grease  manufactured  by  Conoco.  They  were  placed  directly  on 
the  compression  platen  in  various  configurations  as  shown  in  Figure  1.  The  load 
frame  cross-head  was  manually  moved  down  until  the  first  indication  of  load  was 
noted  and  then  the  electronic  cross-head  positioner  was  zeroed  so  upon  unloading 
the  cross-head  would  come  back  to  that  same  position.  Using  this  technique, 
permanent  deformation  would  be  evident  on  unloading  or  in  subsequent  cycles  of 
the  same  stack.  The  cross -head  was  programmed  to  move  with  a  speed  of  0.05  to 
0.50  inches/minute.  After  the  desired  load  was  reached  the  reverse  was  activated 
and  the  stack  was  unloaded  with  the  same  cross-head  velocity  (0.05  to  0.50 
in./min.).  Compliance  was  measured  on  the  system  by  loading  the  compression 

platens  without  any  washers  present  (See  Chart  16).  The  compliance  of  the 
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machine  can  probably  be  neglected  since  it  was  measured  to  be  1.25  x  10' 
inches/pound. 


B-3 


The  same  set  of  washers  were  used  for  the  dry  and  lubed  tests  of  configuration 
(a)  and  (b),  but  for  all  other  configurations,  a  new  set  of  as-received  washers 
were  used  with  one  exception:  configuration  (1)  used  four  washers  which  had  also 
been  tested  in  the  configuration  (k)  stack.  The  load  and  deflection  were 
recorded  on  an  X-Y  plotter  with  the  load  range,  load  calibration,  and  deflection 
calibration  noted  on  each  chart.  The  test  temperature  was  ambient  laboratory 
temperature  (approximately  22-24°C) . 

Results :  •  The  loading  and  unloading  behavior  of  the  various  stack  configurations 
shown  in  Fig.l  is  described  by  the  X-Y  plots.  The  matrix  of  tests  conducted  is 
listed  in  Table  1  which  also  lists  the  maximum  load  and  deflection  for  each 
configuration  stack. 
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Configuration  (a) 
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Configuration  (b) 
Stack  of  4 
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k  of  12 


Configuration  (d) 
Stack  of  16 


Figure  1.  Stack  configurations  used  in  testing  Belleville  washers 


Configuration  (e),  Stack  of  20 


Configuration  (f),  Stack  of  12 
Configuration  (g),  Stack  of  16 
Configuration  (h),  Stack  of  20 


Configuration  (i),  Stack  of  20  Configuration  (j),  Stack  of  12 

Configuration  (k),  Stack  of  16 
Configuration  (1),  Stack  of  20 


Figure  1.  (cont'd)  Stack  configurations  used  in  testing  Belleville  washers 
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Table  1.  Belleville  washer  load-deflection  stack  tests. 


Stack 
Configura¬ 
tion 
(  Fig.l) 

Configuration 

Description 

Lube 

Max. 

Load 

(lbs) 

Max. 

Defl. 

10’^  in. 

No. 

Cycles 

Chart 

No. 

a 

Single 

N 

3,500 

14.5 

3 

1 

a 

Single 

Y 

3,500 

14.5 

3 

1 

b 

4  Parallel 

N 

14,000 

16.1 

3 

2 

b 

4  Parallel 

Y 

14,000 

17.3 

10 

2,3 

c 

12(Two  6's) 

N 

21,000 

34.2 

3 

4 

c 

12(Two  6's) 

Y 

21,000 

38.8 

10 

5,6 

d 

16(Two  8's) 

Y 

45,000 

39.4 

6 

7,8 

e 

20 (Two  10' s) 

Y 

70,000 

65.0 

5 

9 

f 

12{Alternating) 

Y 

4,500 

176.5 

5 

10 

9 

16(Alternating) 

Y 

4,500 

228.5 

5 

11 

h 

20(Alternating) 

Y 

4,500 

264.0 

5 

12 

i 

20(Ten  sets  of 
alternat.  2's) 

Y 

14,000 

146.5 

5 

13 

j 

12  Parallel 

Y 

70,000 

29.3 

5 

14 

k 

16  Parallel 

Y 

80,000 

42.2 

5 

15 

1 

20  Parallel 

Y 

100,000 

74.0 

5 

16 
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APPENDIX  C 


BELLEVILLE  WASHER  DAMPING  CHARACTERIZATION 


Unloaded  initial  displacement  ringdown  of  the 
BelleviUe  washers  is  examined  herein  to  ascertain 
the  damping  behavior  of  the  Eaton  Throat- Valve 
Element  Prototype  (ETVE).  The  ideal  spring 
behavior  assumed  is  that  described  in  Appendix  A 
of  this  report.  I  wrote  a  simple  FORTRAN  pro¬ 
gram  to  use  the  ETVE  mechanical  motion  subrou¬ 
tines  described  in  Appendix  D.  Ringdown  behav¬ 
ior  of  both  ETVE  opening  and  closing  BelleviUe 
washer  stacks  is  observed;  we  see  that  overload  of 
either  washer  significantly  reduces  its  effective 
damping. 


where 

y(0)  =  The  initial  displacement 
Fc  =  The  Coulomb  damping  force 
M  =  The  mass  of  the  system 
T  =  The  natural  vibration  period. 

Thus,  if  the  mass  is  known  and  the  period  is 
measured  from  the  actual  decay  curve,  the  equiva¬ 
lent  Coulomb  damping  force  may  be  found 
through  substituting  y(t)  and  a  relative  maxima  and 
corresponding  time  into  the  above  relationship. 


Ringdown  Model 

Figure  C-1  iUustrates  the  mechanical  motion 
ringdown  model  used.  A  single  degree-of-freedom 
mass  is  aUowed  to  osciUate  between  two  identical 
ideal  BeUeviUe  washers.  Mass  and  washer  stiff¬ 
ness  properties  are  those  defined  in  Appendix  A  of 
this  report.  An  initial  displacement,  y(0),  is 
imposed,  foUowed  by  free  osciUation  of  the  mass. 
Figure  C-2  iUustrates  a  typical  ringdown  displace¬ 
ment  history  for  an  idealized  BeUeviUe  washer 
subjected  to  moderate  initial  displacement. 
Observe  that  the  absolute  values  of  succeeding  dis¬ 
placement  relative  minima  and  maxima  decrease 
with  respect  to  the  absolute  values  of  time. 

Damping  Models 

Again,  referring  to  Figure  C-2,  the  envelope 
curve  tangent  to  the  displacement  history  relative 
maxima  is  defined  to  be  Y(t).  If,  for  a  simple  wash¬ 
er-mass  system,  Y(t)  is  Unear  with  respect  to  time, 
the  damping  behavior  is  termed  Coulomb  or  friction 
damping,  having  the  foUowing  description,  Yc(t): 

Y,(t)  =  y(0)  - 


An  exponentiaUy  decaying  definition,  Yv(t),  of 
the  envelope  defines  classical  viscous  damping  of 


a  simple  washer-mass-dashpot  system: 
where 


^  =  the  viscous  critical  damping  ratio 

T  =  the  natural  vibration  period. 


In  this  case,  z  is  easily  found  by  means  of  a 
technique  termed  logarithmic  decrement:  if  the 
natural  logarithm  of  the  ratio  of  adjacent  relative 
maxima  or  minima  of  y(t)  is  found. 


y(t)  ~ 
yCt+T. 


History  Calculations 

A  simple  FORTRAN  program,  listed  in  Table 
C-1,  numericaUy  determines  the  ringdown  behavior 
of  an  ideal  BeUeviUe  washer.  Force-displacement 
behavior  of  the  ideal  washer,  as  described  in 
Appendix  D  of  this  report,  is  defined  in  BLOCK 


a.  Figures  and  tables  are  found  in  a  section  beginning  on  page  5. 
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DATA  (respectively  listed  in  Tables  C-2  and  C-3  for 
the  opening  and  closing  washers).  Variables  input 
during  program  execution  are  initial  displacement, 
calculation  duration,  and  file  output  interval. 
ASCII  file  output  records  consist  of  the  following: 
time,  displacement,  velocity,  and  acceleration. 

Opening  Belleville  Washer 

Figures  C-3  and  C-4  illustrate  results  of  two 
ringdown  calculations  for  the  opening  Belleville 
washer  (12  washers).  The  first  figure  represents 
behavior  where  the  initial  displacement  is  set  to  be 
d3  (idealized  characterization  data.  Appendix  D) 
which  represents  normal  loading  of  the  washer. 
The  second  figure  represents  behavior  where  initial 
displacement  is  greater  than  d3  (overloaded  condi¬ 
tion),  It  is  seen  that  effective  damping  is  signifi¬ 
cantly  reduced  for  the  overload  case.  Also,  it  is 
seen  that  the  damping  envelope  is  essentially  lin¬ 
ear,  which  indicates  that  the  equivalent  simple 
washer-mass  damping  is  predominantly  frictional 
(Coulomb).  If  the  mass  of  the  moving  parts  of  the 
ETVE  is  taken  to  be  0.0324  Ibf-s2fin.,  then,  from 
Figure  C-3  (T  =  0.0023  s,  y(0)  =  0.0327  ia, 
Yc(0.0043  s)  =  0.0159  in.),  the  equivalent 
Coulomb  damping  force  for  normal  Belleville 
washer  ringdown  is 

Fen  =  540  Ibf. 

However,  from  Figure  C-4,  if  the  initial  dis¬ 
placement  of  the  Belleville  washer  is  increased  to 
0.050  in.,  damping  is  significantly  reduced,  hi  this 
case,  the  equivalent  damping  force  is  approximately 


Fco  =  3101bf. 

Logarithmic  decrement  of  the  latter  relative 
maxima  and  minima  from  Rgure  F-3  results  in  an 
average  equivalent  critical  viscous  damping  ratio: 

Cn  =  0.066  . 

Closing  Belleville  Washer 

Rgures  C-5  and  C-6  present  the  results  of  ring- 
down  calculations  for  the  ETVE  closing  Belleville 
washer  (4  washers).  Again,  the  first  figure  illus¬ 
trates  results  for  initial  displacement,  representing 
normal  loading  on  the  washer.  The  second  figure 
depicts  results  of  an  overload  condition  is.  In  this 
case,  the  ratio  of  the  overload  initial  displacement 
to  the  normal  initial  displacement  is  not  as  great  as 
this  ratio  was  for  the  opening  washer  comparison. 
We  see  that  the  reduction  of  equivalent  damping  is 
not  as  great  as  in  the  opening  washer  comparison. 
Coulomb  damping  force  calculations,  similar  to 
those  previously  performed,  indicate 

Fen  =  360  Ibf 

Fco  =  2401bf. 

Logarithmic  decrement  calculations  for  the 
normally  loaded  ringdown  case  result  in  an  aver¬ 
age  equivalent  viscous  critical  damping  ratio: 

Cn  =  0.055 . 
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Figure  C-4.  Overload  ringdown — opening  Belleville  washer. 
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Figure  C-5.  Normal  ringdown — closing  Belleville  washer. 
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Table  C-1.  Belleville  Spring  Ringdown  Program  Listing  - 


Program  damp  ---  Calculates  Belleville  spring  ringdown  history. 
An  initial  displacement  is  specified  for  a  mass  located 
between  two  identical  Belleville  springs  and  allowed  to 
C  oscillate  with  no  external  applied  force. 

C  Displacement,  velocity,  and  acceleration  are  calculated  and 

C  output  to  file  for  subsequent  analysis. 


C  Developed  in  support  of  ETVE  behavior  coupled  hydrodynamic  - 
C  mechanical  motion  prediction  code  for  U.S.  Army. 


C  J.  G.  Arendts,  INEL/EG&G  Idaho,  Dec.  1993. 


PROGRAM  damp 
IMPLICIT  REAL*8  (a-h,  o 
CHARACTER*12  z_fout 


COMMON  /mtn/  tmass , 

&  f_cdc , 

&  d_cur, 

COMMON  /bspr/  f l_bsc , 

&  dl_bsc , 

&  sl_bsc, 

&  sl_bso. 


z) 

gap_sc,  gap_so,  f_act, 
f_cdb,  s_csp,  c_csp, 
d_upd,  v_cur,  v_upd, 
f2_bsc,  f3_bsc,  fl_bso, 
d2_bsc,  d3_bsc,  dl_bso, 
s2_bsc,  s3_bsc,  s4_bsc, 
s2_bso,  s3_bso,  s4_bso 


c_bsp , 

a_cur,  a_upd 
f2_bso,  f3_bso, 
d2_bso,  d3_bso. 


C  ***  Calculate  Belleville  spring  characteristics  *** 


CALL  belspr 

C  ***  Open  output  file  *** 

PRINT  *,  ‘Output  File  Name:  ‘ 

READ  *,  z_fout 

OPEN  (UNIT=2,  FILE=z_fout) 

C  ***  Input  data  *** 

PRINT  *,  'Initial  Displacement:  ‘ 
READ  *,  d_upd 

PRINT  *,  ‘Calculation  Duration:  ‘ 
READ  *,  t_end 

PRINT  *,  ‘Output  Interval:  ‘ 

READ  *,  t_del 
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Table  C-1.  (continued) 


C  ***  Initialize  *** 

iu_flg  =  0 

C  ***  Time  loop  *** 

time  =  -  t_del 
10  time  =  time  +  t_del 
d_eur  =  d_upd 
v_cur  =  v_upd 
a_cur  =  a_upd 

WRITE  (2,FMT=’  (4E20.6)’ )  time,  d_cur,  v_cur,  a_cur 
CALL  motion  (time,  t_del ,  0.0,  iu_flg) 

C  ***  End  time  loop  ** 

IF  (time.LT.t_end)  GO  TO  10 
CLOSE  (2) 

STOP 

END 


Table  C-2.  Opening  Spring  Ringdown  Block  Data  Listing 


BLOCK  DATA 

IMPLICIT  REAL*8  (a-h,  o-z) 


COMMON  /mtn/ 

tmass , 

gap.sc. 

gap_so. 

f_act. 

f_cdc , 

f_cdb, 

s_csp. 

c_csp. 

c_bsp , 

t 

d_cur , 

d_upd. 

v_cur. 

v_upd, 

a_cur , 

a_upd 

COMMON  /bspr/ 

f l_bsc , 

f2_bsc , 

f 3_bsc , 

fl_bso,  f2_bso. 

f 3_bso , 

dl_bsc , 

d2_bsc , 

d3_bsc , 

dl_bso,  d2_bso, 

d3_bso. 

sl_bsc, 

s2_bsc , 

s3_bsc. 

s4_bsc. 

sl_bso, 

s2_bso, 

s3_bso. 

s4_bso 

DATA 

tmass, 

gap_sc. 

gap_so. 

f_act 

/  0.0324, 

0.0, 

,  l.OE-6, 

,  0.0 

/ 

DATA 

f  cdc, 

f_cdb. 

s_csp. 

c_csp , 

c_bsp 

/  0.0, 

,  0.0 

00 

1,  0.0, 

0.1 

0  / 

DATA 

d_cur , 

d_upd. 

v_cur, 

v_upd , 

a_cur 

,  a_upd 

O 

O 

,  0.0 

o 

o 

1,  0.0, 

0.0 

,  0.0/ 

C  ***  Initialize  Belleville  spring  force-displacement  data  here  *** 
DATA  fl_bsc,  f2_bsc,  f3_bsc,  fl_bso,  f2_bso,  f3_bso 
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/  8.5E+3,  5.5E+3,  2.1E+4,  8.5E+3,  5.5E+3,  2.1E+4/ 
dl_bsc,  d2_bsc,  d3_bsc,  dl_bso,  d2_bso,  d3_bso 
/  0.0280,  0.0285,  0.0328,  0.0280,  0.0285,  0.0328/ 


& 

DATA 

& 

Q  ***************************************************************** 

/ 
/ 


DATA  sl_bsc,  s2_bsc,  s3_bsc,  s4_bsc 

&  /  0.0,  0.0,  0.0,  0.0 

DATA  sl_bso,  s2_bso,  s3_bso,  s4_bso 

&  /  0.0,  0.0,  0.0,  0.0 

END 


Table  C-3.  Closing  spring  ringdown  block  data  list. 


BLOCK  DATA 

IMPLICIT  REAL*8  (a-h,  o-z) 


COMMON  /mtn/ 

tmass,  gap_sc,  gap_so. 

Oct, 

& 

f_cdc,  f_cdb,  s_csp. 

c_csp,  c_bsp. 

& 

d_cur,  d_upd,  v_cur. 

v_upd,  a_cur,  a_upd 

COMMON  /bspr/ 

fl_bsc,  f2_bsc,  f3_bsc. 

fl_bso,  f2_bso,  f3_bso. 

& 

dl_bsc,  d2_bsc,  d3_bsc. 

dl_bso,  d2_bso,  d3_bso. 

& 

sl_bsc,  s2_bsc,  s3_bsc. 

s4_bsc , 

& 

sl_bso,  s2_bso,  s3_bso. 

s4_bso 

DATA 

tmass,  gap_sc,  gap_so. 

f_act 

& 

/  0.0324,  0.0,  l.OE-6, 

0.0 

/ 

DATA 

f_cdc,  f_cdb,  s_csp. 

c_csp ,  c_bsp 

& 

/  0.0,  0.0,  8.0 

,  0.0,  0.0 

/ 

DATA 

d_cur,  d_upd,  v_cur. 

v_upd,  a_cur,  a_upd 

& 

/  0.0,  0.0,  0.0 

,  0.0,  0.0,  0.0/ 

Initialize  Belleville  spring  force-displacement  data  here  *** 

DATA 

fl_bsc,  f2_bsc,  f3_bsc. 

fl_bso,  f2_bso,  f3_bso 

& 

/  5.6E+3,  4.1E+3,  1.4E+4, 

5.6E+3,  4.1E+3,  1.4E+4/ 

DATA 

dl_bsc,  d2_bsc,  d3_bsc. 

dl_bso,  d2_bso,  d3_bso 

& 

/  0.0127,  0.0135,  0.0155, 

0.0127,  0.0135,  0.0155/ 

£  ***************************************************************** 


DATA  sl_bsc,  s2_bsc,  s3_bsc,  s4_bsc 

&  /  0.0,  0.0,  0.0,  0.0  / 

DATA  sl_bso,  s2_bso,  s3_bso,  s4_bso 

&  /  0.0,  0.0,  0.0,  0.0  / 

END 
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Appendix  D 


Mechanical  Dynamics  Subroutine 
for  the  Eaton  Throat-Vaive  Eiement  Prototype 


This  appendix  presents  the  analytical  develop¬ 
ment  of  the  mechanical  dynamic  model  represent¬ 
ing  the  Eaton  Throat-Valve  Element  (ETVE)  pro¬ 
totype  moving  parts.  The  ETVE’s  moving  parts 
(sleeve,  piston,  piston  rod,  and  Belleville  washers) 
are  analyzed,  from  which  the  equation  of  motion 
for  the  ETVE’s  configuration  is  derived.  The 
resulting  solution  of  the  nonlinear  equation  of 
motion  is  presented  as  explicit  time-difference 
equations.  The  paper  also  presents  FORTRAN 
implementation  of  the  motion  solution,  with  sup¬ 
porting  subroutines. 

Analytical  Model 

Since  all  of  the  moving  parts  of  the  ETVE  are 
rigidly  connected  and,  thus,  move  together,  we 
assume  a  single  degree-of-freedom  rigid  body 
dynamic  model.  Figure  D-1  represents  the  assumed 
model.  We  see  that  nonlinearities  are  introduced 
owing  to  the  slide  opaiing  and  closing  gaps  (Ggo  and 
Ggc.  respectively)  and  Coulomb  damping  friction 
forces  acting  on  the  sliding  mechanism  and 
Belleville  washers  (F^-g  and  F^j,,  respectively).  Stiff¬ 
ness  and  viscous  damping  of  the  coil  spring  and  slide 
mechanism  are  represented,  respectively,  as  Kg  and 
Cg.  Since  the  closing  and  opening  rebound  BeUeviUe 
washers  have  differing  stiffness  properties,  their 


respective  stiffoess  properties  are  represented  as 
and  K^q.  Viscous  damping  of  the  BeUeviUe  washers 
is  represented  as  C5.  RnaUy,  extemaUy  appUed 
force,  resulting  from  time-dependent  actuator  pres¬ 
sure  and  hydrodynamic  gas  pressure  (both  time-  and 
sUde  position-dependent),  is  represented  as  Fp(y,t). 

Belleville  Spring  Characterization 

Quasi-static  compression  tests  of  representa¬ 
tive  ETVE  BeUeviUe  spring  stacks  indicate  signif¬ 
icantly  different  loading  and  unloading  force-dis¬ 
placement  behavior.  Figure  D-2  represents  ideal¬ 
ized  quadrilinear  BeUeviUe  spring  behavior.  We 
see  that  the  entire  ideal  behavior  may  be  represent¬ 
ed  through  specification  of  three  force-displace¬ 
ment  coordinate  pairs,  (Fj,  dp,  where  i  =  1,  2,  3. 
The  four  stiffnesses  that  represent  the  overaU 
BeUeviUe  spring  behavior  are  then  foimd  as  the 
slopes  of  the  linear  curve  segments.  For  maximum 
compressive  forces  exceeding  F3,  unloading 
behavior  foUows  the  bounding  unloading  curves. 
However,  for  maximum  spring  compressive  forces 
less  than  F3,  initial  unloading  behavior  is  assumed 
to  foUow  a  curve  paraUel  to  the  bounding  initial 
unloading  curve  imtil  intersection  with  the  final 
bounding  unloading  curve,  at  which  time  unload¬ 
ing  behavior  foUows  the  latter  curve. 


Feb 


Figure  D-1 .  ETVE  mechanical  motion  model. 
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Force 


Equivalent  Spring  Force 

Since  the  Belleville  spring  force-displacement 
characterization  curves  do  not  necessarily  pass 
through  the  origin  of  the  force-global  displacement 
coordinate  system,  the  restoring  force  acting  on  the 
moving  mass  caused  by  the  Belleville  spring  com¬ 
pression  is  not  equal  to  the  current  global  displace¬ 
ment  multiplied  by  the  corresponding  stiffness. 
Figure  D-3  illustrates  the  Belleville  spring  restoring 
force  Fp  as  related  to  global  displacement,  y,  for  a 
typical  spring  characterization  curve.  We  see  that  Fj. 
is  equal  to  Ky  less  Fgg,  the  equivalent  spring  force. 
Note  that  when  neither  of  the  gaps  is  closed,  the 
equivalent  spring  force  is  zero,  and  the  coil  spring 
restoring  force  is  equal  to  Kjy. 

Equation  of  Motion 

Equating  restoring  forces  acting  on  the  single 
degree-of-freedom  mass  (inertial,  viscous  damping, 
spring,  and  Coulomb  friction)  with  the  total  disturb¬ 
ing  force  fiom  actuator  and  hydrodynamic  pressures 
yields  the  following  equation  of  motion: 

My-i-  C(y)  y  -i-  K(y)  y  -  (y)  -i-  F^  (y,y)  =  F^  (y,t) 

where 

M  =  the  total  mass  of  all  moving  valve  parts 

C(y)  =  the  viscous  damping  function 

c(,)=c..-e„<y<G„ 

C(y)  =  C;,y<-G^or  y>G,„ 

K(y)  =  the  stiffness  function 

K(y)=  K,,-0„<y<G,„ 

K(y)=K^^.y<-G„ 

K(y)=K^.y>G„ 

F  (y)  =  the  equivalent  spring  force  fimction 
Fjy)  =  0,  -G^^<y<G^^ 

Pes(y>  =  Feso’y>Gso 


FpCy.y)  =  the  Coulomb  damping  force  fimction 
Fc^y.  y)  =  0,  y  =  0 

F,(y,  y)  =  F^^,  <  y  <  G^^,  and  y  >  0 

F,(y.  y)  =  -F^.  <  y  <  G^^,  and  y  <  0 

F,(y,  y)  =  F^^,  +  F^j,,  y  >  G^^,  or  y  <  -G^^, 
andy  >0 

F,(y,  y)  =  -(F^3  +  FJ,  y  >  G^^,  or  y  <  -G^^, 
andy  >0 

Fp(y,t)  =  Fjj(y»t)  +  F^FCt)  is  the  pressure 
induced  load  history 

Fjj(y,t)  is  the  net  hydrodynamic  pressure 
load  history 

F^(t)  is  the  actuation  pressure  load  history. 

Linearized  Equation  of  Motion 

If  C(y),  K(y),  Fjy),  andF^(y,y)  are  constant 
over  some  small  time,  t,  0  <  t  < At,  and  are  equal  to 
C,  K,  Fg^,  respectively;  and  if 
F(t)  =  F^(y,t)  +  F,,-F^ 
then 

M^  +  C^  +  ky  =  Fa). 

Frequency  Description 

If  the  following  frequency  terms  are  defined. 

Undamped  natural  frequency,  ©  = 

Damping  factor  frequency,  P  =  C^M 
Damped  natural  fi'equency,  ©jj  =  (©^  -p2)l/2 
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then, 


df 


+  2p 


(ft 


2  1 
+  ®  y  = 


yO)  =e‘“|^(^iy^+iy^JsinWjl+y„cos(l)jJ 

4[>--"  (f 


-smco  .t  +  cos®  .t 

®d  "  " 


Closed  Solution 

The  above  equation  has  the  following  general 
solution: 

y(t)  =0"^* 


J. 

F(t)  sin®j  (t  -  t)  dx 
0 

=\® 

+yp® 

where  and  y^  are,  respectively,  displacement  and 
velocity  att-  0. 

Particular  integral 

Since  the  present  linearized  solution  is  to  be 
used  for  explicit  time  integration,  the  disturbing 
force  is  set  to  be  constant  over  the  interval,  0  <  t  <  x. 


^  (yo + W  yo 

d 


F(x)  =  Fo. 


dy0) 

dt 


-Pt 


®  p  . 

— y  +— y 
®/o  (d/o^ 
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d^(t) 
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=  e 
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LV  ®d  ®d 
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-  (®  +  2py^)  cos®jt 
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^  2  -pf 
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K 
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or. 


Explicit  Time  Integration 


y?  ®  I,  ® 

=  rr^  f  sin®.0-x)  dx  . 

M®d  Jo  “ 

Upon  integration, 

yp  ®  "  T  [  ^  ■  • 

Upon  combination  of  the  above  homogeneous 
and  particular  solutions  and  differentiation  of  the 
result,  displacement,  velocity,  and  acceleration  are, 
respectively. 


Let 


yo=y(t) 

yo=y(t) 

be  the  previously  calculated  displacement  and 
velocity  at  time,  t.  Then,  from  the  closed  solution, 
where  t  =  6t  , 


y  ^  X  -PSt 

y(t+  5t)  =  e 


+ 


+  y(t)  cos  (®j5t) 


+ 


F(t) 

K 


-P6t 


i. 

®d 


sin  (®j5t)  +  cos  (®j5t) 
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y(t  +  8t)  =  e 


-P6t 


(2  D  ^ 

®  /N  P 


<0^ 
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sin  (C0j8t) 


+  y(t)  cx)s  (C3 .50 


F(t)<ll>^'^'  .  , 
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sin  {(Oj80 


. .  2  -p6t 

F(0  (0  e 
K 


sin  (C0j50  +  cos  (co^jSO  j 


sured  actuator  pressures  (see  A^Jendix  A  of  this 
report).  Also  shown  is  the  ideal  actuator  history 
specified  for  the  mechanical  motion  model.  The 
time  and  force  values  referred  to  in  Figure  D-6  and 
used  in  the  model  are 
tai  =  0.05  sec 
ta2  =  0.50  sec 
Fai  =  1425  Ibf 
Fa2  =  28501bf. 

Overall  Model  Damping  Characterization 

We  obtained  the  system’s  damping  characteris¬ 
tics  by  running  several  algorithm  predictions,  vary¬ 
ing  different  parameters  until  the  axial  displace¬ 
ment  history  matched  the  actual  test  data  obtained 
from  the  tests  described  above  (see  Appendix  A  of 
this  report).  We  found  that  the  following  damping 
parameters  produced  a  reasonable  fit  to  the  test 
data: 


Experimental  Characterization 

To  determine  the  Belleville  washers’  ideal 
stiffness  behavior  and  dynamic  model  overall 
damping  characteristics,  results  from  quasi-static 
Belleville  spring  compression  tests  and  ETVE 
opening  tests  are  used. 

Belleville  Spring  Characterization 

Figures  D-4  and  D-5  illustrate  the  force-dis¬ 
placement  idealized  Belleville  spring  stiffness 
characteristic  curves  used  herein.  The  closing 
spring  is  a  stack  of  four  parallel  Belleville  washers, 
whereas  the  opening  spring  is  two  stacks  of  six 
washers  oriented  in  series.  Also  represented  are 
actual  average  test  data  (see  Appendix  B  of  this 
report).  The  idealized  curves  are  determined  by 
visual  best  fit  to  the  overall  test  data. 

Actuator  Force  History  Characterization 

I  used  the  results  fiom  a  recent  ETVE  opening 
test  that  used  2(XX)-psi  actuator  pressure  (quickly 
applied  and  having  a  long  duration)  but  no  driver 
gas  volume  or  pressure.  Figure  D-6  represents  a 
typical  actuator  force  history  derived  fi'om  mea- 


Slide  Coulomb  damping  force  =  100  Ibf 
Slide  viscous  damping  ratio  =  0.0 
Belleville  spring  Coulomb  damping  force  =  0.0 
Belleville  spring  viscous  damping  ratio  =  0.08. 

Figure  D-7  compares  the  calculated  ETVE  dis¬ 
placement  response  damping  data  with  the  test 
data. 

Subroutine  Descriptions 

Three  FORTRAN  subroutines  that  support  and 
perform  the  ETVE  motion  calculations  depicted 
above  are  described  here.  The  only  significant 
deviation  fiom  ANSI  FORTRAN-77  is  in  the  use 
of  underscored  variable  names  (which  will  be 
included  in  ANSI  FORTRAN-90).  Generic  names 
are  used  for  intrinsic  library  functions,  32  bit  real 
words  are  implicitly  defined,  and  constants  are 
defined  in  either  floating  point  or  “E”  fonnats. 
Three  named  common  blocks  are  used  for  global 
variable  storage  and  a  BLOCK  DATA  segment  is 
defined  for  global  variable  value  initialization. 
The  three  subroutines,  namely  “belspr,”  “fact,”  and 
“motion,”  are  described  below.  BLOCK  DATA  is 
also  listed. 
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ure  D-4.  Belleville  closing  spring  stiffness  characteristics. 
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Figure  D-5.  Belleville  opening  spring  stiffness  characteristics. 
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Rgure  D-6.  Actuator  force  histories. 


Figure  D-7.  ETVE  opening  transient  displacement  histories. 
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Subroutine  “belspr” 

The  purpose  of  this  subroutine  is  to  compute 
the  Belleville  washers’  stifihiess  from  the  force- 
displacement  characterization  values  specified  in 
BLOCK  DATA.  This  subroutine  is  only  called 
from  the  main  program  once  (prior  to  first  call  to 
“motion”).  The  source  list  of  this  subroutine  is 
given  in  Table  D-1,  and  use  of  the  variables  is 
described  as  follows  (see  Figure  D-2): 

fl_bsc  =  Fj  (closing  Belleville  spring) 
f2_bsc  ='  F2  (closing  Belleville  spring) 
f3_bsc  =  F3  (closing  Belleville  spring) 
fl_bso  =  Fj  (opening  Belleville  spring) 
f2J)SO  =  F2  (opening  Belleville  spring) 

O_bso  =  F3  (opening  Belleville  spring) 


dl_bsc  =  dj  (closing  Belleville  spring) 
d2_bsc  =  d2  (closing  Belleville  spring) 
d3_bsc  =  d3  (closing  Belleville  spring) 
dl_bso  =  dj  (opening  Belleville  spring) 
d2_bso  =  d2  (opening  Belleville  spring) 
d3_bso  =  d3  (opening  Belleville  spring) 
sl_bsc  - 

s4_bsc  =  calculated  stiffiiesses  (closing 
Belleville  spring) 

sl_bso  - 

s4_bso  =  calculated  stiffnesses  (opening 
Belleville  spring). 

Subroutine  “fact” 

This  subroutine  calculates  the  current  actuator 
force  acting  on  the  moving  ETVE  parts.  One  argu- 


Table  D-1.  Subroutine  “belspr”  List. _ _ 

C  Subroutine  belspr  --  Calculates  Belleville  spring  stiffness 
C  characteristics  from  force  displacement  values. 

C  Developed  for  inclusion  in  ETVE  behavior  coupled  hydrodynamic  - 
C  mechanical  motion  prediction  code  for  U.S.  Army. 

C  J.  G.  Arendts,  INEL/EG&G  Idaho,  Nov.  1993. 

SUBROUTINE  belspr 
IMPLICIT  REAL*8  (a-h,  o-z) 

COMMON  /bspr/fl_bsc,  f2_bsc,  f3_bsc,  fl_bso,  f2_bso,  f3_bso, 

&  dl_bsc,  d2_bsc,  d3_bsc,  dl_bso,  d2_bso,  d3_bso, 

&  sl_bsc,  s2_bsc,  s3_bsc,  s4_bsc, 

&  sl_bso,  s2_bso,  s3_bso,  s4_bso 

C  ***  Belleville  closing  spring  characteristics  *** 

sl_bsc  =  fl_bsc/dl_bsc 

s2_bsc  =  f2_bsc/d2_bsc 

s3_bsc  =  (f3_bsc-fl_bsc)/(d3_bsc-dl_bsc) 

s4_bsc  =  (f3_bsc-f2_bsc)/(d3_bsc-d2_bsc) 

C  ***  Belleville  opening  spring  characteristics  *** 

sl_bso  =  fl_bso/dl_bso 

s2_bso  =  f2_bso/d2_bso 

s3_bso  =  (f3_bso-fl_bso)/(d3_bso-dl_bso) 

s4_bso  =  (f3_bso-f2_bso)/(d3_bso-d2_bso) 

RETURN 

END 
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ment  is  passed,  the  current  solution  time.  This  sub¬ 
routine  must  be  called  prior  to  each  call  of 
“motion.”  The  source  list  of  this  subroutine  is 
given  in  Table  D-2,  and  use  of  the  variables  is 
described  as  follows  (see  Figure  D-6): 

fl_act  =  Fgj 
f2_act  =  F32 
tl_act  =  tai 
t2_act  =  tg2 

f_acf  =  calculated  current  actuator  force. 

Subroutine  “motion” 

This  subroutine,  listed  in  Table  D-3,  calculates 
the  incremental  motion  of  the  ETVE  movable  parts 
in  terms  of  the  variables:  displacement,  velocity, 
and  acceleration.  Arguments  passed  to  the  subrou¬ 
tine  are  current  calculation  time,  current  time 
increment,  net  fluid  dynamic  force  acting  on  the 
moving  mass,  and  an  integer  flag,  iu_flg,  which 
indicates  whether  a  Belleville  spring  is  loading 
(iu_flg  =  0)  or  unloading  (iu_flg  =  1).  Note  that 
iu_flg  is  not  initialized  in  BLOCK  DATA;  it  must 
be  initialized  to  zero  early  in  the  main  program. 

The  subroutine  first  determines  if  either  of  the 
the  Belleville  washers  are  compressed  upon  sub¬ 
routine  call  or,  if  not,  will  a  Belleville  spring 
become  compressed  at  the  end  of  the  time  step.  If 
either  of  the  above  cases  are  true,  the  motion  cal- 
culational  time  step  is  reduced  to  1.0  x  (10)'^  sec 


and  the  incremental  motion  calculations  are  sum¬ 
marily  performed  imtil  die  argument  passed  time 
increment  is  attained.  The  reason  for  this  small 
incremental  motion  calculation  is  due  to  the  very 
large  changes  of  the  motion  parameters  that  occur 
in  the  transition  between  relatively  free  motion 
(neither  BelleviUe  spring  compressed)  and  very 
stiff  resistance  when  either  of  the  Belleville  wash¬ 
ers  are  compressed. 

The  following  section  of  the  subroutine  con¬ 
sists  of  logic  blocks  wherein  the  applicable  motion 
parameters  are  specified.  Finally,  the  updated 
motion  variables  are  computed.  It  is  important  to 
know  that  the  subroutine  does  not  update  the  cur¬ 
rent  motion  variables;  rather,  the  updated  (end  of 
time  increment)  variables  are  re-computed  in  the 
subroutine.  It  is  the  responsibility  of  the  calling 
program  to  update  the  current  variables,  either 
before  “motion”  is  called,  or  after. 

Parametric  variables  (mass,  stiffnesses,  and 
damping  values),  as  well  as  motion  variables  and 
forces,  are  described  in  the  subroutine  listing  and 
will  not  be  repeated  here. 

Block  Data 

The  purpose  of  BLOCK  DATA,  listed  in  Table 
D-4,  is  to  initialize  all  global  variables.  General 
motion  and  spring  parameters  ate  specified  here. 
Changes  to  these  parameter  values  are  easily  made 
in  the  source  module  (with  recompilation). 


Table  D-2.  Subroutine  “fact”  List. 


C  Subroutine  fact  ---  Calculates  actuator  force  as  a  function 
C  of  current  time. 

C  Developed  for  inclusion  in  ETVE  behavior  coupled  hydrodynamic  - 
C  mechanical  motion  prediction  code  for  U.S.  Army. 

C  J.  G.  Arendts,  INEL/EG&G  Idaho,  Nov.  1993. 
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Table  D-2.  (continued) 


SUBROUTINE  fact  (time) 

IMPLICIT  REAL*8  (a-h,  o-z) 

COMMON  /fac/  f_act,  fl_act,  f2_act,  tl_act,  tZ_act 

C  ***  Current  actuator  force  *** 

IF  (  time.LT.tl_act  )  THEN 

f_act  =  fl_act*tiine/tl_act 
ELSE  IF  (  time.GT.t2_act  )  THEN 
f_act  =  f2_act 

ELSE 

f_act  =  fl_act+(f2_act-fl_act)*(time-tl_act)/(t2_act-tl_act) 
ENDIF 

RETURN 

END 


Table  D-3.  Subroutine  “motion"  list. 


C  Subroutine  motion  ---  Calculates  and  updates  displacement, 

C  velocity,  and  acceleration  of  single  degree-of-freedom, 

C  damped,  and  forced  (constant  over  time  step) 

C  dynamic  mechanical  system.  A  closed-form  solution  of 

C  equation-of-motion  is  utilized  for  explicit  time  integration. 

C  Developed  for  inclusion  in  ETVE  behavior  coupled  hydrodynamic  - 
C  mechanical  motion  prediction  code  for  U.S.  Army. 

C  J.  G.  Arendts,  INEL/EG&G  Idaho,  Nov.  1993. 


C 

C 

C 

C 

C 


Argument  passed  variables: 

time  =  current  absolute  time  (sec) 

t_del  =  time  increment  (sec) 

f_hyd  =  current  net  hydrodynamic  force  (Ibf) 

iu_flg  =  Belleville  spring  unlaoding  flag  (=  1  if  unloading) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 


Named  common 
tmass 
gap_sc  = 
gap_so  = 
f_act 
f_cdc 
f_cdb 
s_csp 
c_csp 
c_bsp 
d_cur 
d_upd 
v_cur 
v_upd 


(motion)  variables: 

total  mass  of  moving  parts  (lbf-sec**2/in) 

slide  closing  gap  (ref.  to  undisplaced  position  -  in.) 

slide  opening  gap  (ref.  to  undisplaced  position  -  in.) 

current  actuator  force  (Ibf) 

coil  spring  Coulomb  damping  force  (Ibf) 

Belleville  spring  (both)  Coulomb  damping  force  (Ibf) 
coil  spring  stiffness  (Ibf /in) 

coil  spring  viscous  critical  damping  ratio  (dim. -less) 

Belleville  spring  viscous  critical  damping  ratio  (“) 

current  displacement  (in.) 

updated  displacement  (in.) 

current  velocity  (in/sec) 

updated  velocity  (in/sec) 
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Table  D-3.  (continued) 


C  a_cur  =  current  acceleration  (in/sec**2) 

C  a_upd  =  updated  acceleration  (in/sec**2) 

C  Calculated  volatile  variables: 

C  fre_u  =  undamped  natural  frequency  (rad/sec) 

C  fre_d  =  damped  natural  frequency  (rad/sec) 

C  beta  =  damping  factor  frequency  (rad/sec) 

C  f_cd  =  velocity  dependent  Coulomb  damping  force  (Ibf) 

C  f_es  =  effective  gap-closure  spring  force  (Ibf) 

SUBROUTINE  motion  (time,  t_del,  f_hyd,  iu_flg) 

IMPLICIT  REAL*8  (a-h,  o-z) 

COMMON  /mtn/  tmass,  gap_sc,  gap_so, 

&  f_cdc,  f_cdb,  s_csp,  c_csp,  c_bsp, 

&  d_cur,  d_upd,  v_cur,  v_upd,  a_cur,  a_upd 

COMMON  /fac/  f_act,  fl_act,  f2_act,  tl_act,  t2_act 
COMMON  /bspr/  fl_bsc,  f2_bsc,  f3_bsc,  fl_bso,  f2_bso,  f3_bso, 

&  dl_bsc,  d2_bsc,  d3_bsc,  dl_bso,  d2_bso,  d3_bso, 

&  sl_bsc,  s2_bsc,  s3_bsc,  s4_bsc, 

&  sl_bso,  s2_bso,  s3_bso,  s4_bso 

SAVE  f2_bst,  d2_bst 

C  ***Net  external  force  *** 

f_ex  =  f_hyd  +  f_act 

C  ***  Initialize  *** 

tt  =  t_del 
tt_sum  =0.0 
dl  =  d_cur 
vl  =  v_cur 
g_sc  =  -gap_sc 
g_so  =  gap_so 

C  ***Set  gap  boundary  and  Belleville  spring  compression  flags  *** 
i_flg  =  0 

IF  (  d_cur.GT.g_so  .OR.  d_cur.LT.g_sc  )  i_flg  =  I 
10  j_flg  =  0 

IF  (  i_flg.EQ.l  .AND.  t_del.GT.1.0E-07  )  tt  =  l.OE-07 

C  ***  Current  stiffness,  damping,  and  effective  forces  *** 

20  IF  (  dl.GT.g_so  )  THEN 

C  ***  Belleville  opening  spring  is  active  *** 

beta  =  c_bsp 
f_cd  =  f_cdc  +  f_cdb 
d_rel  =  dl  -  g_so 
f_so  =  s_csp*g_so 
IF  (  vl.GE.0.0  )  THEN 
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Table  D-3.  (continued) 


C  ***The  spring  is  loading  *** 
iu_flg  =  0 

IF  (  d_rel.GT.d3_bso  )  THEN 
stif  =  s4_bso 

f_es  =  s4_bso*{g_so+d3_bso)  -  f_so  -  f3_bso 
ELSE  IF  (  d_rel.LT.dl_bso  )  THEN 
stif  =  sl_bso 
f_es  =  sl_bso*g_so  -  f_so 
ELSE 

stif  =  s3_bso 

f_es  =  s3_bso*(g_so+dl_bso)  -  f_so  -  fl_bso 
ENDIF 

ELSE 

C  ***The  spring  is  unloading  *** 

IF  (  iu_flg.Ea.O  )  THEN 
C  ***  Do  this  only  for  the  initial  unload  step 
iu_flg  =1 

IF  (  d_rel.LT.dl_bso  )  THEN 

d2_bst  =  d_rel*(s4_bso-sl_bso)/(s4_bso-s2_bso) 
f2_bst  “  s2_bso*d2_bst 
ELSE  IF  (  d_rel.LT.d3_bso  )  THEN 

d2_bst  =  d_rel*(s4_bso-s3_bso)  +  dl_bso*s3_bso 
d2_bst  =  (d2_bst-fl_bso)/(s4_bso-s2_bso) 
f2_bst  =  s2_bso*d2_bst 
ELSE 

d2_bst  =  d2_bso 
f2_bst  =  f2_bso 
ENDIF 
ENDIF 

IF  (  d_rel.GT.d2_bst  )  THEN 
stif  =  s4_bso 

f_es  =  s4_bso*(g_so+d2_bst)  -  f_so  -  f2_bst 
ELSE 

stif  =  s2_bso 
f_es  =  s2_bso*g_so  -  f_so 
ENDIF 
ENDIF 

ELSE  IF  (  dl.LT.g_sc  )  THEN 

C  *** Belleville  closing  spring  is  active  *** 

beta  =  c_bsp 
f_cd  =  f_cdc  +  f_cdb 
d_rel  =  -dl  -  gap_sc 
f_sc  =  s_csp*gap_sc 
IF  (  vl.LE.0.0  )  THEN 
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Table  D-3.  (continued) 


C  ***The  spring  is  loading  *** 
iu_flg  =  0 

IF  (  d_re1.GT.d3_bsc  )  THEN 
stif  =  s4_bsc 

f_es  =  -s4_bsc*(gap_sc+d3_bsc)  +  f_sc  +  f3_bsc 
ELSE  IF  (  d_rel.LT.dl_bsc  )  THEN 
stif  =  sl_bsc 

f_es  =  -sl_bsc*gap_sc  +  f_sc 
ELSE 

stif  =  s3_bsc 

f_es  =  -s3_bsc*(gap_sc+dl_bsc)  +  f_sc  +  fl_bsc 
ENDIF 
ELSE 

C  ***The  spring  is  unloading  *** 

IF  (  iu_flg.EQ.O  )  THEN 
iu_flg  =  1 

IF  (  d_re1.LT.dl_bsc  )  THEN 

d2_bst  =  d_rel*(s4_bsc-sl_bsc)/(s4_bsc-s2_bsc) 
f2_bst  =  s2_bsc*d2_bst 
ELSE  IF  (  d_re1.LT.d3_bsc  )  THEN 

d2_bst  =  d_rel*(s4_bsc-s3_bsc)  +  dl_bsc*s3_bsc 
d2_bst  =  (d2_bst-fl_bsc)/(s4_bsc-s2_bsc) 
f2_bst  =  s2_bsc*d2_bst 
ELSE 

d2_bst  =  d2_bsc 
f2_bst  =  f2_bsc 
ENDIF 
ENDIF 

IF  {  d_rel.GT.d2_bst  )  THEN 
stif  =  s4_bsc 

f_es  =  -s4_bsc*(gap_sc+d2_bst)  +  f_sc  +  f2_bst 
ELSE 

stif  “  s2_bsc 

f_es  =  -s2_bsc*gap_sc  +  f_sc 
ENDIF 
ENDIF 
ELSE 

C  ***0nly  the  coil  spring  is  active  *** 

stif  =  s_csp 
beta  =  c_csp 
f_es  =0.0 
f_cd  =  f_cdc 
ENDIF 

C  ***  Displacement  and  velocity  dependent  properties  and  forces  *** 

fre_u  =  SQRT(stif/tmass) 

beta  =  beta*fre_u 

fre_d  =  SQRT(fre_u**2-beta**2) 

IF  {  vl.GE.0.0  )  f_cd  =  -f_cd 
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Table  D-3.  (continued) 


IF  (  time.LT.tl_act  )  f_cd  =  f_cd*f_act/fl_act 
f_tot  =  f_ex  +  f_cd  +  f_es 

C  ***  Exponential  and  trig  factors  *** 

fac_e  =  EXP(-beta*tt) 
fac_s  =  SIN(fre_d*tt) 
fac_c  =  COS(fre_d*tt) 

C  ***  Updated  or  incremental  displacement  *** 

d2  =  fac_e  *{fac_s  *(dl*beta/fre_d  +  vl/fre_d)  +  fac_c*dl) 
d2  =  d2  +  f_tot*(L0  -  fac_e  *(fac_s*beta/fre_d  +  fac_c))/stif 

C  ***  Check  gap  boundary  crossings  *** 

IF  {  i.flg.EQ.O  )  THEN 

IF  (  d2.GT.g_so  .OR.  d2.LT.g_sc  )  THEN 
i_flg  =  1 
j_flg  =  1 

END  IF 
ENDIF 

IF  {  j.flg.EO.l  )  GO  TO  10 

C  ***  Updated  or  incremental  velocity  *** 

v2  =  -fac_s  *(dl*fre_u**2/fre_d  +  vl*beta/fre_d) 
v2  =  fac_e  *(v2  +  fac_c*vl) 

v2  =  v2  +  f_tot*fac_e*fac_s*fre_u**2/stif/fre_d 

C  ***  Updated  or  incremental  acceleration  *** 

a2  =  dl*beta*fre_u**2/fre_d  -  vl  *(fre_d**2  -  beta**2)/fre_d 
a2  =  fac_e  *(a2*fac_s  -  fac_c  *(dl*fre_u**2  +  2.0*vl*beta)) 
a2  =  a2  +  f_tot*fre_u**2*fac_e  *(-fac_s*beta/fre_d  +  fac_c)/stif 

C  ***  Check  for  end  of  step  *** 

dl  =  d2 
vl  =  v2 

tt_sum  =  tt_sum  +  tt 

IF  (  tt_sum.LT.t_del  )  GO  TO  20 

C  ***  Finish-up  *** 


d_upd  =  d2 
v_upd  =  v2 
a_upd  =  a2 

RETURN 

END 
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Table  D-4.  Block  Data  Listing. 


BLOCK  DATA 

IMPLICIT  REAL*8  (a-h,  o-z) 
COMMON  /mtn/  tmass,  gap_sc, 

gap_so. 

& 

f_cdc,  f_cdb. 

s_csp,  c_csp,  c_bsp. 

& 

d_cur ,  d_upd , 

v_cur,  v_upd,  a_cur. 

a_upd 

COMMON  /fac/ 

f_act,  fl_act, 

f2_act,  tl_act,  t2_act 

COMMON  /bspr/ 

fl_bsc,  f2_bsc. 

f3_bsc,  fl_bso,  f2_bso. 

f3_bso. 

& 

dl_bsc,  d2_bsc, 

d3_bsc,  dl_bso,  d2_bso. 

d3_bso. 

& 

sl_bsc,  s2_bsc. 

s3_bsc,  s4_bsc. 

& 

sl_bso,  s2_bso, 

s3_bso,  s4_bso 

C  ****  Initialize  general  motion  parameters  here  ******************* 

DATA 

tmass,  gap_sc, 

gap_so 

& 

/  0.0324,  0.0, 

1.12 

/ 

DATA 

f_cdc ,  f_cdb , 

s_csp,  c_csp,  c_bsp 

& 

/  100.0,  0.0, 

8.0,  0.0,  0.08 

/ 

DATA 

d_cur,  d_upd. 

v_cur,  v_upd,  a_cur. 

a_upd 

& 

/  0.0,  0.0, 

0.0,  0.0,  0.0, 

,  0.0/ 

C  ****  Initialize  actuator  force  data  here  ************************* 


DATA  f_act,  fl_act,  f2_act,  tl_act,  t2_act 

&  /  0.0,  1425.0,  2850.0,  0.05,  0.5  / 

C  ****  Initialize  Belleville  spring  force -displacement  data  here  *** 


DATA  fl_bsc,  f2_bsc,  f3_bsc,  fl_bso,  f2_bso,  f3_bso 

&  /  5.6E+3,  4.1E+3,  1.4E+4,  8.5E+3,  5.5E+3,  2.1E+4/ 

DATA  dl_bsc,  d2_bsc,  d3_bsc,  dl_bso,  d2_bso,  d3_bso 

&  /  0.0127,  0.0135,  0.0155,  0.0280,  0.0285,  0.0327/ 

Q  Tk-***************************************************************** 


DATA  sl_bsc,  s2_bsc,  s3_bsc,  s4_bsc 

&  /  0.0,  0.0,  0.0,  0.0  / 

DATA  sl_bso,  s2_bso,  s3_bso,  s4_bso 

&  /  0.0,  0.0,  0.0,  0.0  / 


END 
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Intentionally  left  blank. 
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Intentionally  left  blank. 
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APPENDIX  B 


GAS  DYNAMIC  SIMULATIONS 


This  appendix  describes  the  gas  dynamics  of  the  Eaton  Throat- Valve  Element  Prototype 
(ETVE)  with  compressible,  inviscid,  nonheat  conducting,  time-dependent  Euler  equations  formu¬ 
lated  for  axisymmetiic  coordinate  systems.  The  code  is  formulated  with  a  variable  thickness  or 
depth,  denoted  by  A  for  each  mesh  cell.  In  two-dimensional  computations,  the  volume  of  a  cell  of 
width  5x  and  height  6x  is  A6x^,  giving  a  quasi-three-dimensional  effect.  The  physical  assumption 
underlying  the  variable  depth  description  is  that  the  flow  variables  are  depth-averaged  in  the  depth 
direction  to  give  mean  values  with  variation  only  in  the  x  and  y  dimensions,  or,  equivalently,  that 
the  flow  variables  are  independent  of  displacement  in  the  depth  direction.  The  variable-depth 
equations  wwe  used  here  because  they  greatly  increase  the  geometrical  flexibility  of  the  descrip¬ 
tion  without  correspondingly  increasing  its  complexity.  The  variable  A  may  be  used  to  represent 
flow  through  a  two-dimensional  duct  of  gradually  varying  thickness.  Axisymmetiic  coordinates 
are  generated  by  having  A  increase  linearly  with  distance  from  the  axis.  The  use  of  zero  values  of 
A  in  selected  mesh  cells  provided  a  convenient  means  of  including  stationary  internal  obstacles 
(valve  internal  parts)  in  the  flow  region.  In  particular,  a  cell  with  A  =  0  will  allow  no  flow  across 
its  boundaries.  As  presently  coded,  the  A-quantity  must  be  a  constant  in  time. 

The  valve  model  uses  an  Eulerian  description  for  the  gas  dynamics  (with  the  stationary 
valve  internals  represented  with  zero  depth,  A  =  0),  while  the  moving  valve  parts  are  treated  as  a 
moving  material  interface.  Various  problems  can  arise  in  the  treatment  of  material  interfaces:  (a) 
their  discrete  representation,  (b)  their  evolution  in  time,  and  (c)  the  manner  in  which  boundary 
conditions  are  imposed  on  them. 

In  an  Eulerian  representation,  the  computational  grid  remains  fixed,  and  the  identity  of 
individual  fluid  elements  is  not  maintained.  It  is  then  necessary  to  compute  the  flow  of  fluid 
through  the  mesh.  This  flow,  or  convective  flux  calculation,  necessarily  requires  an  averaging  of 
the  flow  properties  of  all  fluid  elements  that  end  up  in  a  given  mesh  cell  after  some  period  of  time. 
This  convective  averaging  results  in  a  smoothing  of  all  variations  in  flow  quantities,  and,  in  par¬ 
ticular,  a  smearing  of  surfaces  of  discontinuity  such  as  material  interfaces.  In  order  to  overcome 
this  loss  of  boundary  resolution,  an  adaptation  of  the  VOF  (fractional  volume  of  fluid)  technique^ 
was  incorporated  that  recognizes  a  discontinuity  and  avoids  averaging  across  it. 

With  the  modified  VOF  technique  used  here,  a  function  F(x,y,t)  is  defined  whose  value  is 
zero  at  any  point  occupied  by  gas  and  unity  elsewhere.  When  averaged  over  the  cells  of  a  com¬ 
putational  mesh,  the  average  value  of  F  in  a  cell  is  equal  to  the  fractional  volume  of  the  cell  occu¬ 
pied  by  a  solid.  In  particular,  a  zero  value  of  F  corresponds  to  a  cell  full  of  fluid,  whereas  a  unit 
value  indicates  that  the  cell  contains  no  fluid  and  is  occupied  totally  by  solid,  moveable  valve 
components.  Cells  with  F  values  between  zero  and  one  contain  a  material  interface. 

The  F  function  is  also  used  to  define  where  fluid  is  located  in  an  interface  cell.  The  nor- 
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mal  direction  to  the  boundary  lies  in  the  direction  in  which  the  value  of  F  changes  most  rapidly. 
Because  F  is  a  step  function,  its  derivative  is  computed  in  a  special  way,  using  a  donor-acceptor 
method.^ 

This  technique  was  used  because  it  provides  a  means  of  following  the  gas-solid  interfaces 
of  the  valve  through  an  Eulerian  mesh  of  stationary  cells,  and  because  it  follows  regions  rather 
than  boundaries;  it  avoids  the  logic  problems  associated  with  explicitly  tracking  interface  sur¬ 
faces. 

Equations  and  Constitutive  Reiations 

llie- Euler  equations  describing  the  dynamics  of  an  inviscid  gas  have  been  cast  in  many 
forms.  In  the  present  case,  is  chosen  the  dependent  variables  are  chosen  to  be  the  mass  density 
p,  the  gas  velocities  u  and  v,  and  the  internal  eneigy  per  unit  volume  E.  In  terms  of  the  chosen 
dependent  variables,  the  basic  two-dimensional  Euler  equations  used  in  the  code  are  the  continu¬ 
ity  equation. 


dp  .  1  r  d(puA)  d(pvA)  I 
dt  A  L  dx  dy  J 

(1) 

the  momentum  equations. 

du  du  du  1  dp 

dt  ^  dx  ^  dy  p  dx 

(2) 

dv  dv  dv  1  dp 

dt  ^  dx  ^  dy  ~  p  dy 

(3) 

and  the  internal  energy  equation 

dE  .  1  raCEuA)  .  d(EvA)l  (Y-l)Erd(uA)  d(vA)l 

dt  ^  A  L  dx  '  dy  J  A  L  dx  dy  J 

(4) 

In  these  equations,  the  independent  variables  are  time  t  and  coordinates  x  and  y.  The  grav¬ 
itational  force  terms  have  been  neglected,  as  have  the  heat  conduction  term,  assuming  they  are 
insignificant.  The  viscous  terms  have  been  dropped  to  simplify  the  model,  with  the  assumption 
that  reasonable  approximations  could  be  attained  without  their  inclusion.  The  quantity  A  is  the 
time-independent  thickness  or  depth  of  the  flow  channel.  Thus,  the  volume  associated  with  an 
area  5x  wide  and  5y  high  is  A5x8y.  Suitably  defined  A  values  were  used  to  represent  cylindrical 
coordinates  (A  =  y,  the  circumferential  area  per  unit  azimuthal  angle)  and  stationary  obstacle 
regions  (A  =  0). 

To  complete  these  equations,  a  constitutive  relations  relating  pressure,  p,  and  internal 
energy,  E,  must  be  defined.  For  this  case,  assume  the  gas  to  obey  the  perfect  gas  relationship 
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P  =  (Y-1)E 


(5) 


where  y  is  the  ratio  of  specific  heats  at  constant  pressiffe  and  at  constant  volume. 


The  time  dependence  of  the  solid  volume  fraction  function  F  is  governed  by  the  equation. 


dF  dF 
dt  "s  ax 


=  0 


(6) 


where  Us(x,t)  is  the  velocity  of  the  moveable  solid  valve  components  in  the  x-coordinate  direction 
since,  or  the  valve  considered,  a  single  degree  of  freedom  in  the  axial  direction  exists.  The  solid 
volume  fraction  function  F(x,y,t)  takes  on  point  values  of  0  or  1-  0  if  the  point  (x,y)  is  occupied 
by  gas  at  time  t,  and  1  if  the  point  is  occupied  by  a  solid  at  time  t.  In  the  numerical  scheme,  the 
mesh  value  of  F  is  an  average  of  the  point  values  of  F  within  a  computational  cell,  which  repre¬ 
sents  the  volume  fraction  of  that  cell  occupied  by  the  solid.  Obviously,  the  volume  fraction  of 
that  cell  occupied  by  gas  is  1-F. 


Dynamical  Coupling 


To  represent  the  valve  opening  dynamics,  the  moving  parts  of  the  valve  must  be  coupled 
to  the  subsequent  fluid  d5mamics  of  the  gas  as  it  flows  through  the  valve.  The  gas  surrounding 
the  moving  valve  parts  exerts  a  force  (in  fact,  the  dominant  force)  on  the  valve  parts  to  produce 
their  motion.  The  valve  sleeve  motion,  in  turn,  influences  the  gas  directly  with  its  motion  and 
indirectly  by  opening  or  closing  the  port  holes  to  alter  the  gas  flow  dynamics.  The  valve  sleeve 
motion  is  also  influenced  by  other  components  and  phenomena  such  as  the  centering  spring, 
sleeve  friction,  Belleville  spring  loading  and  unloading,  and  the  actuating  hydraulics.  The  gas 
dynamic  system  and  valve  slide  train  dynamics  are  solved  together,  in  a  fully  coupled  manner. 
Information  is  exchanged  explicitly  between  the  gas  dynamic  and  sleeve  dynamic  models  every 
gas  dynamic  time  step.  The  sleeve  dynamic  model,  however,  may  subcycle  its  time  step  to  val¬ 
ues  smaller  than  a  gas  dynamic  time  step  when  appropriate  (such  as  near  the  reversal  times,  etc.). 

Numerical  Methods 


This  section  briefly  describes  the  numerical  technique  used  to  integrate  the  gas  dynanoics 
and  advecting  solid  fraction  function.  The  numerical  technique  used  to  integrate  the  valve  slide 
train  dynamics  is  included  in  Appendix  D. 

For  simplicity,  the  solution  domain  is  divided  into  a  two-dimensional  array  of  rectangular 
cells  or  control  volumes  each  of  dimension  5x  by  5y.  To  each  of  these  cells,  a  cell-averaged  den¬ 
sity  pg,  pressure  pg,  internal  energy  Ejj,  depth  Ag,  and  solid  volume  fraction  Fg  are  assigned. 
A  staggered  mesh  is  used  wherein  average  velocities,  ug  and  vg  (respectively  for  the  x-  and  y- 
coordinate  directions),  are  assigned  at  the  cell  edges.  For  cell-centered  quantities,  the  g  subscripts 
indicate  position  (i5x,  j6y);  for  edge-defined  u-velocities,  the  g  subscripts  indicate  position  ((i+ 
V2)5x,  j5y);  for  the  edge-defined  v-velocities,  these  subscripts  indicate  position  (i5x,  (j-n  V2)5y)- 
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Since  we  expected  rapid  valve  motions  to  produce  significant  gas  dynamical  waves,  we 
used  an  explicit  Euler  method  to  integrate  the  conservation  and  F-advection  equations  in  time. 
The  convective  terms  are  differenced  spatially  with  a  fixed  variable  weighted  upwind  differenc¬ 
ing  formula.  By  setting  the  parameter  a  equal  to  unity,  full  upwind  differencing  results;  where¬ 
as,  setting  a  to  zero  produces  centered  differencing.  Setting  a  to  a  value  between  zero  and  one 
gives  a  linearly  weighted  combination  of  upwind  and  centered  spatial  differencing. 

The  first  calculation  in  the  solution  sequence  involves  computing  the  forces  on  the  move- 
able  valve  parts  owing  to  the  pressures  in  the  gas.  This  is  accomplished  by  assigning  to  each  cell 
an  index,  the  value  of  which  depends  on  whether  the  cell  is  an  interface  cell.  If  it  is  an  interface 
cell,  its  index  value  depends  on  the  orientation  of  the  interface.  To  compute  the  gas  dynamic  force 
on  the  sleeve,  an  individual  cell’s  contributions  are  summed  with  appropriated  coefficients  deter¬ 
mined  by  the  cell’s  interface  index.  For  example,  if  a  cell’s  index  does  not  indicate  the  cell  to  be 
an  interface  cell,  the  coefficient  for  this  cell  is  zero,  i.e.,  it  contributes  no  force  to  the  sleeve.  On 
the  other  hand,  if  a  cell’s  index  indicates  it  has  an  interface  with  positive  outward  normal  in  the 
positive  x-direction,  the  cell’s  coefficient  is  -8y,  and  its  contribution  to  the  total  force  on  the  sleeve 
is  -pj  j8y.  Interface  cells  with  orientation  in  the  negative  x-direction  contribute  +Pi  jSy  to  the  total 
force!  Other  cell  interface  orientations  contribute  nothing  to  the  total  sleeve  force  since  a  single 
degree  of  freedom  (in  the  x-direction)  was  assumed  for  the  valve  sleeve  spool. 


Next,  the  valve  dynamics  are  integrated  for  a  time  step,  or,  if  warranted,  several  sub-time  steps 
to  get  the  new  displacement,  velocity,  and  acceleration  of  the  moveable  valve  parts.  With  the  new 
valve  velocity,  the  F-function  is  advected  for  a  time  step  to  get  the  new  gas  dynaimc  flow  con¬ 
figuration.  With  the  flow  configuration,  the  gas  dynamic  solution  is  ready  to  be  updated 
(advanced  one  time  step).  Hrst,  the  u-  and  v- velocity  equations.  Equations  (2)  and  (3),  are  updat¬ 
ed  as 


=  u.  .  -1-  8t 

ij 


-FUX-FUY 


(7) 


V.  .  =  V.  .  4-  St  -r— 7 - r 
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FVX  = 


FVY  = 


m  + "i.j)  (vi,j  -  Vi,j’ 

“  I  Vi.j+1  +  “ijrt  +  “i-i.i  +  "i.i  I  (Vi.j  -  \i  +  '’mjM 
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In  these  discretized  equations  (and  in  all  subsequent  discretized  equations),  the  superscripts  indi¬ 
cate  time  level,  e.g.,  uy"'*'!  indicates  the  u-velocity  at 

n+1 

time  =  ^  8t 
i=l 

where  Sq  is  the  time  step  size,  and  the  absence  of  a  superscript  implies  that  the  variable  is  evalu¬ 
ated  at  the  current  time  level, 

n 

time  =  ^  6l  • 
i=l 


Using  the  updated  velocity  field,  the  mass  density  is  next  updated  to  the  new  time  level  as 
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ABT  = 
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i.j  i.J+1 

A.  .  +  A.  .  , 

i,j  i,j+i 


ABB  = 


2A  A  .  , 
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In  a  similar  manner,  using  the  new  time  velocities  the  energy  equation  is  updated  to  the  new  time 
level,  as  follows: 


e"*‘=e  .-  ®‘ 


ij  i.j  A  -L  6x 

1.J 
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where 


ABRr  n+i, 
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Por  the  cell  boundary  depths,  a  linearly  interpolated  value  is  not  used,  but,  instead,  a  combined 
geometric  and  arithmetic  average  is  used.  Por  example. 


ABR  = 


2A  A.^,  . 

1.J  1+1. j 

A. .  +  A.  .  . 

i.j  1+1, j 


This  choice  has  the  convenient  feature  that  a  cell  boundary  depth  vanishes  when  either  of  the  cell 
depths  on  each  side  of  boundary  is  zero,  i.e.,  for  this  case,  ABR  vanishes  whenever  Ajj  or  Aj+jj 
is  zero. 
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Finally,  with  the  new  time  energy  field,  the  pressure  is  updated  ns 


n+1  , 

Pij 


and  the  new  time  specific  internal  energy,  e  ,  can  be  obtained  as 


_n+l 

n+1 

e. .  =  — =4" 

»0  n+1 


Pi.j 


(11) 
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RESULTS 

Figure  E-1^  shows  the  gas  dynamic  simulation  model  used  for  the  ETVE.  The  stationary 
and  moving  valve  elements  are  approximated  as  rectangular  solid  blocks.  The  mesh  spacing  used 
is  5x  =  8y  =  1  mm.  Constant  pressure  boundary  conditions  were  applied  at  the  inlet  and  oudet 
boundaries  of  the  flow  model.  For  the  first  set  of  runs,  the  high-pressure  boundary  (upstream, 
right  boundary  in  Figure  E-1),  is  set  to  14(X)  psi  while  the  low-pressure  boundary  (downstream, 
left  boundary)  is  set  to  14  psi.  Figures  E-2  and  E-3  show  the  results  for  runs  with  combined  8% 
viscous  and  100  Ibf  coulomb  damping  on  the  Belleville  washers;  Figures  E-4  and  E-5  show  the 
results  for  runs  with  a  combined  16%  viscous  and  100  Ibf  coulomb  damping  on  the  Belleville 
washers.  In  Figures  E-2  and  E-4,  the  “a”  figure  is  a  plot  of  valve  sleeves  displacement  with  time, 
the  “b”  figure  is  a  plot  of  its  velocity  with  time,  and  the  “c”  and  “d”  figures  are  plots  of  its  accel¬ 
eration  with  time.  The  “e”  figures  show  pressure  in  the  gas  as  a  function  of  time  for  a  location  6 
mm  in  from  the  low-pressure  boundary  and  29-mm  radially  outward  from  the  axis  of  symmetry 
(labeled  as  pressure  plot  point  in  Figure  E-1).  In  Figures  E-2  and  E-4,  each  of  the  plots  also  have 
the  total  gas  static  pressure  loading  acting  on  the  valve  sleeve  as  a  function  of  time  shown  as  a 
“ghost”  dotted  line.  Figures  E-3  and  E-5  show  the  pressure  distribution  (spatial)  at  the  instant  in 
time  corresponding  to  the  last  time  step  in  each  run;  isobars  are  at  10-psi  intervals. 

Clearly,  the  gas  dynamic  force  on  the  sleeve  dominates  and  is  responsible,  when  coupled 
with  the  Belleville  spring  force,  for  the  severe  oscillations.  The  other  dynamical  forces  have  lit¬ 
tle  effect  on  the  sleeve  motion.  Examination  of  the  gas  dynamic  force  on  the  sleeve  shows  that 
this  force  predominantly  acts  left  or  toward  the  direction  of  valve  closure.  In  an  attempt  to  reduce 
this  large  gas  dynamic  force  on  the  sleeve,  a  run  was  made  in  which  the  downstream  face  of  each 
of  the  sleeve  ports  was  beveled  along  its  entire  face  at  a  45  degree  angle.  This  had  some  effect, 
but  not  enough;  the  sleeve  still  oscillated  violently. 

Based  on  the  above  observation,  the  question  was  asked,  what  might  happen  if  the  initial 
sleeve  position  was  shifted  to  the  right  with  the  sleeve  ports  on  the  other  side  of  the  stationary 
ports,  with  the  opening  motion  now  being  to  the  left  (reversed).  Then,  perhaps,  the  gas  dynami- 


a.  The  figures  are  found  as  a  section  of  the  report  beginning  on  page  E-12. 
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cal  force  on  the  sleeve  would  act  to  force  the  valve  into  the  open  position.  When  this  change  in 
geometry  was  made  in  the  simulation  model,  that  is  exactly  what  was  observed. 

In  order  to  reposition  the  valve  sleeve  as  described  above  with  the  same  sleeve  travel 
length  with  the  same  dimensions  between  ports,  it  was  necessary  to  lengthen  the  downstream  end 
of  the  sleeve  in  order  that  the  downstream  stationary  port  hole  be  totally  blocked  (closed)  in  the 
initial  configiu^tion.  Figure  E-6  shows  the  gas  dynamical  simulation  model  for  this  new  config¬ 
uration.  With  this  modification,  and  with  8%  viscous  and  100-lbf  coulomb  damping  in  the 
Belleville  springs,  the  simulation  produced  the  results  shown  in  Figures  E-7,  E-8,  and  E-9. 
Figitres  E-7a  and  E-7b  show  the  time  history  of  the  sleeve  displacement  and  velocity,  respective¬ 
ly.  Figures  E-7c  and  E-7d  show  the  valve  sleeve  acceleration  history.  Figure  E-7e  shows  the  gas 
pressure  history,  again  at  the  position  6  mm  in  from  the  left  (low)  pressure  boundary  and  29  mm 
radially  out  from  the  axis  (pressure  plot  point  in  Figure  E-6).  Figure  E-8a  and  E-8b  show  the 
phase  portrait  (displacement  versus  velocity)  for  the  sleeve  motion  with  Figure  E-8b  being  an 
enlargement  of  the  region  around  the  region  of  attraction  to  the  open  position.  Figure  E-9  shows 
the  spatial  pressure  distribution  at  the  time  corresponding  to  the  last  time  step  of  the  problem,  with 
isobars  at  10-psi  intervals.  These  results  shown  that  the  gas  dynamic  forces  do  force  the  valve 
into  a  strongly  open  position  after  3  to  4  small  bounces  from  the  Belleville  spring.  Beyond  this 
initial  opening  transient,  we  observed  that  the  valve  maintains  the  open  position  with  a  growing, 
high-frequency  oscillation  imposed  on  a  fairly  significant  steady  Belleville  spring  compression. 
A  growing,  high-frequency  oscillation  is  seen  in  the  gas  dynamic  force  on  the  sleeve.  However, 
the  average  gas  dynamic  force  on  the  sleeve  is  seen  to  be  strongly  toward  the  open  position.  It  is 
not  known  whether  this  growing  oscillation  of  the  sleeve  could  eventually  grow  to  the  point  of 
blocking  enough  of  the  flow  through  the  portholes  to  force  the  sleeve  motion  into  another  tran¬ 
sient  form.  We  do  not  know  whether  this  growing  oscillation  is  physical  or  just  numerical  arti¬ 
fice. 


For  the  second  set  of  runs,  the  high-pressure  boundary  was  increased  to  1800  psi  with  con¬ 
ditions  corresponding  to  700°F,  and  the  low-pressure  boundary  was  lowered  to  12  psi.  Figure  E- 
10  shows,  in  a  manner  similar  to  the  earlier  set  of  runs,  the  displacement,  velocity,  acceleration, 
and  single  point  pressure  histories  with  total  static  pressure  gas  load  on  the  sleeve  (dotted  line) 
included,  all  for  the  case  of  8%  viscous  damping  plus  100  Ibf  coulomb  damping.  Figure  E-11 
shows  the  corresponding  pressure  distribution  for  this  case  at  the  last  time  step,  with  12  psi  iso¬ 
bar  intervals.  Similar  data  are  shown  in  Figures  E-12  and  E-13  for  the  case  of  16%  viscous  plus 
1(X)  Ibf  coulomb  damping.  Figures  E-14  and  E-15  show  the  same  data  arrangement  displays  for 
the  modified,  backward  opening  valve.  Again,  a  high-frequency,  low-amplitude  wave  form  is 
seen  imposed  on  the  approximate  steady  state  for  later  times,  after  the  valve  is  open.  The  same 
caveat  applies.  We  are  unsure  at  this  time  of  the  physical  reality  of  this  phenomena. 

In  the  pressure  distribution  figures,  we  noted  the  large  pressure  gradients  occurring  m  the 
vicinity  of  the  low-pressure  boundary.  This  may  be  due  to  the  fixed  pressure  boundary  condition. 
If  the  flow  local  to  the  exit  boundary  is  supersonic,  it  is  no  longer  appropriate  to  impose  a  pres¬ 
sure  at  this  boundary.  To  rule  out  the  possibility  of  such  a  boundary  condition  causing  the  high- 
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frequency  oscillations  about  the  open  position,  the  code’s  boundary  condition  subroutine  was 
rewritten  to  impose  an  exit  boundary  condition  as  a  local  extrapolation  from  the  interior  of  the 
flow  field  where  the  velocities  are  locally  supersonic.  The  backward  opening  valve  with  length¬ 
ened  spool  was  executed  (8%  viscous  plus  100  Ibf  Belleville  damping)  with  the  new  boundary 
conditions.  Figme  E-16  shows  the  displacement,  velocity,  acceleration,  and  pressure  point  histo¬ 
ries,  all  plotted  with  the  total  gas  dynamic  force  history  on  the  valve  spool  (dotted  line).  Figure 
E-17  shows  the  pressure  distribution  (~13  psi  between  isobars)  at  the  last  time  step.  The  large 
pressure  gradients  that  were  piling  up  at  the  exit  boundary  are  now  gone.  However,  the  high-fre¬ 
quency  oscillation  is  still  present,  and  the  valve  opening  dynamics  are  the  same.  While  the  new 
boundary  condition  is  certainly  appropriate,  the  oscillating  phenomena  cannot  be  attributed  to  the 
previously  used,  fixed  pressure  boundary  condition. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  coupled  gas  dynamic  and  valve  slide  dynamic  simulation  model  has  been  found  to  be 
an  extremely  valuable  tool  in  assessing  the  nature  of  the  current  valve’s  operation  and  is  indis- 
pensible  in  assessing  the  effects  of  valve  modifications  or  redesigns.  Such  complex  nonlinear 
behavior  is  impossible  to  predict  in  any  way  other  than  numerical  simulation.  Even  if  the  numer¬ 
ical  simulations  are  accepted  only  on  a  qualitative  basis,  they  increase  understanding  of  behavior, 
show  trends  and  sensitivities,  and  suggest  alternate  designs.  On  the  other  hand,  the  alternative, 
cut  and  try  approach  to  design  would  be  a  long,  tedious,  and  expensive  route  to  valve 
design/redesign. 

For  future  design/analysis  work,  we  recommend  that  the  numerical  technique  used  to  sim¬ 
ulate  the  gas  dynamics  be  modified  to  give  a  partially  implicit  time  integration  instead  of  the 
explicit  time  integration.  This  will  reduce  the  CPU  time  required  to  make  runs  and  should 
increase  the  accuracy  by  reducing  the  large  number  of  time  steps  required  for  these  runs.  The 
importance  of  viscous  effects  should  also  be  examined,  especially  in  the  valve  open  position 
where  the  high-speed  flow  though  the  ports  could  possibly  place  significant  shear  forces  on  the 
sleeve,  and  therefore  influencing  its  motion. 
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Figure  E-2a.  Histories  of  valve  spool  displacement  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high- 
pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-2b.  Histories  of  valve  spool  velocity  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high-pressure 
boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E>2c.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high-pressure 
boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-2d.  Histories  of  valve  spool  acceleration  (enlargement)  and  total  gas  dynamic  force  on  valve 
spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  KXMbf  (on  Belleville)  Coulomb  damping. 
The  high-pressure  boundary  is  14(X)  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-2e.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high- 
pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-3.  Spatial  pressure  distribution  at  time  =  0.042  sec.  Isobars  are  at  ~10-psi  intervals. 
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Figure  E-4a.  Histories  of  valve  spool  displacement  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The 
high-pressure  boundary  is  1400  psia;  the  low-pressure  boimdary  is  14  psia. 


E-19 


Gas  Force  (Ibf) 


1 6  %  Viscous  +  1 00  Ibf  Bellville  Damping 


Time  (nsec) 


O.OOOOEO 


-5.0000E3 


-1.0000E4 


Figure  E-4b.  Histories  of  valve  spool  velocity  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high-pressure 
boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-4c.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  16%  viscous  (on  spool)  plus  KXllbf  (on  Belleville)  Coulomb  damping.  The  high-pressure 
boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-4d.  Histories  of  valve  spool  acceleration  (enlargement)  and  total  gas  dynamic  force  on  valve 
spool.  Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping. 
The  high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-4e.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The 
high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  E-5.  Spatial  pressure  distribution  at  time  =  0.0443  sec.  Isobars  are  at  ~10-psi  intervals. 
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Figure  E-6.  Proposed  backward-opening  ETVE  gas  dynamic  simulation  model  All  dimensions  are  mm. 
(The  schematic  is  not  to  scale.) 


Displacement  (in) 


8  %  Viscous  +  1 00  Ibf  Bellville  Damping 
Backward  Opening  /  Lengthened  Spool 


Time  (ixsec) 


Figure  E-7a.  Histories  of  modified  valve  spool  displacement  and  total  gas  dynamic  force  on  modified 
valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb 
Damping.  The  high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  7-b-  Histories  of  modified  valve  spool  velocity  and  total  gas  dynamic  force  on  modified  valve 
spool  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb  Damping. 
The  high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  7-C.  Histories  of  modified  valve  spool  acceleration  and  total  gas  dynamic  force  on  modified 
valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  BeUeviUe)  Coulomb 
Damping.  The  high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  7-d.  Histories  of  modified  valve  spool  acceleration  (enlarged)  and  total  gas  dynamic  force  on 
modified  valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville) 
Coulomb  Damping.  The  high-pressure  boundary  is  1400  psia;  the  low-pressure  boundary  is  14  psia. 
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Figure  7-e.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  modified  valve 
spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb  Damping. 
The  high-pressure  boundary  is  1400  psia;  the  low-pressure  botmdary  is  14  psia. 
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Figure  8-a.  Phase  plane  portrait  for  modified  valve  spool  motion.  Damping  forces  are  set  at  8%  viscous 
(on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb  damping.  The  high-pressure  boundary  is  1400  psia;  the 
low-pressure  boundary  is  14  psia. 
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Figure  E-8b.  Phase  plane  portrait  enlargement  of  node  region  (Figure  E-8a)  around  the  valve  open  positioa 
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Figure  E-9.  Spatial  pressure  distribution  about  modified  valve  spool  at  time  =  0.090  sec.  Isobars  are  at 
~10-psi  intervals. 
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Rgure  E-10a.  Histories  of  valve  spool  displacement  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damjring.  The  high- 
pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Rgure  E-10b.  Histories  of  valve  spool  velocity  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high-pressure 
boimdary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-10c.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  valve  spooL 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high- 
pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-1  Od.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The  high- 
pressure  boundary  is  18(X)  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-11.  Spatial  pressure  distribution  at  time  =  0.0432  sec.  Isobars  are  at  12-psi  intervals. 
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Figure  E-12a.  Histories  of  valve  spool  displacement  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The 
high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 


E-39 


Gas  Force  (Ibf) 


1 6%  Viscous  + 1 00  Ibf  Bellville  Damping 


Time  (usee) 


O.OOOOEO 


-5.0000E3 


-1.0000E4 


-1.5000E4 


Rgure  E'12b.  Histories  of  valve  spool  velocity  and  total  gas  dynamic  force  on  valve  spool.  Damping 
forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  BeUevUle)  Coulomb  damping.  The  high-pressure 
boundary  is  18(X)  psia;  the  low-pressure  boundary  is  12  psia. 
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Rgure  E-12c.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  valve  spool 
Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The 
high-pressure  boundary  is  1800  psia;  the  low-pressure  bound^  is  12  psia. 
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Figure  E-12ci.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  valve  spool. 
Damping  forces  are  set  at  16%  viscous  (on  spool)  plus  lOOlbf  (on  Belleville)  Coulomb  damping.  The 
high-pressure  bormdary  is  18(X)  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-13.  Spatial  pressure  distribution  at  time  =  0.0432  sec.  Isobars  are  at  12-psi  intervals. 
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Figure  E-14a.  Histories  of  modified  valve  spool  displacement  and  total  gas  dynamic  force  on  modified 
valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  1(X)  Ibf  (on  Belleville)  Coulomb 
Damping.  The  high-pressure  boundary  is  18(X)  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-14b.  Histories  of  modified  valve  spool  displacement  (enlargement)  and  total  gas  dynamic  force 
on  modified  valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville) 
Coulomb  Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-14c.  Histories  of  modified  valve  spool  velocity  and  total  gas  dynamic  force  on  modified  valve 
spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb  Damping. 
The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boimdary  is  12  psia. 
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Figure  E-14d.  Histories  of  valve  spool  velocity  (enlargement)  and  total  gas  dynamic  force  on  modified 
valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb 
Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Rgure  E-14e.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  modified  valve  spool. 
Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb  Damping.  The 
high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-14f.  Histories  of  valve  spool  Acceleration  (enlargement)  and  total  gas  dynamic  force  on  mod¬ 
ified  valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb 
Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary  is  12  psia. 
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Figure  E-14g.  Pressure  history  at  the  pressure  histoiy  point  and  total  gas  dynamic  force  on  modified 
valve  spool.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on  Belleville)  Coulomb 
Damping.  The  high-pressure  boundaiy  is  1800  psia;  the  low-pressure  boimdary  is  12  psia. 
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Rgure  E“16a.  Histories  of  modified  valve  spool  displacement  and  total  gas  dynamic  force  on  modified 
valve  spool,  with  modified  boundary  condition.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100 
Ibf  (on  Belleville)  Coulomb  Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boimd- 
ary  is  12  psia. 
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Figure  E-1 6b.  Histories  of  modified  valve  spool  velocity  and  total  gas  dynamic  force  on  modified  valve 
spool,  with  modified  boundary  condition.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  ICX)  Ibf 
(on  Belleville)  Coulomb  Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  boundary 
is  12  psia. 
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Rgure  E-1 6c.  Histories  of  valve  spool  acceleration  and  total  gas  dynamic  force  on  modified  valve  spool, 
with  modified  boundary  condition.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100  Ibf  (on 
Belleville)  Coulomb  Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  botmdary  is  12 
psia. 


E-54 


Gas  Force  (Ibf) 


8%  Viscous  +  100  Ibf  Bellville  Damping 
Backward  Opening  /  Lengthened  Spool  /  New  B.C. 


Time  (usee) 


Figure  E-16cl.  Pressure  history  at  the  pressure  history  point  and  total  gas  dynamic  force  on  modified 
valve  spool,  with  modified  boundary  condition.  Damping  forces  are  set  at  8%  viscous  (on  spool)  plus  100 
Ibf  (on  Belleville)  Coulomb  Damping.  The  high-pressure  boundary  is  1800  psia;  the  low-pressure  bound¬ 
ary  is  12  psia. 
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Figure  E-17.  Spatial  pressure  distribution  about  modified  valve  spool  at  time  =  0.087  sec.  Isobars  are  at 
13-psi  intervals. 
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1.  SUMMARY  OF  RESULTS 


Revision  1  to  this  report  employed  an  updated  acceleration  history  for  the  backward 
opening  piston  and  sleeve.  Section  7  contains  that  data. 

Selected  components  of  the  Eaton  Throat— Valve  Element  were  evaluated  with  respect  to 
the  rules  of  the  ASME  Boiler  and  Pressure  Vessel  Code,  Section  VIII  (Rules  for 
Construction  of  Pressure  Vessels).  An  exception  was  taken  forthe  valve  element  because 
one  of  its  materials  (AMS  5643  in  plate  form)  was  not  approved  for  use  by  the  ASME  Code 
Section  VIII. 

The  body  of  the  valve  contained  three  rows  of  rectangular  portholes  sized  0.76”  x  0.50” 
each.  The  analysis  showed  that  the  valve  body  could  withstand  the  design  conditions  of 
2,000  psi  external  pressure  at  700  °F  with  its  current  porthole  size  and  configuration.  It  also 
showed  that  the  portholes  could  be  enlarged  to  0.76”  x  0.76”  square  with  1/8”  radius 
corners  while  under  design  conditions  and  remain  stable,  maintaining  a  factor  of  safety  of 
2  against  buckling. 

However,  the  dynamic  forces  associated  with  a  backward  opening  piston  and  sleeve  (used 
to  ’’open”  the  portholes)  caused  excessive  stresses  throughout  the  valve.  The  location  and 
magnitude  of  selected  elastically  calculated  stresses  were:  (1)  the  valve  body  around  the 
portholes  (using  the  current  porthole  size)  exceeded  58  ksi,  (2)  the  piston  stop  plate 
exceeded  400  ksi,  and  (3)  the  stop  plate  attachment  bolts  reached  224.9  ksi.  All  above 
elastically  calculated  stresses  exceeded  the  allowable  stress  level  of  32.0  ksi  defined  by 
the  ASME  Boiler  and  Pressure  Vessel  Code,  Section  VIII  (Rules  for  Construction  of 
Pressure  Vessels). 

The  large  dynamic  (impulse)  forces  created  by  the  deceleration  of  the  piston  and  sleeve 
were  primarily  due  to  the  stiffness  of  the  Belleville  washer  stack.  The  stackwas  clearly  very 
’’stiff”,  and  had  a  short  travel  before  bottoming  out.  The  piston  impulse  forces  would  be 
greatly  reduced  if  either  of  the  following  modifications  were  made: 

(1 )  Arrange  (or  add  to)  the  stack  of  Belleville  washers  so  that  the  effective  travel  length 
(before  bottoming  out)  is  much  longer.  When  the  washers  bottomed  out  their  stiffness 
increased  tremendously,  thus  increasing  the  piston/sleeve  decelerations. 

(2)  Employ  Belleville  washers  (or  other  springs)  that  have  a  lower  stiffness  value  and  a 
longer  travel  length.  This  will  reduce  the  magnitude  of  the  piston/sleeve  decelerations. 
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2.  INTRODUCTION 

2.1.  General:  The  Eaton  Throat-Valve  Element  Prototype''  is  a  quick-opening  valve 
designed  to  vent  gas  at  2,000  psig.  The  valve  element  consisted  of  a  cylindrical  body 
with  one  end  blocked  and  the  other  open.  The  body  had  three  rows  of  slotted  portholes 
around  the  circumference.  An  inner  cylindrical  sleeve  also  had  three  rows  of  portholes 
through  which  pressurized  passed  when  the  sleeve  was  positioned  so  that  its  portholes 
coincided  with  those  of  the  body.  The  sleeve  position  was  controlled  by  a  piston,  which 
was  in  turn  located  by  a  control  system.  Two  components  of  the  valve  were  evaluated 
in  this  analysis:  (1 )  the  valve  body  at  the  location  of  the  portholes,  and  (2)  the  end  plate 
that  stops  the  travel  of  the  piston  when  the  sleeve  is  in  the  ’’open  porthole”  position. 

2.2.  Design  Code:  The  design  code  followed  was  the  ASME  Boiler  and  Pressure 
Vessel  Code  Section  VIII,  Rules  for  Construction  of  Pressure  Vessels^  (ASME  Code). 
Strict  observance  of  the  ASME  Code  requirements  was  not  always  kept  -  specific 
exceptions  will  be  noted  in  this  report  as  they  were  employed. 

2.3.  Design  Conditions:  The  design  conditions  on  the  valve  were:  2,000  psi  external 
pressure,  and  700  °F  material  temperature. 

2.4.  Materials:  All  components  that  were  evaluated  were  constructed  of  the 
precipitation-hardened  stainless  steel  AMS  5643  (SAE  Standard)'^.  This  material  had 
an  equivalent  in  the  ASTM  Standards^,  A-564  for  bars,  A-693  for  plates,  and  A-705  for 
forgings.  The  ASME  Code  gave  the  properties  for  this  material  as  listed  in  Table  1 . 


Table  1.  Material  Properties  for  ASTM  A-564,  A-693,  &  A-705 


Temperature  °F 


Property 

100 

200 

300 

400 

500 

600 

650 

700 

Yield  Strength,  ksi 

115.0 

106.3 

101.9 

98.3 

95.2 

92.8 

91.5 

*[90.0] 

Modulus  of  Elas¬ 
ticity  X  10“®,  psi 

28.3 

27.6 

27.0 

26.5 

25.8 

25.3 

* 

[25.1] 

24.8 

ASME  Code  Al¬ 
lowable  Stress,  ksi 

35.0 

35.0 

35.0 

34.1 

33.3 

32.8 

32.6 

★ 

[32.0] 

‘Properties  not  specified  at  this  temperature.  Number  in  brackets  was  either  interpolated  or  extrapolated. 


The  minimum  yield  strength  specified  by  the  ASME  Code  (shown  in  Table  1)  was 
somewhat  lower  than  that  specified  by  the  AMS  5643  Standard.  This  appeared  to  be 
because  the  ASME  Code  applied  a  reduction  factor  (factor  of  safety)  to  the  yield  strength 
before  it  calculated  allowable  stress  levels  (from  that  yield  strength).  It  was  noted  that  the 
ASME  Code  only  approved  the  bar  material  (A564)  for  Section  VIII  Code  use.  This  design 
used  the  approved  bar  as  well  as  unapproved  plate  material.  An  exception  was  taken  from 
the  ASME  Code  for  the  plate  material. 
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3.  VALVE  BODY  PORTHOLE  EVALUATION 


3.1.  Purpose:  The  purpose  of  the  valve  body  porthole  evaluation  was  to  determine  if 
and  (if  so)  how  much  the  body  portholes  could  be  widened,  given  the  design  conditions. 
This  was  done  by:  (1)  checking  the  body  to  determine  the  minimum  required  thickness 
at  design  conditions,  (2)  checking  the  portholes  to  determine  the  wall  thickness  around 
them  required  for  reinforcement,  and  (3)  determining  the  maximum  allowable  size  of 
the  portholes  based  on  the  actual  valve  body  wall  thickness.  One  additional  condition 
that  was  evaluated  on  the  portholes  was  due  to  the  force  produced  by  the  piston, 
operating  in  a  backward  opening  mode,  on  the  stop  plate.  That  force  would  be  carried 
from  the  stop  plate  through  its  attachment  bolts  to  the  valve  body. 

3.2.  Valve  Body  Minimum  Thickness  Under  Design  Conditions:  The  ASME  Code 
developed  charts  to  indicate  the  minimum  required  thickness  for  cylindrical  components 
subject  to  external  pressure.  However,  since  the  plate  material  that  the  valve  body  was 
made  of  was  not  approved  by  the  Code,  no  external  pressure  chart  was  available. 

(The  ASME  Code  did  not  set  up  a  chart  for  the  approved  bar  material  because  it  was 
not  expected  to  be  used  in  vessel  shell  constructions).  Therefore,  the  following  process 
was  followed  to  establish  a  minimum  required  thickness  for  the  cylindrical  wall  of  the 
valve  body. 

S.  P.  Timoshenko®  developed  a  method  for  determining  the  thickness  required  for  a 
cylindrical  vessel  subject  to  external  radial  and  end  pressure.  It  employed  the  modulus  of 
elasticity  and  Poisson’s  Ratio  of  the  material,  the  cylinder  dimensions,  and  a  variable 
buckling  node  term.  The  buckling  node  term  accounted  for  the  fact  that  the  cylinder  could 
collapse  flat  (2  nodes),  collapse  in  a  three  corner  shape  (3  node),  and  so  on.  An  initial 
assumption  for  his  solution  was  that  the  cylinder  edges  were  simply  supported. 
Weingarten®  further  simplified  Timoshenko’s  equation  for  cylinders  that  met  the  following 
additional  requirement: 

•  60  <  (l/r)2  *  (r/t)  <  2.5  *  (r/t)2 

The  parameters  for  this  cylinder  were:  I  =  unsupported  length  =  10.3”,  r  =  outside  radius 
=  5.872  /  2  =  2.936”,  and  t  =  thickness  of  wall  (t  nominal  =  0.339”).  This  gave: 

•  60  <  (1 0.3/2.936)2  *  (2.936/0.339)  =  107  <  2.5  *  (2.936/0.339)2  =  188 
Therefore,  the  simplified  equation  was  used,  and  was  as  follows: 

•  q’  =  (0.92  *  E)  /  ({l/r}  *  {r/t}2-5),  where  E  =  the  modulus  of  Elasticity  at  700  °F 

The  term  q’  was  the  critical  buckling  pressure  for  the  cylinder.  It  is  important  to  note  that 
actual  tests  have  shown  that  a  cylinder  may  buckle  at  pressures  that  are  20%  higher  or 
lowerthat  the  calculated  critical  buckling  pressure^.  Table  2  shows  the  required  valve  body 
thickness  for  various  factors  of  safety  against  the  2,000  psi  external  pressure. 
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Table  2.  Minimum  Valve  Body  Wall  Thickness  Required  vs.  Factors  of  Safety 
(External  Pressure  of  2,000  psi) 

Factor  of  Safety  Against  Buckling 

Thickness  Required,  in. 

1.0 

0.115 

1.25 

0.126 

1.43 

0.133 

1.67 

0.141 

2.0 

0.152 

-  10.0 

0.290 

Since  the  actual  valve  body  wall  was  0.339”  thick  Table  2  shows  that  a  factor  of  safety  of 
10  against  buckling  was  maintained  under  a  2,000  psi  external  pressure  load. 


:  Valve  Body  Thickness  was  Adequate! 

3.3.  Required  Reinforcement  for  Portholes  Under  Design  Conditions:  The  ASME 
Section  VIII  Code  defined  the  required  reinforcement  for  holes  in  the  cylindrical  shell  of 
the  valve  body.  Paragraph  UG-37  (d)  (1 )  stated  that  the  required  reinforcement  around 
each  porthole  was  50%  of  that  required  by  paragraph  (c)  which  was: 

•  Reinforcement  Area  Required  =  d  *  tr  *  F 

■  where:  d  =  cord  length  of  porthole  opening  =  0.760” 

tr  =  required  thickness  of  Table  2 
F  =  correction  factor  =1.0  for  shells 

The  dimensioned  portholes  are  shown  in  Figure  1  in  solid  lines.  The  valve  body  has  three 
circumferential  rows  of  1 6  portholes,  each  circumferentially  spaced  as  shown  on  the  figure. 
The  smaller  circles  indicate  an  equivalent  hole  with  a  diameter  equal  to  the  cord  length 
of  the  portholes.  The  areas  between  the  larger  and  smaller  circles  are  those  where 
reinforcement  must  be  in  order  to  qualify  as  material  available  to  reinforce  the  portholes. 
Since  the  portholes  were  (area  wise)  smaller  than  the  smaller  circles,  the  calculations  for 
reinforcement  available  will  be  somewhat  conservative.  However,  the  ASME  Code  does 
not  define  that  magnitude  of  conservatism. 

The  areas  required  for  reinforcement  of  the  portholes  overlap  as  shown  on  Figure  1  in  the 
circumferential  direction  only  (each  row  is  longitudinally  spaced,  from  center  to  center,  at 
2.214”  apart,  which  equated  to  about  3  times  the  inner  circle  diameter).  On  the  outside 
diameter  of  the  valve  body  the  areas  do  not  cross  into  the  portholes,  but  they  do  on  the 
inside  diameter.  However,  it  was  assumed  that  the  material  in  the  reinforcement  areas  lost 
to  the  portholes  was  made  up  by  the  extra  material  available  in  the  longitudinal  direction. 
The  metal  available  for  reinforcing  the  portholes  was  calculated  using  ASME  Code  Figure 
UG-37. 1 ,  as  follows: 

•  Metal  Available  =  d  *  (Ei*t  -F*tr)  -  2*tn  *  (Ei*t  -  F*tr)  *  (1  -  fri) 

where:  fri  =1.0  (no  nozzle  reduction) 

E-|  =  1 .0  (no  welds  at  porthole) 

F  =  1 .0  correction  factor 
which  reduced  Metal  Available  =  d  *  (t  -  tr) 
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Date;  December  1993 


Porthole  Dimensions  on  Outside  Diameter  of  Valve  Body 


Circumferential  Direction 


Porthole  Dimensions  on  Inside  Diameter  of  Valve  Body 


Figure  1.  Porthole  Dimensions  on  Valve  Body  Outside  &  Inside  Diameters 
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Equating  the  required  area  (reduced  by  50%  as  allowed)  and  the  available  area  (dividing 
the  available  area  of  metal  between  two  adjacent  portholes  equally)  gave; 

•  1/2*d*tr  <  d*  1/2*(t-tr) 

-or- 

2tr<  t 

The  thickness  of  the  valve  body  around  the  portholes  had  to  be  twice  that  required  to 
maintain  stability  of  the  wall.  Table  3  compares  the  required  thickness  for  stability  listed 
in  Table  2,  doubles  it  for  that  required  for  reinforcement,  and  lists  the  actual  wall  thickness. 


Table  3.  Valve  Body  Wall  Thickness  Required  at  Portholes 
(2,000  psi  External  Pressure) 


Factor  of  Safe¬ 
ty  Against 
Buckling 

Thickness  Required 
for  Stability,  in. 

Thickness  Required 
for  Reinforcement,  in. 

Actual  Valve  Body 

Wall  Thickness,  in. 

1.0 

0.115 

0.230 

0.339 

1.25 

0.126 

0.252 

0.339 

1.43 

0.133 

0.266 

0.339 

1.67 

0.141 

0.282 

0.339 

2.0 

0.152 

0.304 

0.339 

10.0 

0.290 

0.580 

0.339 

:  The  reinforcement  provided  by  the  valve  body  wall  at  the  portholes  allowed 
a  2,000  psi  external  pressure  with  a  factor  of  safety  of  2  against  buckling! 

3.4.  Maximum  Allowable  Porthole  Size  Under  Design  Conditions:  The  evaluation 
reported  in  the  previous  subsection  used  a  ”d”  value  equal  to  the  cord  length  of  the 
porthole  opening.  That  ”d”  value  would  be  the  same  if  the  hole  were  circular,  with  a 
diameter  equal  to  ”d”.  If  the  hole  were  square  (with  a  corner  radius)  at  the  current 
dimension  of  0.760”  as  shown  on  Figure  1 ,  the  cord  length  would  still  be  that  which  was 
used  above.  A  porthole  with  larger  dimensions  (than  0.760”  square)  would  not  maintain 
a  factor  of  safety  of  2  against  the  theoretical  buckling  pressure  on  the  valve  body. 

:  Maximum  recommended  porthole  size  was  0.760”  square,  with  a  corner 
radius  value  of  0.1 25”!  This  maintained  a  factor  of  safety  of  2  against  buckling  under 
a  2,000  psi  external  pressure. 

3.5.  Valve  Body  Under  Piston  Loads:  Section  4.4  calculates  a  maximum  equivalent 
static  force  of  103,518  lbs  that  will  be  applied  to  the  valve  body  via  the  stop  plate.  That 
force  was  produced  by  the  dynamics  of  the  sleeve  being  positioned  by  the  piston 
operating  in  a  backward  opening  mode.  The  force  is  compressive  on  the  body  of  the 
valve,  and  will  produce  higher  stresses  on  the  row  of  "webs”.  Each  web  is  that  material 
between  portholes  in  the  circumferential  direction.  The  total  area  of  the  webs,  and  the 
resulting  stresses  in  them  would  be: 


8 


F-9 


siihjfir:tr  Analysis  of  Selected  Eaton  Throat-Valve  Element  Components  Date:  December  1993 

Prepared  by:  S.  D.  Snow _ Checked  by  G.  K.  Miller _ 

Revision  1 


•  Web  Area  =  1 6  Webs  *  Area  per  Web  (^=«{(0.26  +  0.0665)  *  0.339}) 

=  16*0.111 
=  1.77  ln.2 

•  Compressive  Stress  in  Webs  =  Force  /  Area 

=  lbs/ 1.77  in.2 

=  58,485  psi 

From  Table  1  the  allowable  stress  level  of  the  material  that  the  valve  body  was  constructed 
of  was  32.0  ksi.  The  actual  compressive  stress  caused  by  the  maximum  force  of  the 
backward  opening  piston  {58,485  psi)  exceeded  the  allowabie  stress  level.  It  was  also 
suspected  that  the  valve  body  wall  would  be  unstable  under  the  above  piston  load. 
However,  no. further  calculations  were  performed  since  the  ailowable  compressive  stress 
was  unacceptable. 

:  Valve  Body,  in  its  Current  Porthole  Configuration  Could  Not  Carry  the  Piston 
Impact  Loading  from  Backward  Opening! 
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4.  PISTON  STOP  PLATE 


4.1.  General:  The  end  plate  that  stops  the  travel  of  the  sleeve-positioning  piston  was 
made  of  AMS  5643  steel.  The  stop  plate  consisted  of  a  3.75”  diameter  disc,  1 .0”  thick. 

A  central  hole  allowed  the  piston  guide  rod  to  be  attached  to  the  valve  sleeve.  Two 
180°  opposed  holes  were  provided  for  piston  actuating  pressure  ports,  along  with  six 
0.33”  holes  for  mounting  the  stop  to  the  valve  body  with  5/1 6”  bolts. 

4.2.  Load  Conditions:  The  loading  on  the  stop  plate  contained  several  components: 

(1)  the  actuation  pressure  on  the  piston,  (2)  the  force  caused  by  the  nitrogen  gas 
flowing  through  the  portholes,  and  (3)  the  dynamics  of  the  piston  during  its  actuation. 

4.3.  Stop  Plate  Dynamics:  A  finite  element  model  of  the  stop  plate  was  developed  in 
order  to  determine  its  fundamental  frequency.  This  was  needed  to  see  how  that 
frequency  compared  to  the  frequencies  of  the  input  forces.  The  stop  plate  was 
modeled  and  analyzed  using  the  l-DEAS®  software.  Parabolic  tetrahedron  elements 
(1 ,355  total)  simulated  the  stop  itself  and  spring  elements  took  the  place  of  the 
attachment  bolts.  Two  boundary  conditions  were  considered  possible  for  the  stop 
plate:  (1)  the  outside  edge  of  the  stop  plate  was  simply  supported  (this  assumed  that 
the  bolts  always  kept  the  stop  plate  in  contact  with  the  valve  body),  and  (2)  the  stop 
plate  was  unrestrained  except  by  the  springs  that  took  the  place  of  the  attachment  bolts 
(this  assumed  that  some  small  gap  existed  between  the  stop  plate  and  the  valve  body). 
Figures  2  &  3  show  the  top  and  isometric  views  of  the  model  (only  half  of  the  stop 
needed  due  to  symmetrical  geometry).  Material  properties  from  at  700  °F  used  were: 
Modulus  of  elasticity  of  24.8x10®  psi,  Poisson’s  ratio  of  0.29. 

Figure  4  shows  the  fundamental  frequency  of  the  stop  plate  with  simply  supported  edges. 
This  mode  shows  the  stop  plate  center  flexing  in  and  out  at  12,367  Hz.  The  valve  body 
will  allow  the  outward  movement  but  will  restrain  the  inward  movement  of  the  stop  plate. 
The  second  natural  frequency,  shown  in  Figure  5,  shows  the  stop  plate  shifting  from 
side— to— side.  The  geometry  of  the  valve  body  should  impede  this  mode. 

Figures  6  &  7  show  the  first  two  natural  frequencies  of  the  stop  plate  that  was  restrained 
only  by  the  attachment  bolts  (represented  by  nodal  rigid  restraints).  Their  frequencies 
were  at  7,1 36  and  9,256  Hz.  Much  of  their  motion  will  also  be  restrained  by  the  valve  body. 

One  other  possible  motion  of  the  stop  plate  not  evaluated  in  the  finite  element  model  was 
that  of  the  mass  of  the  stop  plate  riding  on  six  attachment  bolts  (i.e.,  rigid  motion  of  the  stop 
plate  on  springs).  The  spring  constant  of  each  bolt  was: 

•  Spring  Constant  K  =  Modulus  of  Elasticity  *  Bolt  Area  /  Spring  Length  of  bolt 

=  (24.8x1 0®)  *  0.077  in.2  / 1 .0  in. 

=  1 ,909,600  Ibs/in.  per  bolt 

•  Total  Spring  Constant  Kt  =  6  bolts  *  K  per  bolt  =  11,457,600  Ibs/in. 

The  weight  of  the  stop  plate  was  about  3  lbs.  The  fundamental  frequency  of  this 
mass-spring  system  was 

•  fi  =  (1  /  2jc)  *  (Kt  /  mass  of  stop  plate) ‘'^2 

=  (1  /  2n)  *  (1 1 ,457,600  /  {3  /  386.4})l/2 

=  6,114  Hz 
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Figure  2.  Model  of  Stop  Plate  (Top  View) 


Figure  3.  Model  of  Stop  Plate  (Isometric  View) 
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Figure  4.  Mode  Shape  #1  for  Simply  Supported  Edges  at  12,367  Hz 


Figure  5.  Mode  Shape  #2  for  Simply  Supported  Edges  at  12,895  Hz 
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Figure  6.  Mode  Shape  #1  for  Stop  Plate  with  Free  Edges  at  7,136  Hz 
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Therefore,  of  the  three  conditions  on  the  stop  plate  the  minimum  calculated  natural 
frequency  was  6,114  Hz. 

4.4.  Forces  on  Stop  Plate:  The  actual  forces  applied  to  the  stop  plate  were: 

•  Piston  Actuation  Force:  The  actuation  force  on  the  piston  started  at  0, 
linearly  increased  to  1 ,425  lbs  at  0.05  seconds,  continued  to  increase  at  a 
somewhat  slower  rate  to  2,850  lbs  at  0.5  seconds  where  it  remained  constant 
(See  Reference  9). 

•  Gas  Forces  on  Sleeve:  The  gas  forces  on  the  portholes  in  the  sleeve 
began  at  zero  at  time  0  and  remained  there  until  the  piston  and  sleeve  slid 
just  enough  to  allow  flow  through  the  portholes.  At  that  time  (roughly  0.01 7 
seconds)  the  force  fluctuated  until  the  piston  stopped  moving  (bouncing  on 
the  washers  and  the  stop,  at  about  0.027  seconds).  Then  the  gas  force 
oscillated  in  a  pattern  that  increased  in  magnitude  with  time  (See  Ref.  1 0). 

•  Piston  impact  Loads:  The  piston  actuation  force  and  the  gas  force  on  the 
sleeve  caused  the  motion  of  the  piston.  When  the  piston  contacted  the  stop 
plate  a  spring-type  washer  stack  (Belleville  washers)  compressed.  The 
dynamics  of  the  piston  from  the  initial  contact  to  when  it  came  to  rest  on  the 
stop  plate  were  given  in  Reference  1 0.  The  force  history  associated  with  the 
piston  impacting  the  stop  was  determined  by  multiplying  the  acceleration 
history  of  the  piston  by  its  mass  (inc.  the  mass  of  the  attached  sleeve  and 
etc.).  Those  forces  were  small  until  the  piston  moved  the  sleeve  to  ’’just 
begin  to  open”  the  portholes,  at  which  time  the  force  increased  until  the 
piston  impacted  the  washers  and  the  stop.  Four  piston  impact  loads 
occurred  until  the  piston  settled  onto  the  stop,  as  follows: 


Table  4.  Piston  Impulse  Loads 


Piston  Impulse 
or  Impact 

Max.  Accelera¬ 
tion,  in./sec2 

Associated 
Force,  lbs 

Frequency  of 
Impulse,  Hz 

Freq.  between 
this  and  next 
Impulse,  Hz 

First 

3.0E6 

97,200 

1,515 

382 

Second 

1.85E6 

59,940 

1,316 

649 

Third 

1.0E6 

32,400 

1,000 

794 

Fourth 

4.3E5 

13,932 

1,667 

(last  impulse) 

Note  that  the  piston  impact  frequencies  were  all  less  that  the  lowest  natural 
frequency  (6,11 4  Hz)  of  the  stop  plate.  The  magnitude  of  the  piston  impact  force 
that  was  actually  transmitted  into  the  stop  plate  as  an  equivalent  static  force  was 
dependent  on  the  piston  impact  frequency  and  the  natural  frequency  of  the  stop 
plate.  The  transmission  ratio  for  the  first  and  most  severe  impulse  was 
calculated  as  follows: 

Transmission  Ratio’’"'  =  1  /  {1  -  (Forcing  Freq.^  /  Natural  Freq.^)} 

=  1  /  {I  -(7,5152/6,1142)} 

=  1.065 
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This  meant  that  during  the  application  of  the  first  impulse  load  the  equivalent 
static  force  applied  to  the  stop  plate  was  97,200  lbs  *  1.065  =  103,518  lbs\  All 
other  impulse  loads  were  much  smaller. 

Note  that  if  a  triangular  impulse  approach^^  was  taken  to  determine  the 
’’Dynamic  Load  Factor”  (equivalent  to  transmission  ratio)  the  results  would  give 
about  the  same  value  (1 .06). 

•  Total  Load  on  Stop  Plate:  The  two  cases  of  interest  were: 

(1)  Maximum  equivalent  static  load  on  stop  plate  before  the  piston  came  to 
rest  (from  maximum  piston  impulse  force)  =  103,518  lbs 

..  (2)  Maximum  load  on  stop  plate  after  impulse  loads  ceased 
=  actuation  force  +  gas  force 
=  2,850  lbs  +  Increasing  Gas  Forces 
Gas  forces  and  their  associated  frequencies  (<  1 ,000  Hz)  did  cause 
amplification  in  the  stop  plate  response  since  the  stop  plate’s  lowest  natural 
frequency  was  so  much  higher  (6,114  Hz).  This  load  could  not  be  evaluated 
since  it  was  not  known  if  the  gas  forces  continued  to  increase,  or  stabilized 
with  time.  However,  it  was  not  expected  that  the  gas  loads  would  approach 
those  of  the  piston  impulses. 

Only  the  maximum  equivalent  piston  impact  static  load  of  103,518  lbs  was  used  for 
calculating  stresses. 

4.5.  Stresses  in  Stop  Plate  Bolts:  The  stresses  in  the  stop  plate  bolts  were  tensile, 
due  to  the  103,518  lb  maximum  equivalent  static  load  from  the  piston  maximum 
impulse  force.  Six  5/16”  attachment  bolts  were  machined  from  the  same  material  as 
the  valve  body  and  stop  plate,  AMS  5643.  The  elastically  calculated  stresses  in  the 
bolts  were: 

•  Tensile  Stress  =  Force  /  Area 

=  103,518  lbs  /  {6  *  (jt/4)*(5/16/2} 

=  224.9  ksi 

The  elastically  calculated  stress  in  the  bolts  {224.9  ksi)  exceeded  the  tensile  strength  of 
the  material(190  ksi  per  Ref.  3). 

:  Bolts  Could  Not  Carry  the  Piston  Impact  Loading! 


What  bolt  size  and  quantity  would  be  needed  to  adequately  carry  the  loads?  Table  1 
showed  that  the  allowable  tensile  stress  in  the  specified  material  was  32.0  ksi.  The 
required  area  of  the  bolts  to  give  a  tensile  stress  at  32.0  ksi  was: 

•  Bolt  Area  Required  =  Force  /  Allowable  Stress 

=  103,519  lbs  /  32,000  psi 

=  3.23  ln.2 

Table  5  shows  sizes  and  quantities  of  bolts  required  to  limit  the  tensile  stress  to  32.0  ksi. 
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Table  5.  Bolts  Required  to  Carry  Maximum  Piston  Impulse  Load 


Bolt  Size,  in. 

Tensile  Area  per  Bolt,  in.2 

Number  of  Bolts  Required 

1 

0.606 

6 

7/8 

0.462 

7 

3/4 

0.334 

10 

5/8 

0.226 

15 

9/16 

0.182 

18 

4.6.  Stresses  in  the  Stop  Plate:  The  finite  element  model  of  the  stop  plate  was 
employed  in  a  //near  elastic  analysis.  The  analysis  used  the  equivalent  static  load  of 
92,728  lbs  (load  used  in  the  original  release  of  this  report).  That  force  was  converted 
into  a  pressure  load  on  the  face  that  contacts  the  piston.  The  piston  (specifically  the 
spring-type  Belleville  washers)  did  not  contact  the  entire  raised  back  face  of  the  stop 
plate,  shown  cross-hatched  in  Figure  8  below,  but  that  was  assumed  for  a  first 
approximation. 

•  Pressure  Load  =  Force  /  Area  of  Application 

=  92,728  lbs  /  {{jt/4)*  (1 .426^  -  0.7802)} 

=  82,850  psi 

Revision  1  Note:  Since  a  linear  elastic  analysis  was  used  in  the  original  analysis,  those 
results  were  scaled  to  determine  the  results  for  the  higher  load  of  103,518  lbs. 

Scale  Factor  =  103,518/92,728  =  1.116 


Back  Raised  Face  of 
Stop  Plate  (do  = 

1 .426”,  di  =  0.780”) 


Figure  8.  Assumed  Piston  Contact  on  Stop  Plate 
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The  above  equivalent  pressure  load  of  (82,850  *  1.116  =)  92,461  psLwas  also  the  bearing 
stress  at  the  face  where  the  force  was  applied.  That  stress  exceeded  the  level  allowed  by 
the  ASME  Code,  listed  in  Table  1  of  32,000  psi.  If  the  load  was  distributed  over  the 
somewhat  smaller  area  that  the  piston  washers  actually  contact  that  value  of  stress  would 
be  even  greater. 

The  finite  element  model  with  restraints  at  the  attachment  bolts  only  was  stressed  as 
shown  in  Figure  9.  The  maximum  elastically  calculated  von  Mises  stress  level  was  located 
under  the  bolts  in  the  stop  plate  at  a  magnitude  of  (362.0  *  1.116  =)  404.0ks\. 

:  Stop  Plate  Could  Not  Carry  the  Maximum  Piston  Impact  Force! 


Eaton  Throat-Valve  Stop 

LOAD  SET:  1  -  PISTON  LOAD 
FRAME  OF  REF:  GLOBAL 

STRESS  -  VON  MISES  MIN:  2199.83  MAX:  3 . 62E+05 


Max.  Stress  of 
362.0  ksi 


Stress  fpsi) 

3.1 1E5 

2.59E5 

2.08E5 

1.56E5 

1.05E5 

53,608 


Figure  9.  Stress  Contour  in  Stop  Plate  Under  92,728  lbs  Piston  Load 
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5.  COMPUTER  PROGRAM  CONFIGURATION  DOCUMENTATION 


The  following  documentation  presents  the  traceability  for  computer  programs  used 
in  the  analysis  reported  here.  This  documentation  should  accompany  any  and  all 
analysis  reports  transmitted  by  Applied  Mechanics. 


Task:  Analysis  of  Selected  Eaton  Throat-Valve  Element  Components 

Charge  No.  or  Work  Package:  820253914 
Report  Title:  Same  as  task 

Author:  S.  D.  Snow  Date:  December  1993 

Program  Used:  I— DEAS  Version:  VI. i  Module:  NA 

Computer  Used:  WCUSV2  Model:  DECstation  5000/200 

Verification  Manual/Test  Problem  Manual/Example  Manual: 

W.  D.  Richins  letter  to  Applied  Mechanics,  ’’Verification  of  SDRC  l-DEAS  VI  Software,” 
WDR-28-91 

Author:  W.  D.  Richins  Date:  December  18, 1991 

Comments:  ”SDRC  l-DEAS  V  Verification  Manual”  is  also  available  within  the  Applied 
Mechanics  Unit. 
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6.  FINITE  ELEMENT  MODEL  INPUT  &  OUTPUT 


(The  finite  element  model  input  end  output  dete  from  the  original  release  of  this  report  did 
not  change  in  this  revision.  In  the  interest  of  conserving  pater,  that  data  was  not  included 
in  this  revision.) 
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7.  ACCELERATION  HISTORY  FOR  BACKWARD  OPENING  PISTON 
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SEAL  STUDY 


INTRODUCTION 

A  meeting  was  held  with  M.R.Stacy  (Requester)  on  Dec  7, 1993  to  develop  scope  of 
assignment,  review  problem,  and  identify  design  and  operating  conditions;  and  a  follow  up 
meeting  was  held  with  Jim  Hall  (Test  Technician)  to  determine  test  scenario,  test  conditions 
and  results. 

o  Scope  of  Assignment 

Review  Valve  operation  and  conduct  scoping  study  to  identify  available  seal  materials  and 
configurations  capable  of  withstanding  the  conditions  experienced  by  operation  of  the 
LB/TS  Throat  Valve  Element.  The  study  shall  also  identify  Manufacturers,  cost,  and  delivery. 
Schedule:  Complete  draft  study  by  December  17, 1993 
Budget:  60  mhrs 
Charge  No.  820253914 

o  Background 

EG&G  was  previously  requested  by  the  US  Army  to  evaluate  the  performance  of  the  LB/TS 
Throat  Valve  Element,  which  was  supplied  to  them  by  Eaton  Corporation.  A  series  of  tests 
were  conducted.  However,  design  conditions  were  never  achieved  because  of  seal  failure 
and  excessive  leakage.  The  seals  were  completely  destroyed  during  one  of  these  tests. 

(see  figure  1  ).  Close  examination  of  this  photo,  which  was  taken  after  the  performance  test, 
reveals  that  only  three  of  the  six  seals  were  damaged.  These  are  the  three  seals  that  cross 
the  ports  locaited  on  the  valve  body.  Since  the  other  three  seals  appear  to  be  undamaged,  it 
is  reasonable  to  assume  that  damage  was  due  to  pinching  rather  than  temperature, 
however,  high  temperature  may  be  a  significant  contributor. 

o  Test  Scenario,  Test  Conditions,  and  Results 
Three  series  of  tests  were  conducted  as  follows: 

First  test  using  original  seal  design 

Two  different  configurations  of  seals  were  used.  The  design  of  the  seals  crossing  the 
cylinder  ports  was  comprised  of  a  teflon  15%  graphite  filled  split  overlapping  joint  ring 
(piston  ring  type)  with  a  rectangular  cross  section,  a  radial  metal  spring  (approx  0.010  thk  x 
.25  wide),  and  a  metal  spacer  (approx  0.010  thk  x  .25  wide)  between  the  spring  and  the 
gland  diameter.  The  design  of  the  seals  that  were  not  required  to  cross  the  ports  was 
comprised  of  a  teflon  15%  graphite  filled  guide  ring  with  a  lip  on  the  leading  edge.  This 
guide  ring  was  supported  by  a  silicon  o-ring  with  two  segmented  half  washers  0.030  thk 
that  backed  up  the  o-ring.  See  figure  2.  Experienced  significant  leakage.  Unable  to  reach 
design  temperatures  or  pressures.  Leaked  at  approximately  400  psi  and  350  degrees  F. 
Installation  of  the  teflon  guide  ring  required  use  of  an  expanding  mandrel  type  tool  to 
expand  the  ring  over  the  piston  diameter  and  the  use  of  a  swaginging  type  tool  to  then 
reduce  the  seal  diameter  after  installation  into  the  seal  gland. 

Second  Test 

The  same  seal  design  was  used  for  all  six  seals.  This  design  was  of  the  same  configuration 
used  previously  for  those  seals  not  crossing  the  ports  (guide  ring  with  o— ring  energizer  and 
segmented  backup  rings).  The  seals  crossing  the  ports  were  destroyed.  The  lip  of  the  teflon 
guide  ring  was  sheared  off  during  valve  cycling.  See  figure  1.  Unable  to  get  any  where  near 
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design  pressure  or  temperature. 


Third  series  of  tests 

All  six  seals  were  of  the  same  configuration.  A  thicker  teflon  15%  graphite  filled  guide  ring 
with  a  rectangular  cross  section  was  backed  up  by  a  silicon  o-ring  energizer  and  the  two 
segmented  metal  half  rings.  Tests  were  conducted  at  room  temperature.  Leakage  was 
experienced  at  approximately  375  psi. 

NOTE:  Many  of  the  test  details  were  recalled  from  memory  and  there  was  some 

difficulty  in  developing  tie  actual  test  scenario.  Hence,  no  real  evidence  of 
actual  sea!  loading  exists. 

o  Valve  Operating  Conditions.  Materials,  Operation,  and  Configuration 
The  requester  was  consulted  and  drawings  were  reviewed  to  determine  Valve  operating 
conditions,  materials.  Operation,  and  configuration  as  follows: 

Valve  Operating  Conditions 

inlet  pressure  :  2000  psig 
outlet  pressure  :  0  psig 
inlet  temperature  :  700  degree  F 
outlet  temperature :  70  degree  F  (room) 
flow  media :  nitrogen  gas 
ambient  humidity  :  10-1 00% 
operating  speed  (close  to  open) :  45  ms 
piston  travel  (close  to  open) :  1.12  inches 

acceptable  leakage  (valve  closed) :  87lbm/hr  at  2000  psig  and  59  degree  F 
seal  design  life  :  150  cycles  between  replacements 

Materials 

Valve  Body  :  CRES  17-4  PH  Cond  H925 

Valve  Sleeve  :  CRES  17-4  PH  Cond  HI  025  (OD  chrome  plated  by  Electrolize  process) 

Seal :  Teflon  15%  graphite  filled  ( drawing  not  available  to  EG&G)  Seal  configuration  as 
described  above. 

Valve  Operation 

Ail  components  at  room  temperature,  valve  closed,  exterior  pressurized  to  2000  psig  with 
700  degree  F  nitrogen  gas  and  held  for  15  minutes.  Valve  is  then  opened.  Full  rest  to  full 
open  (1.12  inches)  in  45  ms.  Sleeve  oscillates  at  open  position  until  damped  out. 

Valve  Configuration 

Valve  Body  :  Cylinder  5.1 97  dia.  with  3  sets  of  16  siots  (ports)  equally  spaced  around  the 
circumference  of  the  cylinder.  Each  port  is  0.497  wide  by  0.750  long.  Each  set  of  ports  is 
spaced  2.20  inches  apart. 

Valve  Sleeve  :  Cylinder  5.1 91  OD  with  6  0-Ring  glands  and  3  set  of  slots  (ports).  Each  set 
of  ports  consists  of  8  slots  0.449  in  wide  equally  spaced  around  the  circumference  of  the 
cylinder  and  separated  by  0.100  in  wide  webs.  Each  set  of  ports  is  spaced  2.20  inches 
apart. 

ASSESSMENT  OF  PROBLEM 

The  0— Ring  seal  glands  are  located  on  the  OD  of  the  Valve  sleeve,  hence,  a  continuous 
ring  seal  must  be  stretched  over  the  OD  of  the  Sleeve  for  seal  installation.  During  valve 
operation,  the  seals  are  moved  at  high  speed  across  the  cylinder  ports.  The  present  seal  is 
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a  fairly  hard  polymer  (15%  teflon  graphite  filled)  at  room  temperature,  however,  teflon  will 
cold  flow  and  as  the  seal  passes  the  cylinder  port,  it  is  left  unsupported.  The  seal  is  then 
extruded  into  the  port  and  pinched  or  stripped. 

SURVEY  OF  MATERIALS 

A  literature  search  was  conducted  to  identify  candidate  materials  for  this  seal  application.  The 
most  commonly  used  seal  materials  for  piston/cylinder  applications  are  the  elastomers. 
However,  silicon,  which  has  the  highest  temperature  limit  for  continued  service  is  limited  to 
600  degree  F.  Table  1  provides  temperature  limits  for  these  materials 


Table  1  -Temperature  for  Continued  Service  (Elastomers) 


elastomer 

Temp  (F) 

Nitrile  rubber  NBR 

300 

Natural  Isoprene  NR 

250 

Ethylene  Propylene  EPM,EPDM 

350 

Flurosiiicon  FVMQ 

400 

Flurocarbon,  FPM 

>550 

Silicon  VMQ 

600 

Ref:  Materials  Engineering,  Material 
Selector- 1987.  pi 77- 179 _ 


Because  of  the  low  temperature  limit,  no  manufacturers  of  elastomer  seals  other  than  silicons 
and  flurocarbons  were  contacted. 

Supplier  literature  from  various  Seal  Manufactures  was  also  reviewed  to  select  candidate 
suppliers.  These  candidate  suppliers  were  then  contacted  to  determine  their  interest  and 
capability.  If  they  expressed  interest,  an  information  package  consisting  of  drawings  and  fact 
sheet  describing  materials,  operation,  and  operating  conditions  was  sent  to  them  by  FAX. 

No  manufacturers  of  silicon  or  flurocarbon  seals  have  expressed  interest  in  submitting  a 
proposal  for  this  application.  See  Appendix  (Memos  of  conversation).  They  have  expressed 
the  following  concerns: 

1  Temperature-  Seals  made  of  Teflon,  a  flurocarbon,  are  fabricated  by  the 
sintering  (powder  metallurgy)  process  at  approximately  700  degree  F.  At 
approximately  720  degree  F,  teflon  has  the  consistency  of  jelly.  [Ref  2  and  3] 

2  Ports  -  Teflon  will  cold  flow  and  with  the  high  design  pressure,  the  seals  can 
be  expected  to  expand  into  the  ports  and  be  pinched  as  they  pass  over  the 
port.  Teflon  has  a  high  temperature  coefficient  of  expansion,  hence,  it  can  be 
expectected  that  under  high  temperature,  it  will  expand  to  fill  the  gland  space 
and  when  unsupported  at  the  outer  diameter,  especially  with  high  load  ,  will 
expand  into  the  port. 

Information  provided  by  the  Flurocarbon  Mechanical  Seal  Division  indicates  that  temperature 
limits  for  their  seal  normally  rated  at  550  degree  (F)  is  approximately  325  degree  (F)  at  2000 
psi  radial  dynamic  seal  with  reciprocating  motion.  See  appendix  (Data). 

Manufacturers  contacted  include: 

Bal  Seal  Engineering  Co.  Inc. 

620  West  Warner  Ave 
Santa  Ana,  Ca  92707  3396 

Tel  (71 4) -557-51 92  Fax  (714)-241  -0185  Michael  Binder 
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Hydra  Pak  Seal  Div  (Parker  Hannafin) 

Salt  Lake  City,  Ut 

Tel  (801 ) - 973-7325  Jeff  Walker 


FURON  Mechanical  Seal  Div 
(formerly  FLUROCARBON) 

PO  box  520 

441 2  Corporate  CenterDr 
Los  Alamitos,  Ca  90720 

Tel  (71 4) -995- 181 8  Fax  (71 4) -761 -1270  Mark  Maloney 
Row,  Inc 

703-T  Annoreno  Dr 
Addison,ll.601 01 

Tel  (708) -628- 9221  Fax  (708) -628-9229  Ores  Kozinczuk 

Du  Pont,  manufacturers  of  VESPEL,  a  polyimide,  claim  that  VESPEL  will  retain  it’s  properties  at 
temperatures  from  cryogenic  to  500  degrees  F,  with  excursions  to  900  degrees  F.  Therefore, 
they  were  contacted  to  determine  if  VESPEL  could  be  used  as  a  piston/cylir-  Jer  seal  in  our 
application  and  to  determine  if  Du  Pont  would  be  interested  in  submitting  a  proposal.  Mr.  Ed 
Miller  reaffirmed  Du  Pont’s  claim  and  expressed  interest  in  submitting  a  proposal.  An 
information  package,  including  drawings  and  fact  sheet  were  sent  by  Fax  to  Mr.  Miller. 

Du  Pont  Co. 

Tel  (302)-733-8120  Fax  (302) -733-  8362  Ed  Miller 

Temperature  limitations  for  various  metals  and  coatings  are  provided  in  tables  2  and  3.  Piston 
rings  fabricated  from  cast  iron,  commonly  used  in  the  hydraulic  and  automotive  industries  at 
temperatures  of  up  to  850  degrees  are  not  recommended  for  high  pressure  nitrogen  service 
above  450  degree  F.  Severe  leakage  would  be  expected  until  the  rings  wear  in. 


Table  2  -Metal  0-Ring  Seals 

matl 

max  temp  (C) 

max  temp  (F) 

copper 

400 

752 

mild  steel 

550 

1022 

cupra- nickel 

600 

1112 

monel 

600 

1112 

nickel 

700 

1292 

stainless  steel 

800 

1472 

inconel 

850 

1562 

RefrU 
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Table  3  - 

Plating  materials  for  metal  0 

-Ring  seals 

matl 

max  temp  (C) 

max  temp  (F) 

Indium 

140 

284 

cadmium 

200 

392 

silver 

800 

1472 

gold 

850 

1562 

PTFE 

300 

572 

Various  suppliers  of  metal  piston/seal  rings  were  contacted  to  determine  their  capabilities  and 
interest.  Information  packages  were  sent  via  Fax  to  those  expressing  interest.  These  include: 

Precision  Ring  Inc. 

5611  Progress  Road 
PO  box  418187 
Indianapolis,  In.  46241 

Tel  (31 7) -247-4786  Fax  (31 7) -248-  9781  Dave  Dyer 

C.Lee  Cook/Dover  Corp 
PO  box  1038 

Louisville  KY  40201-1038 

Tel  (5020-587-6783  Fax  (502)587-6786  Paul  Hanlon 

Grover  Piston  Ring  Co. 

PO  box  340080 
Milwaukee.  Wl  53234-0080 

Tel  (41 4) -384- 9472  Fax  (41 4) -384- 0201  Paul  Cero  Mike  Knoebel 

Superior  Piston  Ring  Co. 

6427  Epwort  Blvd 
Detroit,  Ml  48210 

Tel  1-800-  899-1827  Fax  (313)-361  -0530 

A  quotation  has  been  received  from  C.Lee  Cook/Dover.  They  quoted  $1 15.02  each  for 
quantiteis  of  6-1 1  and  $78.21  each  for  quanties  of  12-49  with  delivery  of  1 0+  weeks.  Their 
quotation  did  not  include  costs  for  any  special  quality  control  specifications.  Any  efforts 
required  to  meet  these  would  be  billed  at  cost. 

SEAL  CONFIGURATIONS 

Seal  rings  are  provided  in  various  configuration,  however,  piston  ring  seals  where  the  seal 
gland  is  machined  on  the  piston  require  an  expandable  seal  for  installation  over  the  piston 
diameter.  Elastomer  seals,  which  are  expandable,  are  provided  in  a  continous  ring 
confiuration  while  metal  seals  must  have  a  split  joint  for  this  type  of  application.  Flurocarbons 
(Teflon)  are  expandable  to  a  small  extent  and  with  tooling,  thin  teflon  rings  can  be  stretched 
over  the  piston  and  swaged  back  to  original  diameter  after  installed  in  the  gland.  However, 
thicker  rings  must  be  split  for  installation.  Elastomeric  seal  can  be  provided  with  various  cross 
sections  including  round,  square,  X,  T,  V,  rectangular,  and  lip  configurations.  They  are  also 
sometimes  used  with  back  up  face  rings  to  prevent  extrusion  or  blow  out  and  with  spring 
retainers  or  energizers.  Flurocarbon  seals  can  also  be  provided  with  a  teflon  outer  layer  and  a 
silicon  or  viton  core.  Metal  seals  can  also  be  provided  with  various  cross  sections  including  O, 
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X  T,  K,  and  C  onfigurations,  however,  continuous  metal  rings  are  generally  used  as  face  seals 
or  where  they  can  be  installed  without  expansion.  They  are  are  also  sometimes  used  with  back 
up  rings  and  spring  energizers  or  pressurized  in  some  cases.  However,  a  metal  seal  for  this 
application  would  have  to  have  a  split  joint.  Joint  configuration  for  Piston/Seal  rings  include  45 
degree  split  joint,  stepped,  straight  cut,  stepped  tongue,  angled  tongue,  and  stepped  hook 
joint.  They  are  also  sometimes  used  in  combination.  See  Appendix  (Data-  Design  Update  - 
Dover  Corporation). 

CONCLUSIONS 

1  If  Testing  would  have  continued,  there  is  fair  probability  that  both  the  teflon  guide  ring  and 
the  silicon  seal  energizer  ring  would  have  failed  since  their  operating  temperature  limit  is 
600  degree  F. 

2  Seal  failure  can  be  attributed  to  pinching  of  the  seal  while  passing  over  the  ports 
due  to  unsupported  section  of  a  material  that  will  cold  flow  under  high  load 
conditions. 

RECOMMENDATIONS 

1  Consideration  should  be  given  to  performance  of  a  temperature  test  to  determine  behavior 
of  teflon  (15%  graphite  filled)  at  700  degree  F. 

2  Conduct  Seal  Test  Program  to  develop  seals  capable  of  meeting  design  conditions. 

Program  should  include  testing  of  split  metal  piston  ring  type  seals  and  VESPEL  seals. 

3  Continue  efforts  to  obtain  proposals  from  fabricators  of  split  metal  seals. 

4  Continue  to  follow  up  DuPont  on  proposal  of  VESPEL  seal. 
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APPENDIX  H: 

PROPOSED  CONFIGURATION  DRAWINGS 

DRAWINGS 

MEMOS  OF  CONVERSATION 

MATERIAL  DATA/SUPPLIER  LITERATURE 

CALCULATIONS 


This  appendix  has  been  sent  to  R.  J.  Pearson  file  at  the  U.S.  Army  Research  Laboratory, 
Aberdeen  Proving  Ground,  Maryland. 
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CONTENTS 


•  2  Bernoulli  disks  containing  data  outputs  from  the  model  for  the 
different  runs: 

Disk  1  contains  runs  at  1 400  psi  driver  gas  pressure 
Disk  2  contains  runs  at  1 800  psi  driver  gas  pressure. 

•  1  Bernoulli  disk  in  the  diskette  mailer  envelope,  containing  test 
data  from  the  1991  tests 

•  2  floppy  diskettes  with  ASCII  files  of  the  data  obtained  from 
testing  the  ETVE  with  the  LVDT  Installed  during  1993 

•  The  LVDT  and  associated  instrumentation 

•  Updated  sheets  for  the  camera-ready  report 

•  Relevant  information  for  future  reference  regarding  ETVE 
hardware  (seals,  seals  configurations  used,  Belleville  washers, 
solenoid  valves,  etc.) 

•  ETVE  as  built  drawings 

•  ETVE  as  built  drawings  with  LVDT  attached 

•  Final  project  status  letter 
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3  PHILLIPS  LABORATORY  (AFWL) 

ATTNNTE 

NTED 

NTES 

KIRTLAND  AFB  NM  87117-6008 


1  AFTT 

ATTN  TECHNICAL  LIBRARY 
BLDG  640  B 

WRIGHT  PATTERSON  AFB  OH  45433 

1  FTDNIIS 

WRIGHT  PATTERSON  AFB  OH  45433 
4  DIR 

ATTN  R  GUENZLER  MS  3505 
R  HOLMAN  MS-3510 
R  A  BERRY 
WCREED 

IDAHO  NATIONAL  ENGINEERING  LABORATORY 
EG&G  IDAHO  INC 
PO  BOX  8757  BWI  AIRPORT 
BALTIMORE  MD  21240 

3  KAMAN  SCIENCES  CORPORATION 
ATTN  LIBRARY 
PA  ELLIS 
F  H  SHELTON 
P  O  BOX  7463 

COLORADO  SPRINGS  CO  80933-7463 

2  DIR 

ATTN  TH  DOWLER  MS  F602 

DOC  CONTROL  FOR  REPORTS  LIBRARY 

LOS  ALAMOS  NATIONAL  LABORATORY 

PO  BOX  1663 

LOS  ALAMOS  NM  87545 

1  DIR 

ATTN  DOC  CONTROL  FOR  TECH  LIB 
SANDIA  NATIONAL  LABORATORIES 
LIVERMORE  LABORATORY 
P  O  BOX  %9 
UVERMORE  CA  94550 

1  DIR 

ATTN  TECHNICAL  LIBRARY 
NASA  LANGLEY  RESEARCH  CENTER 
HAMPTON  VA  23665 

1  ADA  TECHNOLOGIES  INC 

ATTN  JAMES  R  BUTZ 
HONEYWELL  CENTER  SUITE  110 
304  INVERNESS  WAY  SOUTH 
ENGLEWOOD  CO  80112 
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1  ALLIANT  TECHSYSTEMS  INC 

ATTN  ROGER  A  RAUSCH  MN48  3700 
7225  NORTHLAND  DRIVE 
BROOKLYN  PARK  MN  55428 

1  AEROSPACE  CORPORATION 

ATTN  TECH  INFO  SERVICES 
P  O  BOX  92957 
LOS  ANGELES  CA  90009 

1  THE  BOEING  COMPANY 

ATTN  AEROSPACE  LIBRARY 
P  O  BOX  3707 
SEATTLE  WA  98124 

1  CALIFORNIA  RES  &  TECH  INC 

ATTN  M  ROSENBLATT 
20943  DEVONSHIRE  STREET 
CHATSWORTH  CA  91311 

1  DYNAMICS  TECHNOLOGY  INC 

ATTN  D  T  HOVE 
GP  MASON 

21311  HAWTHORNE  BLVD  SUITE  300 
TORRANCE  CA  90503 

1  EATON  CORPORATION 

ATTN  J  WADA 

DEFENSE  VALVE  &  ACTUATOR  DIV 
2338  ALASKA  AVE 
EL  SEGUNDO  CA  90245-4896 

5  DIR 

ATTN  DOC  CONTROL  3141 
C  CAMERON  DIV  6215 
A  CHABAI  DIV  7112 
D  GARDNER  DIV  1421 
J  MCGLAUN  DIV  1541 
SANDIA  NATIONAL  LABORATORIES 
P  O  BOX  5800 

ALBUQUERQUE  NM  87185-5800 

1  BLACK  &  VEATCH 

ENGINEERS  -  ARCHTTECTS 
ATTN  H  D  LAVERENTZ 
1500  MEADOW  LAKE  PARKWAY 
KANSAS  CITY  MO  641 14 


1  DIRECTOR 

ATTN  DR  T  HOLTZ  MS  202-14 
NASA-AMES  RESEARCH  CENTER 
APPLIED  COMPUTATIONAL  AERO  BRANCH 
MOFFETT  FIELD  CA  94035 

2  APPLIED  RESEARCH  ASSOCIATES  INC 
ATTN  J  KEEFER 

NH  ETHRIDGE 
P  O  BOX  548 
ABERDEEN  MD  21001 

1  CARPENTER  RESEARCH  CORPORATION 

ATTN  H  JERRY  CARPENTER 
27520  HAWTHORNE  BLVD  SUTTE  263 
ROLLING  HILLS  ESTATES  CA  90274 

1  GOODYEAR  CORPORATION 
ATTN  R  M  BROWN  BLDG  1 
SHELTER  ENGINEERING 
UTCHFEELD  PARK  AZ  85340 

2  FMC  CORPORATION 
ATTN  J  DROTLEFF 
C  KREBS  MDP95 
ADVANCED  SYSTEMS  CENTER 
BOX  58123 

2890  DE  LA  CRUZ  BLVD 
SANTA  CLARA  CA  95052 

1  SVERDRUP  TECHNOLOGY  INC 

ATTN  B  D  HEKKINEN 
SVERDRUP  CORPORATION  AEDC 
MS  900 

ARNOLD  AFB  TN  37389-9998 

1  KTECH  CORPORATION 

ATTN  DR  E  GAFFNEY 
901  PENNSYLVANIA  AVE  NE 
ALBUQUERQUE  NM  87111 

4  KAMAN  AVTDYNE 

ATTN  R  RUETENDC  (2  CYS) 

S  CRISCIONE 
R  MILLIGAN 
83  SECOND  AVENUE 
NORTHWEST  INDUSTRIAL  PARK 
BURLINGTON  MA  01830 
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2  KAMAN  SCIENCES  CORPORATION 
ATTN  DASIAC  (2  CYS) 

P  O  DRAWER  1479 

816  STATE  STREET 

SANTA  BARBARA  CA  93102-1479 

1  LOCKHEED  MISSILES  &  SPACE  CO 

ATTN  J  J  MURPHY 
DEPT  81  11  BLDG  154 
P  O  BOX  504 
SUNNYVALE  CA  94086 

1  ORLANDO  TECHNOLOGY  INC 
ATTN  D  MATUSKA 

60  SECOND  STREET  BLDG  5 
SHALIMAR  FL  32579 

2  THE  RALPH  M  PARSONS  COMPANY 
ATTN  T  M  JACKSON 

LB  TS  PROJECT  MANAGER 
100  WEST  WALNUT  STREET 
PASADENA  CA  91124 

1  SAIC 

ATTNNSINHA 

501  OFFICE  CENTER  DRIVE  APT  420 
FT  WASHINGTON  PA  19034  3211 

1  SAIC 

ATTN  J  GUEST 
2301  YALE  BLVD  SE 
SUTTEE 

ALBUQUERQUE  NM  87106 

2  S  CUBED 

A  DIVISION  OF  MAXWELL  LABS  INC 
ATTN  C  E  NEEDHAM 
L KENNEDY 
2501  YALE  BLVD  SE 
ALBUQUERQUE  NM  87106 

1  TRW 

BALLISTIC  MISSILE  DIVISION 
ATTN  H  KORMAN 
MAIL  STATION  526  614 
P  OBOX  1310 
SAN  BERNADINO  CA  92402 

1  THERMAL  SCIENCE  INC 

ATTN  R  FELDMAN 
2200  CASSENS  DR 
ST  LOUIS  MO  63026 


2  MCDONNELL  DOUGLAS  ASTRNTCS  CORP 

ATTN  ROBERT  W  HALPRIN 
K  AHEINLY 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH  CA  92647 

1  MDA  ENGINEERING  INC 
ATTN  DR  DALE  ANDERSON 
500  EAST  BORDER  STREET 
SUITE  401 

ARLINGTON  TX  07601 

2  PHYSICS  INTERNATIONAL  CORPORATION 
P  O  BOX  5010 

SAN  LEANDRO  CA  94577-0599 

1  R&D  ASSOCIATES 

ATTN  G  P  GANONG 
P  O  BOX  9377 
ALBUQUERQUE  NM  87119 

1  SCIENCE  CENTER 

ROCKWELL  INTERNATIONAL  CORPORATION 
ATTN  DR  S  CHAKRAVARTHY 
DR  D  OTA 

1049  CAMINO  DOS  RIOS 
THOUSAND  OAKS  CA  91358 

3  S  CUBED 

A  DIVISION  OF  MAXWELL  LABS  INC 

ATTN  TECHNICAL  LIBRARY 

RDUFF 

KPYATT 

P  O  BOX  1620 

LA  JOLLA  CA  92037-1620 

1  SUNBURST  RECOVERY  INC 

ATTN  DR  C  YOUNG 
P  O  BOX  2129 

STEAMBOAT  SPRINGS  CO  80477 

1  SVERDRUP  TECHNOLOGY  INC 

ATTN  R  F  STARR 
P  O  BOX  884 
TULLAHOMA  TN  37388 

1  SRI  INTERNATIONAL 

ATTN  DR  G  R  ABRAHAMSON 
DR  J  GRAN 
DR  B  HOLMES 
333  RAVENWOOD  AVENUE 
MENLO  PARK  CA  94025 


DIST-6 


NO.  OF 

COPffiS  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  BATTELLE 
TWSUAC 

505  KING  AVENUE 
COLUMBUS  OH  43201-2693 

2  THINKING  MACHINES  CORPORATION 
ATTN  G  SABOT 

RFERREL 

245  FIRST  STREET 

CAMBRIDGE  MA  02142-1264 

1  CALIFORNIA  INSITTUTE  OF  TECHNOLOGY 

ATTN  T  J  AHRENS 
1201  E  CALIFORNIA  BLVD 
PASADENA  CA  91109 

1  UNIVERSITY  OF  MINNESOTA 
ARMY  OF  HIGH  PERF  COMP  RES  CTR 
ATTN  DR  TAYFUN  E  TEZDUYAR 
1100  WASHINGTON  AVE  SOUTH 
MINNEAPOUS  MN  55415 

2  CDR 

ATTN  SSCNC  YSD  J  ROACH 
SSCNC  WST  A  MURPHY 
US  ARMY  NRDEC 
KANSAS  STREET 
NATICK  MA  10760-5018 

3  SOUTHWEST  RESEARCH  INSHTUTE 
ATTN  DR  C  ANDERSON 
SMULLIN 

AB  WENZEL 

P  O  DRAWER  28255 

SAN  ANTONIO  TX  78228-0255 

1  STATE  UNIVERSITY  OF  NEW  YORK 
MECHANICAL  &  AEROSPACE  ENGINEERING 
ATTN  DR  PEYMAN  GTVI 

BUFFALO  NY  14260 

2  DENVER  RESEARCH  INSTITUTE 
ATTN  J  WISOTSKI 
TECHNICAL  LIBRARY 

P  O  BOX  10758 
DENVER  CO  80210 

2  UNIVERSITY  OF  MARYLAND 

INSTITUTE  FOR  ADV  COMPUTER  STUDIES 
ATTN  L  DAVIS 
G  SOBIESKI 

COLLEGE  PARK  MD  20742 


1  NORTHROP  UNIVERSITY 

ATTN  DR  F  B  SAFFORD 
5800  W  ARBOR  VITAE  STREET 
LOS  ANGELES  CA  90045 

1  STANFORD  UNIVERSITY 

ATTN  DR  D  BERSHADER 
DURAND  LABORATORY 
STANFORD  CA  94305 


ABERDEEN  PROVING  GROUND 

1  CDR  USATACOM 

ATTN  AMSTE-TE-F  (L  TELETSKl) 

1  CDR  USATHAMA 

ATTN  AMSTH-TE 

1  CDR  USACSTA 

ATTN  STEC-U 

5  DIR  USARL 

ATTN  AMSRL-SL-CM  (BLDG  E3331) 

AMSRL-WT-NC  (R  LOTTERO)(2  CP) 
AMSRL-WT-NC  (A.  MIHALCIN)(2  CP) 


Intentionally  left  blank. 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-CR-239 _ Date  of  Report  July  1995 _ 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


DEPARTMENT  OF  THE  ARMY 


OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001,APG,MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


DIRECTOR 

U.S.  ARMY  RESEARCH  LABORATORY 
ATTN:  AMSRL-WT-NC 

ABERDEEN  PROVING  GROUND,  MD  21005-5066 


NO  POSTAGE 
NECESSARY 
IF  MAILED 
IN  THE 

UNITED  STATES 


